For more than five decades, hydraulic fracturing has been widely used to enhance oil and gas production. Hydraulic fracturing in solid materials (e.g., rock) has been studied extensively. The main goal of this thesis is a comprehensive study of the physical mechanisms of hydraulic fracturing in cohesionless sediments. For this purpose, experimental techniques are developed to quantify the initiation and propagation of hydraulic fractures in dry particulate materials. We have conducted a comprehensive experimental series by varying such controlling parameters as the properties of particulate materials and fracturing fluids, boundary conditions, initial stress states, and injection volumes and rates. In this work, we suggest principle fundamental mechanisms of hydraulic fracturing in particulate materials and determine relevant scaling relationships (e.g., the interplay between elastic and plastic processes). The main conclusion of this work is that hydraulic fracturing in particulate materials is not only possible, but even probable if the fluid leak-off is minimized (e.g., high flow rate, high viscosity, low permeability). Another important conclusion of this work is that all parts of the particulate material are likely to be in compression. Also, the scale effect (within the range of the laboratory scales) appears to be relatively insignificant, that is, the observed features of fractures of different sizes are similar. Based on the observed fracture geometries, and injection pressures we suggested three models of hydraulic fracturing in particulate materials. In the cavity expansion or ??e driving model, the fracturing fluid is viewed as a sheet pile (blade) that disjoints the host material, and the cavity expansion occurs at the fracture (blade) front. The shear banding model is also consistent with a compressive stress state everywhere in the particulate material and explains the commonly observed beveled fracture front. The model of induced cohesion is based on the fluid leak-off ahead of the fracture front. The induced cohesion may be caused by the tensile strain near the fracture tip (where the stress state is also compressive), which, in turn, induces the cavitation of the leaked-off fluid and hence capillary forces.

Multistage hydraulic fracturing of horizontal wells continues to be a major technological tool in the oil and gas industry. Creation of multiple transverse fractures in shale gas has enabled production from very low permeability. The strategy entails the development of a Stimulated Reservoir Volume (SRV), defined as the volume of reservoir, which is effectively stimulated to increase the well performance. An ideal model for a shale gas SRV is a rectangle of length equal to horizontal well length and width equal to twice the half length of the created hydraulic fractures. This project focused on using the Multistage Transverse Fractured Horizontal Wells (MTFHW) for two novel applications. The first application considers using the SRV of a shale gas well, after the gas production rate drops below the economic limit, for low grade geothermal heat extraction. Cold water is pumped into the fracture network through one horizontal well drilled at the fracture tips. Heat is transferred to the water through the fracture surface. The hot water is then recovered through a second horizontal well drilled at the other end of the fracture network. The basis of this concept is to use the already created stimulated reservoir volume for heat transfer purposes. This technique was applied to the SRV of Haynesville Shale and the results were discussed in light of the economics of the project. For the second application, we considered the use of a similarly created SRV for producing hydrocarbon products from oil shale. Thermal decomposition of kerogen to oil and gas requires heating the oil shale to 700 degrees F. High quality saturated steam generated using a small scale nuclear plant was used for heating the formation to the necessary temperature. Analytical and numerical models are developed for modeling heat transfer in a single fracture unit of MTFHW. These models suggest that successful reuse of Haynesville Shale gas production wells for low grade geothermal heat extraction and the project appears feasible both technically and economically. The economics of the project is greatly aided by eliminating well drilling and completion costs. The models also demonstrate the success of using MTFHW array for heating oil shale using SMR technology.

