Anisotropy of hydraulic properties of the subsurface geologic formations is an essential feature that has been established as a consequence of the different geologic processes that undergo during the longer geologic time scale. With respect to subsurface reservoirs, in many cases, anisotropy plays significant role in dictating the direction of flow that becomes no longer dependent only on driving forces like the pressure gradient and gravity but also on the principal directions of anisotropy. Therefore, there has been a great deal of motivation to consider anisotropy into the subsurface flow and transport models. In this dissertation, we present subsurface flow modeling in single and dual continuum anisotropic porous media, which include the single-phase groundwater flow coupled with the solute transport in anisotropic porous media, the two-phase flow with gravity effect in anisotropic porous media, and the natural gas flow in anisotropic shale reservoirs. We have employed the multipoint flux approximation (MPFA) method to handle anisotropy in the flow model. The MPFA method is designed to provide correct discretization of the flow equations for general orientation of the principal directions of the permeability tensor. The implementation of MPFA method is combined with the experimenting pressure field approach, a newly developed technique that enables the solution of the global problem breaks down into the solution of multitude of local problems. The numerical results of the study demonstrate the significant effects of anisotropy of the subsurface formations. For the single-phase groundwater flow coupled with the solute transport modeling in anisotropic porous media, the results shows the strong impact of anisotropy on the pressure field and the migration of the solute concentration. For the two-phase flow modeling with gravity effect in anisotropic porous media, it is observed that the buoyancy-driven flow, which emerges due to the density differences between the phases, migrates upwards and the anisotropy aligns the flow directions closer to the principal direction of anisotropy. Lastly, for the gas flow modeling in anisotropic shale reservoirs, we observe that anisotropy affects the pressure fields and the velocity fields of the matrix and fracture systems as well as the production rate and cumulative production. It is observed from the results that all of the anisotropic cases produce higher amount of gas compared to isotropic case during the same production time. Furthermore, we have also examined the performance of MPFA with respect to mixed finite element (MFE) method over the lowest-order Raviart-Thomas (RT0) space and the first-order Brezzi-Douglas-Marini (BDM1) space. From the comparison of the numerical results we observe that MPFA method show very good agreement with the BDM1 than RT0. In terms of numerical implementation, however, MPFA method is easier than BDM1 and it also offers explicit discrete fluxes that are advantageous. Combining MPFA with the experimenting pressure field approach will certainly adds another advantage of implementing MPFA method as compared with RT0 and BDM1. Moreover, the computational cost (CPU cost) of the three different methods are also discussed.

A large percentage of oil and gas reservoirs in the most productive regions such as the Middle East, South America, and Southeast Asia are naturally fractured reservoirs (NFR). The major difference between conventional reservoirs and naturally fractured reservoirs is the discontinuity in media in fractured reservoir due to tectonic activities. These discontinuities cause remarkable difficulties in describing the petrophysical structures and the flow of fluids in the fractured reservoirs. Predicting fluid flow behavior in naturally fractured reservoirs is a challenging area in petroleum engineering. Two classes of models used to describe flow and transport phenomena in fracture reservoirs are discrete and continuum (i.e. dual porosity) models. The discrete model is appealing from a modeling point of view, but the huge computational demand and burden of porting the fractures into the computational grid are its shortcomings. The affect of natural fractures on the permeability anisotropy can be determined by considering distribution and orientation of fractures. Representative fracture permeability, which is a crucial step in the reservoir simulation study, must be calculated based on fracture characteristics. The diagonal representation of permeability, which is customarily used in a dual porosity model, is valid only for the cases where fractures are parallel to one of the principal axes. This assumption cannot adequately describe flow characteristics where there is variation in fracture spacing, length, and orientation. To overcome this shortcoming, the principle of the full permeability tensor in the discrete fracture network can be incorporated into the dual porosity model. Hence, the dual porosity model can retain the real fracture system characteristics. This study was designed to develop a novel approach to integrate dual porosity model and full permeability tensor representation in fractures. A fully implicit, parallel, compositional chemical dual porosity simulator for modeling naturally fractured reservoirs has been developed. The model is capable of simulating large-scale chemical flooding processes. Accurate representation of the fluid exchange between the matrix and fracture and precise representation of the fracture system as an equivalent porous media are the key parameters in utilizing of dual porosity models. The matrix blocks are discretized into both rectangular rings and vertical layers to offer a better resolution of transient flow. The developed model was successfully verified against a chemical flooding simulator called UTCHEM. Results show excellent agreements for a variety of flooding processes. The developed dual porosity model has further been improved by implementing a full permeability tensor representation of fractures. The full permeability feature in the fracture system of a dual porosity model adequately captures the system directionality and heterogeneity. At the same time, the powerful dual porosity concept is inherited. The implementation has been verified by studying water and chemical flooding in cylindrical and spherical reservoirs. It has also been verified against ECLIPSE and FracMan commercial simulators. This study leads to a conclusion that the full permeability tensor representation is essential to accurately simulate fluid flow in heterogeneous and anisotropic fracture systems. / text

