Spontaneous imbibition plays a very important role in the displacement mechanism of non-wetting fluid in naturally fractured reservoirs. We developed a new 2D two-phase finite element numerical model, as available commercial simulators cannot be used to model small-scale experiments with different boundary conditions as well as complex boundary conditions such as fractures and vugs. Starting with the basic equation of fluid flow, we derived the non-linear diffusion saturation equation. This equation cannot be put in weighted-integral weak variational form and hence Rayleigh-Ritz finite element method (FEM) cannot be applied. Traditionally, the way around it is to use higher order interpolation functions and use Galerkin FEM or reduce the differentiability requirement and use Mixed FEM formulation. Other FEM methods can also be used, but iterative nature of those methods makes them unsuitable for solving large-scale field problems. But if we truncate the non-linear terms and decouple the dependent variables, from the spatial as well as the temporal domains of the primary variable to solve them analytically, the non-linear FEM problem reduces to a simple weighted integral form, which can be put into its corresponding weak form. The advantage of using Rayleigh-Ritz method is that it has immediate effect on the computation time required to solve a particular problem apart from incorporating complex boundary conditions. We compared our numerical models with the analytical solution of this diffusion equation. We validated the FDM numerical model using X-Ray Tomography (CT) experimental data from the single-phase spontaneous imbibition experiment, where two simultaneously varying parameters of weight gain and CT water saturation were used and then went ahead and compared the results of FEM model to that of FDM model. A two-phase field size example was taken and results from a commercial simulator were compared to the FEM model to bring out the limitations of this approach.

Spontaneous imbibition plays a very important role in the displacement mechanism of non-wetting fluid in naturally fractured reservoirs. We developed a new 2D two-phase finite element numerical model, as available commercial simulators cannot be used to model small-scale experiments with different boundary conditions as well as complex boundary conditions such as fractures and vugs. Starting with the basic equation of fluid flow, we derived the non-linear diffusion saturation equation. This equation cannot be put in weighted-integral weak variational form and hence Rayleigh-Ritz finite element method (FEM) cannot be applied. Traditionally, the way around it is to use higher order interpolation functions and use Galerkin FEM or reduce the differentiability requirement and use Mixed FEM formulation. Other FEM methods can also be used, but iterative nature of those methods makes them unsuitable for solving large-scale field problems. But if we truncate the non-linear terms and decouple the dependent variables, from the spatial as well as the temporal domains of the primary variable to solve them analytically, the non-linear FEM problem reduces to a simple weighted integral form, which can be put into its corresponding weak form. The advantage of using Rayleigh-Ritz method is that it has immediate effect on the computation time required to solve a particular problem apart from incorporating complex boundary conditions. We compared our numerical models with the analytical solution of this diffusion equation. We validated the FDM numerical model using X-Ray Tomography (CT) experimental data from the single-phase spontaneous imbibition experiment, where two simultaneously varying parameters of weight gain and CT water saturation were used and then went ahead and compared the results of FEM model to that of FDM model. A two-phase field size example was taken and results from a commercial simulator were compared to the FEM model to bring out the limitations of this approach.

Igneous and metamorphic rock units have long been considered marginal aquifers yet they are a significant source for potable drinking water in many areas worldwide. Additionally, use of these systems is on the rise due to many factors including, contamination and overuse of surficial systems, as well as expanding population and drought. The Nashoba Terrane is a fault-bounded block of high-grade, steeply dipping metavolcanic and metasedimentary rock located in eastern Massachusetts, U.S.A. The Nashoba is northeast trending, extending from Oxford, MA to the Gulf of Maine south of Newburyport, MA. Seventeen previously drilled wells throughout the Nashoba were selected for use in this study. The goal of this study was to characterize the hydrogeologic system of the Nashoba Terrane. Wells studied were in three bedrock types: granite, schist and amphibolite. Three fracture types were identified: FPF, subhorizontal unloading joints and tectonic joints. Several major fracture orientation sets were also identified including northeast trending FPF, east-west trending and north-south trending tectonic joints as well as northwest trending tectonic joints. Dominant sets varied in the three rock types and the frequency of fractures was found to decrease with depth. Only four percent of all fractures measured in this study were flowing. Approximately 32% of the flowing fractures were northeast trending, 17% of subhorizontal fractures were flowing and the remaining 51% were of variable orientation and dip. In general, the orientation of fractures was not found to determine whether a fracture flows, nor was rock type a significant determinant of flow. There was no flow identified below 170 meters and the majority of flow in the Nashoba Terrane is constrained to the upper 100 meters. This is most likely due to decreased fracture frequency and permeability with depth. This study is significant to the search for a sustainable groundwater source in bedrock because results show that the few fractures are actually contributing to flow and that flow is primarily occurring near the surface.

Les travaux de cette thèse portent sur la modélisation et la discrétisation des écoulements diphasiques dans les milieux poreux fracturés. On se place dans le cadre des modèles dits dimensionnels hybrides couplant l'écoulement dans la matrice 3D à l'écoulement dans un réseau de fractures modélisées comme des surfaces 2D. La discrétisation s'appuie sur le cadre abstrait des schémas gradients. Dans cette étude nous présentons deux classes de schémas de types Gradient Schemes sur ces modèles en monophasique et en diphasique. Les objectifs sont motivés par l'application cible de la thèse qui concerne les procédés de récupération assistée de gaz par fracturation hydraulique dans les réservoirs de très faibles perméabilités / This thesis presents the work on modelling and discretisation of two-phase flows in the fractured porous media. These models couple the flow in the fractures represented as the surfaces of codimension one with the flow in the surrounding matrix. The discretisation is made in the framework of Gradient schemes which accounts for a large family of conforming and nonconforming discretizations. The test cases are motivated by the target application of the thesis concerning the gas recovery under the hydraulic fracturing process in low-permeability reservoirs