Tight gas and shale oil play an important role in energy security and in meeting an increasing energy demand. Hydraulic fracturing is a widely used technology for recovering these resources. The design and evaluation of hydraulic fracture operation is critical for efficient production from tight gas and shale plays. The efficiency of fracturing jobs depends on the interaction between hydraulic (induced) and naturally occurring discrete fractures. In this work, a coupled reservoir-fracture flow model is described which accounts for varying reservoir geometries and complexities including non-planar fractures. Different flow models such as Darcy flow and Reynold's lubrication equation for fractures and reservoir, respectively are utilized to capture flow physics accurately. Furthermore, the geomechanics effects have been included by considering a multiphase Biot's model. An accurate modeling of solid deformations necessitates a better estimation of fluid pressure inside the fracture. The fractures and reservoir are modeled explicitly allowing accurate representation of contrasting physical descriptions associated with each of the two. The approach presented here is in contrast with existing averaging approaches such as dual and discrete-dual porosity models where the effects of fractures are averaged out. A fracture connected to an injection well shows significant width variations as compared to natural fractures where these changes are negligible. The capillary pressure contrast between the fracture and the reservoir is accounted for by utilizing different capillary pressure curves for the two features. Additionally, a quantitative assessment of hydraulic fracturing jobs relies upon accurate predictions of fracture growth during slick water injection for single and multistage fracturing scenarios. It is also important to consistently model the underlying physical processes from hydraulic fracturing to long-term production. A recently introduced thermodynamically consistent phase-field approach for pressurized fractures in porous medium is utilized which captures several characteristic features of crack propagation such as joining, branching and non-planar propagation in heterogeneous porous media. The phase-field approach captures both the fracture-width evolution and the fracture-length propagation. In this work, the phase-field fracture propagation model is briefly discussed followed by a technique for coupling this to a fractured poroelastic reservoir simulator. We also present a general compositional formulation using multipoint flux mixed finite element (MFMFE) method on general hexahedral grids with a future prospect of treating energized fractures. The mixed finite element framework allows for local mass conservation, accurate flux approximation and a more general treatment of boundary conditions. The multipoint flux inherent in MFMFE scheme allows the usage of a full permeability tensor. An accurate treatment of diffusive/dispersive fluxes owing to additional velocity degrees of freedom is also presented. The applications areas of interest include gas flooding, CO₂ sequestration, contaminant removal and groundwater remediation. / text

This Master’s thesis seeks to develop a better understanding of how Indigenous voices can be included in water governance. As a starting point, social learning theory, collaborative governance and Indigenous and Canadian relations were carefully studied. Despite the large body of research on collaborative governance with First Nations and on social learning in water governance, little is known specifically about social learning processes in colonial contexts. Using grounded theory and Indigenous methodologies, this research investigates how the current approach to implementing the Province’s legal constitutional, “duty to consult” affects social learning processes and the inclusion of Indigenous people in water governance. Findings indicate that the laws and policies that have been created based on the Crown’s interpretation of Treaty 8, an agreement signed between the Fort Nelson First Nation and Canada in 1899. This duty to consult constrains social learning, as it does not allow for the flexibility needed for a reframing process that might bring the actors to a common understanding of Treaty 8, the treaty relationship and its application as such today, as a basis for future collaboration. Without reframing processes, the consultation process is perceived by the Fort Nelson First Nation, a Treaty 8 nation, to lack legitimacy and neutral facilitation. Subsequently, consultation is seen as a checkbox that must be completed, but fails to include First Nations’ knowledge, interests and concerns about impacts from development and appropriate accommodation. This research also investigates a new governance arrangement emerging in northeast BC, which changes the way Fort Nelson First Nation voices are included in decision-making. Processes of nation building and capacity building contribute to Fort Nelson First Nation’s New Approach to governance. The New Approach sees changes to the sites of authority, revenue and norms and beliefs, resulting in a governance innovation that circumvents the provincial government’s role in governance by creating a closer working relationship between industry and the Fort Nelson First Nation. The results are development planning and decisions that better reflect the Fort Nelson First Nation’s concerns and interests in the near future. / Graduate