Les circulateurs et les isolateurs à ferrite sont couramment utilisés dans l’électronique hyperfréquence en raison de leur forte résistivité électrique et de leur aimantation spontanée élevée. La conception et l’optimisation des dispositifs micro-ondes à ferrites nécessitent d’une part la connaissance de leurs propriétés dynamiques, permittivité complexe et tenseur de perméabilité, et d’autre part le contrôle de la propagation de l’onde électromagnétique (EM) qui conditionne leurs performances. Les logiciels commerciaux de simulation utilisent différents modèles théoriques pour décrire le tenseur de perméabilité en fonction de l’état d’aimantation. Cependant la plupart de ces simulateurs EM restent limités à des états particuliers d’aimantation en raison des hypothèses simplificatrices des modèles de perméabilité utilisés. Dans ce travail de thèse, nous présentons un outil prédictif pour l’étude des propriétés EM des ferrites quel que soit leur état d’aimantation et qui tient compte de l’inhomogénéité des champs internes de polarisation. Cette modélisation combine des techniques expérimentales de détermination des paramètres physiques des ferrites et un modèle théorique qui utilise ces paramètres pour décrire le comportement dynamique des ferrites quel que soit l’état d’aimantation. Dans la première partie de la thèse nous présentons une méthode large bande en ligne coaxiale pour la mesure du coefficient d’amortissement. Les paramètres S théoriques sont calculés à partir d’une analyse EM (problème direct) de la cellule de mesure. Pour le problème inverse, une optimisation numérique a été développée pour calculer le coefficient d’amortissement (α) par comparaison des paramètres S calculés avec ceux mesurés. Dans la seconde partie de la thèse, nous présentons un outil théorique de modélisation EM qui combine une analyse magnétostatique, un modèle du tenseur de perméabilité généralisé (GPT) et le simulateur Ansys HFSSTM. La majorité des paramètres d’entrée comme l’aimantation à saturation ou le champ d’anisotropie peuvent être mesurés à l’aide de techniques standards de caractérisation statique. Seul le paramètre dynamique, le coefficient d’amortissement, est déterminé à l’aide de la technique en ligne coaxiale proposée dans la première partie de la thèse. L’outil théorique développé est ensuite validé par la modélisation et la réalisation d’un circulateur micro-ruban à jonction Y. Grâce à la prise en compte de l’inhomogénéité des champs internes de polarisation, l’outil théorique proposé permet de mieux prédire le comportement dynamique des dispositifs à ferrites et cela pour tout état d’aimantation. / Ferrites are widely used in microwave electronics, particularly for circulators and insulators, because of their high electrical resistivity and high spontaneous magnetization. Design and optimization of microwave devices using ferrites requires realistic knowledge of its dynamic response, namely complex permittivity and permeability tensor and, on the other hand, control of wave propagation that condition their performance. Commercial simulation software use different theoretical models to describe the permeability tensor according to the state of magnetization. However, most of the electromagnetic (EM) simulators remain limited to certain states of magnetization, due to the simplified assumptions on which their permeability models are based upon.In this thesis work, we presented a predictive electromagnetic tool to study the EM properties of ferrites, whatever their magnetization state is, and takes into account the inhomogeneity of the internal polarization fields. This theoretical modeling approach combines experimental techniques to find the physical parameters of the ferrites, and a theoretical model which will use these parameters to describe the dynamic behavior of ferrites at any magnetization state.In the first part of the thesis, we presented a broadband coaxial line method for damping factor measurement. Theoretical S parameters are calculated using the EM analysis (direct problem) of the measurement cell. In the inverse problem, a numerical optimization procedure is developed to compute the damping factor (α) by matching theoretical S parameters with measured S parameters.During the second part of the thesis, we developed a theoretical EM modeling tool which combines a magneto-static solver, generalized permeability tensor model and commercial simulation software Ansys HFSSTM. Most of the input parameters like saturation magnetization, anisotropy field, etc. can be measured using standard characterization methods, except the damping factor used to represent the dynamic losses. Static input parameters of this theoretical tool are determined using standard material characterization methods.Dynamic input parameter, damping factor is calculated using the coaxial line technique proposed in the first part of this thesis. Theoretical EM tool is validated by modeling, and realizing a microstrip Y-junction circulator. By taking into account the inhomogeneity of the internal polarizing fields, proposed theoretical tool can predict the dynamic behavior of ferrite devices more accurately, at all magnetization states.