L'objectif de cette thèse est de mettre en place un cadre numérique qui permet de simuler la réponse de structures en béton armé de taille industrielle soumises à des impacts sévères. Notre modèle repose sur trois ingrédients principaux : 1) la modélisation éléments discrets (ED) du béton, permettant de profiter de la nature « discrète » de cette formulation et de décrire facilement l'apparition et la propagation des discontinuités fortes de la matière ; 2) la modélisation éléments finis (EF) poutre des armatures, donnant la possibilité de représenter toute la complexité du ferraillage que l'on rencontre dans les structures industrielles en béton armé ; 3) un modèle de liaison acier-béton original, que nous avons proposé et mis en œuvre dans le code EUROPLEXUS, et qui constitue le principal apport de cette thèse. Nous avons effectué une étude théorique et numérique sur la stabilité et la précision de ce modèle de liaison, et avons calibré ses paramètres en simulant un essai d'arrachement. En modélisant l'essai de traction d'un tirant (barre en béton contenant une armature), nous avons montré la capacité de notre modèle à reproduire le transfert des efforts entre le béton et l'armature au niveau de l'interface acier-béton. Testé et validé sur des cas-tests élémentaires, le modèle mixte ED-EF a ensuite été appliqué à la simulation de vraies structures en béton armé. Nous avons simulé de manière détaillée quelques essais connus (l'impact mou sur poutre CEA, l'essai Meppen II-4), ce qui a permis de valider l'ensemble de notre approche numérique et de recueillir des éléments pour définir les pistes d'amélioration des modèles que nous avons mis en œuvre. / The aim of this work is to set up a numerical framework to simulate the behaviour of industrial size reinforced concrete structures subjected to severe impacts. Our model is based on three main features : 1) modeling of the concrete with a discrete method to handle easily strong material discontinuities such as initiation and propagation of macro-cracks ; 2) modeling of the reinforcement bars with finite element method to be able to represent complex reinforcement cages of industrial structures ; 3) an original steel-concrete bond model that we proposed and implemented in the dynamic explicit code EUROPLEXUS and that constitutes the main task of this work. We studied theoretically and numerically this model to guarantee its stability and precision during the time integration. We calibrated its parameters by simulating pull-out tests. To verify our model, we simulated a tie-test (a long concrete column which contains a reinforcement bar) and showed that our model ensures the correct transfert of forces between steel and concrete. After testing our model on simple benchmarks, we simulated real reinforced concrete structures subjected to impacts (soft impact on a beam, Meppen test n°II-4), allowing us to validate our numerical approach and to define some perspectives to improve the models we have developped.

Diss. Linköping : Univ., 2003.