Les phénomènes de transport en milieux poreux sont étudiés depuis près de deux siècles, cependant les travaux concernant les milieux fortement poreux sont encore relativement peu nombreux. Les modèles couramment utilisés pour les poreux classiques (lits de grains par exemple) sont peu applicables pour les milieux fortement poreux (les mousses par exemple), un certain nombre d’études ont été entreprises pour combler ce manque. Néanmoins, les résultats expérimentaux et numériques caractérisant les pertes de charge dans les mousses sont fortement dispersés. Du fait des progrès de l’imagerie 3D, une tendance émergente est la détermination des paramètres des lois d’écoulement à partir de simulations directes sur des géométries reconstruites. Nous présentons ici l’utilisation d’une nouvelle approche cinétique pour résoudre localement les équations de Navier-Stokes et déterminer les propriétés d’écoulement (perméabilité, dispersion, ...). / A novel kinetic scheme satisfying an entropy condition is developed, tested and implemented for the simulation of practical problems. The construction of this new entropic scheme is presented. A classical hyperbolic system is approximated by a discrete velocity vector kinetic scheme (with the simplified BGK collisional operator), but applied to an inviscid compressible gas dynamics system with a small Mach number parameter, according to the approach of Carfora and Natalini (2008). The numerical viscosity is controlled, and tends to the physical viscosity of the Navier-Stokes system. The proposed numerical scheme is analyzed and formulated as an explicit finite volume flux vector splitting (FVS) scheme that is very easy to implement. It is close in spirit to Lattice Boltzmann schemes, but it has the advantage to satisfy a discrete entropy inequality under a CFL condition and a subcharacteristic stability condition involving a cell Reynolds number. The new scheme is proved to be second-order accurate in space. We show the efficiency of the method in terms of accuracy and robustness on a variety of classical benchmark tests. Some physical problems have been studied in order to show the usefulness of both schemes. The LB code was successfully used to determine the longitudinal dispersion of metallic foams, with the use of a novel indicator. The entropic code was used to determine the permeability tensor of various porous media, from the Fontainebleau sandstone (low porosity) to a redwood tree sample (high porosity). These results are pretty accurate. Finally, the entropic framework is applied to the advection-diffusion equation as a passive scalar.