Η ρύπανση του εδάφους από βιομηχανικά και αστικά απόβλητα αποτελεί παγκοσμίως ένα από τα σημαντικότερα περιβαλλοντικά προβλήματα. Ανάμεσα στους πιο επικίνδυνους ρύπους συγκαταλέγονται και οι υδρογονάνθρακες πετρελαιοειδών. Αναπόφευκτη συνέπεια της διαφυγής των υγρών ρύπων στο έδαφος είναι η μεταφορά τους προς στην ακόρεστη ζώνη του εδάφους μέσω μιας σειράς διεργασιών (ροή, διαλυτοποίηση, εξάτμιση, διασπορά, ρόφηση, κτλ) και η ρύπανση των υποκείμενων ταμιευτήρων υπογείων υδάτων. Οι πλέον κλασσικές μέθοδοι απομακρυσμένης (ex situ) απορρύπανσης εδαφών περιλαμβάνουν εκσκαφή και μεταφορά του εδάφους σε κατάλληλους χώρους (π.χ. αποτεφρωτήρες, βιο-σωροί, κλπ). Τα τελευταία χρόνια έχει αναπτυχθεί έντονο ενδιαφέρον για την ανάπτυξη και εφαρμογή τεχνολογιών επιτόπιας απορρύπανσης των εδαφών (in situ soil remediation) με χαμηλό κόστος και ελάχιστο περιβαλλοντικό αποτύπωμα. Στην παρούσα εργασία μελετώνται δύο σχετικά νέες μέθοδοι επιτόπιας απορρύπανσης: α) ο βιοαερισμός και, β) η έγχυση ατμού. Και στις δύο περιπτώσεις, για να αυξηθεί η ακτίνα δράσης της απορρύπανσης, το έδαφος διεγέρθηκε με την δημιουργία οριζόντιων υδραυλικών ρωγμών που λειτουργούν ως οριζόντια φρεάτια διαβίβασης και εξαγωγής ρευστών. Το πεδίο μελέτης είναι ένα πρώην στρατιωτικό αεροδρόμιο της Βόρειο-Δυτικής Πολωνίας, το Kluczewo, το οποίο έχει ρυπανθεί εκτεταμένα με κηροζίνη για μεγάλο χρονικό διάστημα (1935-1992). Το έδαφος της περιοχής αυτής παρουσιάζει ρωγμές ενώ η πορώδης μήτρα έχει σχετικά χαμηλή διαπερατότητα. Οι ίδιες γεωλογικές συνθήκες επικρατούν σε μεγάλο μέρος του υπεδάφους της Βόρειας Ευρώπης. Λόγω της ισχυρά ετερογενούς φύσης αυτών των εδαφών από την κλίμακα των πόρων στην κλίμακα του πεδίου και της δημιουργίας προτιμητέων μονοπατιών ροής, είναι αρκετά δύσκολο να σχεδιαστούν αποδοτικές μέθοδοι απορρύπανσης. Κύριος στόχος της εργασίας είναι η αξιολόγηση της απόδοσης των δύο μεθόδων απορρύπανσης μετά την εφαρμογή τους στην ακόρεστη ζώνη ετερογενούς εδάφους που έχει ρυπανθεί εκτενώς με υδρογονάνθρακες πετρελαιοειδών (κηροζίνη). Για να προσδιοριστεί η απόδοση κάθε μεθόδου, συλλέχθηκαν δείγματα εδάφους από ένα μεγάλο αριθμό σημείων και πραγματοποιήθηκαν χημικές αναλύσεις μέτρησης της συγκέντρωσης και της σύστασης των υδρογονανθράκων με GC-MS και GC-FID. Προκειμένου να διευκρινιστoύν καλύτερα οι κύριοι μηχανισμοί απομάκρυνσης του ρύπου και να εκτιμηθεί η αποδοτικότητα της έγχυσης ατμού και του βιοαερισμού, διεξήχθησαν πειράματα έγχυσης ατμού και βιοαερισμού και σε εργαστηριακή κλίμακα (oρθογώνιο κελί από PMMA με διαστάσεις 55 cm x 50 cm x 12 cm και υπό ακόρεστες συνθήκες). Λόγω των ετερογενειών της πορώδους δομής, στον βιοαερισμό ο κύριος μηχανισμός απομάκρυνσης ρύπου ήταν η εξάτμιση των υδρογονανθράκων και η σχετικά γρήγορη μεταφορά των ατμών μέσω διάχυσης και λόγω της μεγάλης βαθμίδας συγκέντρωσης από την μικροπορώδη μήτρα προς μονοπάτια προτιμητέας ροής (preferential flow paths) αέρα (αερισμός). Αντιστοίχως, λόγω της χαμηλής διαπερατότητας του εδάφους, στην περίπτωση της έγχυσης ατμού ο κύριος μηχανισμός απομάκρυνσης των ημι-πτητικών και μη πτητικών συστατικών ήταν η απόσταξη ατμού (δηλαδή η μείωση του σημείου ζέσεως των υδρογονανθράκων λόγω της παρουσίας μη αναμίξιμης υδατικής φάσης). Η σύγκριση των δύο μεθόδων απορρύπανσης με διάνοιξη οριζόντιων υδραυλικών ρωγμών έδειξε ότι και οι δύο μέθοδοι είναι αρκετά αποδοτικές όσον αφορά στη μείωση της μάζας του ρύπου (~72%) σε ισχυρά ετερογενή εδάφη με ρωγμές. Η έγχυση ατμού όμως είναι πολύ πιο γρήγορη (3 μήνες) από το βιοαερισμό (12 μήνες) αλλά ταυτόχρονα και αρκετά πιο ακριβή από αυτόν (14-25%). Όσον αφορά το ποσοστό μείωσης επικινδυνότητας του υπολειπόμενου ρύπου, ο βιοαερισμός δίνει με μεγάλη διαφορά καλύτερα αποτελέσματα (93%) από ότι η έγχυση ατμού (8-68%) όπου παρατηρείται και μεγάλη διακύμανση τιμών. / Soil contamination by industrial and urban wastes is nowadays one of the most important pollution problems worldwide. The petroleum hydrocarbons are included in the list of toxic pollutants that have contaminated extensive areas all over the world. An inevitable consequence of liquid pollutants release on the ground is their transport by a variety of mechanisms (e.g. gravity flow, dissolution, volatilization, dispersion, sorption, etc) in the unsaturated zone of soil, and subsequent pollution of the underlying aquifers. Conventional methods of ex-situ soil remediation are the soil excavation and its transportation in appropriate places (e.g. incinerator, bio-piles, etc). During the last years, there is a growing interest for the development and application of low-cost and sustainable (low environmental impact) in-situ soil remediation technologies. In the present work, two relatively new in-situ remediation technologies, bioventing and steam injection, were tested on the vadose zone of a low permeability and fractured glacial till sediment that was contaminated by jet fuel. The experimental site is situated in an abandoned military airport (Kluczewo) in North-Western Poland. The area was polluted extensively by jet fuel over a long period (1935-1992). The geological characterization revealed the existence of vertical desiccation fractures at the upper layers and horizontal/sub-horizontal tectonic fractures at the deeper ones, and a relatively low permeability and heterogeneous micro-porous matrix. The same geologic conditions dominate in a major part of the subsurface in Northern and Central Europe. Due to the multi-scale heterogeneities, ranging from the pore-scale to the field-scale, and the creation of preferential flow paths in such soils, it is very difficult to design successful remediation strategies based on vertical wells. For this reason, in both field experiments (bioventing and steam injection), the soil was stimulated by opening hydraulic fractures which acted as horizontal wells of fluid injection/extraction, and enhanced the influence radius of remediation. The main goal of the work is to evaluate, under field conditions, steam injection and bioventing as sustainable and efficient technologies for the removal of petroleum hydrocarbons from highly heterogeneous soils. In order to determine the efficiency of each remediation method, soil samples were collected from twelve wells and seven depths, and placed inside specific flasks pre-filled with dichloromethane (DCM). After the accelerated extraction of non-aqueous phase liquid (NAPL) from the soil and its dissolution in DCM, the composition and concentration of hydrocarbons (NAPL) was performed by using GC-MS and GC-FID. In order to clarify the main NAPL removal mechanisms and evaluate the effectiveness of steam injection and bioventing, lab-scale experiments were also conducted in soil tanks by using synthetic NAPL (PMMA cell with 55 cm x 50 cm x 12 cm dimensions) and keep the conditions comparable to the field ones. Due to pore structure heterogeneities, the main NAPL removal mechanism in bioventing was the volatilization of hydrocarbons and the fast vapour transfer through diffusion (due to the high concentration gradient) from the porous matrix to the high hydraulic conductivity preferential flow paths of injected air (ventilation). Respectively, due to the low soil permeability, in steam injection the main removal mechanism of semi- and non-volatile substances was steam distillation (namely the reduction of NAPL compounds bubble point because of the coexistence of immiscible water). The comparison of the two remediation technologies, steam injection and bioventing, indicate that both methods are efficient with respect to the reduction of the pollutant mass (~72%) in highly heterogeneous and fractured soils. With respect to cost issues, steam injection albeit faster, is 14-25% more expensive than bioventing. Finally, concerning the risks associated with the reduction of residual NAPL in groundwater, bioventing was more effective (~93%) than steam injection (8-68%).

[pt] O desenvolvimento e exploração das formações shale gas a nível mundial é relativamente recente, tendo seu início nos Estados Unidos no final da década 1990. A partir dos resultados obtidos com a aplicação do procedimento de fraturamento hidráulico, como método de estimulação, o estudo e avaliação de outros prospectos shale gas em outras regiões do mundo foi estimulado. No entanto, a análise, estudo e caracterização deste tipo de reservatórios são dificultados uma vez que devem ser levados em conta múltiplos fatores relacionados a geologia, mineralogia, petrofísica, geoquímica, entre outros. O fraturamento hidráulico é um processo hidromecânico acoplado com alto grau de dificuldade especialmente em reservatórios shale gas, onde existe a presença de fraturas. Um estudo numérico é conduzido neste trabalho com o objetivo de investigar o comportamento hidromecânico de uma fratura natural durante a injeção de fluido. O software UDEC (Universal Distinct Element Code) baseado no método dos elementos discretos foi empregado no desenvolvimento da modelagem numérica. O UDEC tem a capacidade de modelar o comportamento hidromecânico de uma fratura incluindo os fenômenos de abertura/fechamento das fraturas, cisalhamento e dilatação sob condições de contato ou separação. Nesta investigação numérica foram criados modelos de reservatório com fraturas naturais aleatórias por meio de uma distribuição voronoi visando representar a complexidade da rede de fraturas naturais presentes nas formações shale gas. Os resultados do estudo numérico mostram que o comportamento hidromecânico do sistema é fortemente dependente da variação de parâmetros como tensões in situ, viscosidade do fluido de fraturamento e taxa de injeção de fluido. Portanto, estes resultados proporcionam um melhor entendimento dos mecanismos de fraturamento e resposta da pressão de um tratamento de fraturamento hidráulico em um reservatório de gás não-convencional naturalmente fraturado. / [en] The development and deployment of shale gas formations around the world are relatively recent, starting in the United Stated in the late 1990. From the results obtained with the application of hydraulic fracturing as a method of stimulation, the study and evaluation of other prospects of shale gas in others places in the world was encouraged. However, the analysis, study and characterization of this type of reservoirs are difficult, because it must be taken into account several factors such as geology, mineralogy, petrophysics, geochemistry among others. Hydraulic fracturing is a complicated hydro-mechanical coupled process, with high difficulty degree especially in shale gas reservoir, where natural fractures exist. A numerical study is conducted to investigate the hydromechanical behavior of a natural fracture during fluid injection. UDEC (Universal Distinct Element Code) software based on discrete elements method was employed to numerical modeling development. UDEC has the ability to model the hydro-mechanical behavior of a fracture including phenomena like fracture enlargement, closure, slippage, and dilation under contact or separation condition. In this numerical investigation, numerical reservoir models, with random natural fractures through a distribution voronoi were created aiming to represent the network complexity of natural fractures present in shale gas formations. The numerical study results show that the hydromechanical system behavior is strongly dependent on the parameters variation such as in situ stress, fluid fracturing viscosity and fluid injection rate. Therefore, these results provide a better understanding of fracturing mechanisms and pressure response of a hydraulic fracturing treatment in a non- conventional naturally fractured reservoir.

L'objectif de cette thèse est de mettre en place un cadre numérique qui permet de simuler la réponse de structures en béton armé de taille industrielle soumises à des impacts sévères. Notre modèle repose sur trois ingrédients principaux : 1) la modélisation éléments discrets (ED) du béton, permettant de profiter de la nature « discrète » de cette formulation et de décrire facilement l'apparition et la propagation des discontinuités fortes de la matière ; 2) la modélisation éléments finis (EF) poutre des armatures, donnant la possibilité de représenter toute la complexité du ferraillage que l'on rencontre dans les structures industrielles en béton armé ; 3) un modèle de liaison acier-béton original, que nous avons proposé et mis en œuvre dans le code EUROPLEXUS, et qui constitue le principal apport de cette thèse. Nous avons effectué une étude théorique et numérique sur la stabilité et la précision de ce modèle de liaison, et avons calibré ses paramètres en simulant un essai d'arrachement. En modélisant l'essai de traction d'un tirant (barre en béton contenant une armature), nous avons montré la capacité de notre modèle à reproduire le transfert des efforts entre le béton et l'armature au niveau de l'interface acier-béton. Testé et validé sur des cas-tests élémentaires, le modèle mixte ED-EF a ensuite été appliqué à la simulation de vraies structures en béton armé. Nous avons simulé de manière détaillée quelques essais connus (l'impact mou sur poutre CEA, l'essai Meppen II-4), ce qui a permis de valider l'ensemble de notre approche numérique et de recueillir des éléments pour définir les pistes d'amélioration des modèles que nous avons mis en œuvre. / The aim of this work is to set up a numerical framework to simulate the behaviour of industrial size reinforced concrete structures subjected to severe impacts. Our model is based on three main features : 1) modeling of the concrete with a discrete method to handle easily strong material discontinuities such as initiation and propagation of macro-cracks ; 2) modeling of the reinforcement bars with finite element method to be able to represent complex reinforcement cages of industrial structures ; 3) an original steel-concrete bond model that we proposed and implemented in the dynamic explicit code EUROPLEXUS and that constitutes the main task of this work. We studied theoretically and numerically this model to guarantee its stability and precision during the time integration. We calibrated its parameters by simulating pull-out tests. To verify our model, we simulated a tie-test (a long concrete column which contains a reinforcement bar) and showed that our model ensures the correct transfert of forces between steel and concrete. After testing our model on simple benchmarks, we simulated real reinforced concrete structures subjected to impacts (soft impact on a beam, Meppen test n°II-4), allowing us to validate our numerical approach and to define some perspectives to improve the models we have developped.

La fracturation hydraulique est au cœur d'un certain nombre de phénomènes naturels et induits et est cruciale pour un développement durable de la production de ressources énergétiques. Compte tenu de son rôle crucial, ce phénomène a été pris en compte au cours des trois dernières décennies par le monde académique. Néanmoins, un certain nombre d'aspects très importants de ce processus ont été systématiquement négligés par la communauté. Deux des plus remarquables sont l'incapacité de la grande majorité des modèles existants à aborder la propagation des fractures hydrauliques dans les massifs rocheux fracturés où l'injection de fluide peut à la fois conduire à la fracturation de la roche intacte et à la réactivation de fractures préexistantes. Un autre aspect essentiel de ce processus est qu'il est intrinsèquement tridimensionnel, ce qui est souvent négligé par les modèles actuellement disponibles. Pour aborder ce problème essentiel, un modèle hydro-mécanique couplé basé sur la méthode des éléments discrets a été développé. La masse rocheuse est ici représentée par un ensemble d'éléments discrets interagissant à travers des lois de contact cohésifs qui peuvent se casser pour former des fissures à l'intérieur du milieu simulé. Ces fissures peuvent se coalescer pour former des fractures. Une méthode de volume fini est utilisée pour simuler l'écoulement de fluide entre les éléments discrets. L'écoulement est calculé en fonction de la déformation de l'espace poreux dans le milieu intact et de l'ouverture des fissures dans les fractures. De plus, les fractures naturelles sont modélisées explicitement de sorte qu'elles peuvent présentées des comportements mécanique et hydraulique différents de ceux de la matrice rocheuse intacte. La simulation des processus de fracturation hydraulique dans un milieu initialement intact en considérant plusieurs points d'injection plus ou moins espacés a permis de mettre en évidence l'évolution spatio-temporelle des fractures hydrauliques et de quantifier l'impact des différentes stratégies d'injection sur des indices représentatifs du volume fracturé, de l'intensité et de la densité des fractures ou encore sur la pression de fluide au niveau du puits. De plus, l'injection dans une fente de perforation non alignée sur le plan de contrainte minimum a génère des fractures hydrauliques non planaires percolantes si la connectivité est faible, ce qui peut être gênant pour la mise en place du proppant. En outre, des interactions fortes prennent place entre des fractures hydrauliques étroitement espacées ont été mises en évidence grâce au le suivi de la orientation de contrainte principale locale et ont révélé l'importance des effets d'ombre de contrainte. Des solutions sont proposées pour optimiser les traitements multiples à partir d'un puits de forage non parfaitement aligné. Enfin, l'interaction entre une seule fracture hydraulique et une seule fracture naturelle de propriétés et d'orientations variables a été étudiée à l'aide du modèle proposé. L'évolution de la fracture hydraulique et la réponse globale de l'échantillon ont été enregistrées d'une manière comparable aux données expérimentales existantes pour établir un pont entre les résultats expérimentaux et numériques. Les fractures naturelles persistantes semblent être des barrières pour la fracture hydraulique si leur conductance est élevée par apport a celle de la matrice ou si leur raideur est faible par rapport a la rigidité du milieu environnant. D'autre part, une faible rigidité dans les discontinuités non persistantes pourrait provoquer une bifurcation de la fracture hydraulique principale. De plus, des angles d'approche élevés et des contraintes différentielles fortes semblent favoriser le croisement de la fracture naturelle alors que des angles faibles engendrent plutôt un glissement ou une dilatation par cisaillement de la partie du plan qui n'est pas affectée par la perturbation de la contrainte. / Hydraulic fracturing is at the core of a number of naturally occurring and induced phenomena and crucial for a sustainable development of energy resource production. Given its crucial role this process has been given increasing attention in the last three decades from the academic world. Nonetheless a number of very significant aspects of this process have been systematically overlooked by the community. Two of the most notable ones are the inability of the vast majority of existing models to tackle at once the propagation of hydraulic fractures in realistic, fractured rocks-masses where hydraulic fracturing is a competing dipole mechanism between fracturing of the intact rock and re-activation of exiting fracture networks. Another essential aspect of this process is that it is intrinsically three-dimensional which is neglected by most models. To tackle this vital problem taking into account these pivotal aspects, a fully coupled hydro-mechanical model based on the discrete element method has been developed. The rock mass is here represented by a set of discrete elements interacting through elastic-brittle bonds that can break to form cracks inside the simulated medium. Theses cracks can coalesce to form fractures. A finite volume scheme is used to simulate the fluid flow in between these discrete elements. The flow is computed as a function of the pore space deformation in the intact medium and of the cracks' aperture in the fractures. Furthermore, the natural fractures are modelled explicitly and present mechanical and hydraulic properties different from the rock matrix. Employing this model in an intact numerical specimen, single fluid injection and multiple closely spaced sequential injections, enabled the description the full spatio-temporal evolution of HF propagation and its impact on quantitative indexes used in description of hydraulic fracturing treatments, such as fractured volume, fracture intensity and down-the-hole pressure for different control parameters and in-situ stress-fields. Moreover, injections from perforation slots which are not well aligned to the minimum stress plane showed possible creation of percolating non-planar hydraulic fractures of low connectivity, which can be troublesome for proppant placement. Also, strong interactions between closely spaced HF were highlighted by tracking the local principal stress rotation around the injection zones, emphasizing the importance of stress shadow effects. Optimization solutions are proposed for multiple treatments from a non-perfectly aligned wellbore. Finally, interaction between a single hydraulic fracture and a single natural fracture of varying properties and orientations was studied using the proposed model. The evolution of the hydraulic fracture and the global response of the specimen were recorded in a way comparable to existing experimental data to bridge the experimental and numerical findings. Persistent natural fractures appeared to be barriers for the hydraulic fracture if their conductance is high compared to the matrix conductivity or if their stiffness is significantly low compared to the rock matrix rigidity. Low stiffness in non-persistent defects might also cause a bifurcation of the main hydraulic fracture due to the local stress field perturbation around the defect and ahead of the hydraulic fracture tip. Furthermore, high approach angles and differential stresses seemed to favour crossing of the natural fracture while low angles enable shear slippage or dilation on the part of the plane which is not affected by the local stress perturbation.