Orientador: Caio Gorla Nogueira / Resumo: Este trabalho apresenta um conjunto de técnicas estatísticas básicas aplicadas à modelagem de problemas de fluxo em meios porosos fraturados e de estabilidade de taludes, com o objetivo de identificar as variáveis explicativas mais influentes sobre a variabilidade das variáveis resposta. Diferentes planejamentos de experimentos foram utilizados para possibilitar a construção de metamodelos polinomiais representativos dos fenômenos estudados. Uma modificação do planejamento do tipo Box-Behnken é apresentada e foi proposta pelo autor para analisar problemas que envolvem elevado número de variáveis explicativas (e.g. 30). Os metamodelos, obtidos via método dos mínimos quadrados, são também chamados de superfícies de resposta ou modelos de regressão e são indispensáveis à verificação da sensibilidade das variáveis explicativas. O conjunto de técnicas mostrou- se muito eficaz na identificação das variáveis explicativas que provocaram efeitos mais significativos sobre a variável resposta. Evidenciou-se também, por meio dos exemplos de estabilidade de taludes tratados, a possibilidade de se quantificar incertezas com o uso de metamodelos suficientemente adequados, uma opção que pode ser bastante útil no processo de quantificação de incertezas de problemas que não possuem soluções analíticas simples. / Abstract: This paper presents a set of basic statistical techniques applied to the modeling of flow problems in fractured porous and slope stability media, aiming to identify the most influential explanatory variables on the response variables variability. Different designs of experiments were used to enable the construction of polynomial metamodels representative of the studied phenomena. A Box-Behnken type design modification is presented and was proposed by the author to analyze problems involving high number of explanatory variables (e.g. 30). The metamodels, obtained by the least squares method, are also called response surfaces or regression models and are indispensable for verifying the sensitivity of the explanatory variables. The set of techniques was very effective in identifying the explanatory variables that had the most significant effects on the response variable. It was also evidenced, through the examples of stability of treated slopes, the possibility of quantifying uncertainties using sufficiently adequate metamodels, an option that can be very useful in the process of quantifying uncertainty of problems that do not have simple analytical solutions. / Mestre

Nous construisons une méthode numérique fiable pour simuler un écoulement dans un milieu poreux modélisé par une équation elliptique. La simulation est rendue difficile par les hétérogénéités du milieu, la taille et la géométrie complexe du domaine de calcul. Un maillage d'hexaèdres réguliers ne permet pas de représenter fidèlement les couches géologiques du domaine. Par conséquent, nous sommes amenés à travailler avec un maillage de cubes déformés. Il existe différentes méthodes de volumes finis ou d'éléments finis qui résolvent ce problème avec plus ou moins de succès. Pour la méthode que nous proposons, nous nous imposons d'avoir seulement un degré de liberté par maille pour la pression et un degré de liberté par face pour la vitesse de Darcy, pour rester au plus près des habitudes des codes industriels. Comme les méthodes d'éléments finis mixtes standards ne convergent pas, notre méthode est basée sur un élément fini mixte composite. En deux dimensions, une maille polygonale est découpée en triangles en ajoutant un point au barycentre des sommets, et une expression explicite des fonctions de base a pu être obtenue. En dimension 3, la méthode s'étend naturellement au cas d'une maille pyramidale. Dans le cas d'un hexaèdre ou d'un cube déformé quelconque, la maille est divisée en 24 tétraèdres en ajoutant un point au barycentre des sommets et en divisant les faces en 4 triangles. Les fonctions de base de l'élément sont alors construites en résolvant un problème discret. Les méthodes proposées ont été analysées théoriquement et complétées par des estimateurs a posteriori. Elles ont été expérimentées sur des exemples académiques et réalistes en utilisant le calcul parallèle. / We develop a reliable numerical method to approximate a flow in a porous media, modeled by an elliptic equation. The simulation is made difficult because of the strong heterogeneities of the medium, the size together with complex geometry of the domain. A regular hexahedral mesh does not allow to describe accurately the geological layers of the domain. Consequently, this leads us to work with a mesh made of deformed cubes. There exists several methods of type finite volumes or finite elements which solve this issue. For our method, we wish to have only one degree of freedom per element for the pressure and one degree of freedom per face for the Darcy velocity, to stay as close to the habits of industrial software. Since standard mixed finite element methods does not converge, our method is based on composite mixed finite element. In two dimensions, a polygonal mesh is split into triangles by adding a node to the vertices's barycenter, and explicit formulation of the basis functions was obtained. In dimension 3, the method extend naturally to the case of pyramidal mesh. In the case of a hexahedron or a deformed cube, the element is divided into 24 tetrahedra by adding a node to the vertices's barycenter and splitting the faces into 4 triangles. The basis functions are then built by solving a discrete problem. The proposed methods have been theoretically analyzed and completed by a posteriori estimators. They have been tested on academical and realistic examples by using parallel computation.

Experimental data convincingly show that viscosity of a fluid may change significantly with pressure. This observation leads to various generalizations of well-known models, like Darcy's law, Stokes' law or the Navier-Stokes equations, among others. This thesis investigates three such models in a series of three published papers. Their unifying topic is development of existence theory and finding a weak solution to systems of partial differential equations stemming from the considered models.

The theory of homogenization allows to find for a given system of partial differential equations governing a model with a very complicated internal struc- ture a system governing a model without this structure, whose solution is in a certain sense an approximation of the solution of the original problem. In this thesis, methods of the theory of homogenization are applied to three sys- tems of partial differential equations. The first one governs a flow of a class of non-Newtonian fluid through a porous medium. The second system is utilized for modeling of a flow of a fluid through an electric field wherein the viscosity depends significantly on the intensity of the electric field. For the third system is considered an elliptic operator having growth and coercivity indicated by a general anisotropic inhomogeneous N-function. 1

No description available.

[pt] A injeção de emulsão em meio poroso como método de recuperação avançada pode se tornar realidade na operação de campos de petróleo devido à maior rigidez no descarte de água produzida e aos potenciais ganhos na produção de óleo. Para entender o comportamento macroscópico desta técnica de EOR, é necessário compreender o fenômeno microscópico. Com este objetivo, fez-se a modelagem numérica do escoamento de uma gota imersa em fase contínua escoando em capilar restrito utilizando o método de Lattice-Boltzmann. Este método foi escolhido devido à sua facilidade de ser aplicado em geometrias complexas de rocha e fluido e ao bom compromisso na representação dos fenômenos de microescala. Para entender a influência de cada parâmetro, foram realizadas diversas simulações em domínio tridimensional, alterando a velocidade do fluxo, a razão de viscosidades dos fluidos, a relação entre os diâmetros da gota e do tubo e a magnitude da tensão interfacial. Observou-se que a passagem da gota pela restrição causa uma redução na mobilidade do escoamento, representada por um aumento na perda de carga, pela conjunção dos efeitos viscosos e capilares. Obteve-se correlação negativa do fator redutor de mobilidade 𝑓 com a razão de viscosidade e com o tamanho da gota, assim como fora determinado numericamente por Roca-Reyes (2011) com o método level-set. Foi notada uma pequena sensibilidade de 𝑓 ao número de capilaridade, assim como estabelecido experimentalmente por Robles-Castillo (2011). Verificou-se a importância de se determinar o conjunto adequado de parâmetros do sistema para ter sucesso na implantação de injeção de emulsões. / [en] Emulsion injection in porous medium as an Enhanced Oil Recovery method can turn out to be reality in the operation of onshore and offshore fields, due to increasing rigidity in the disposal of produced water and also due to the potential additional oil production. In order to understand macroscopic behavior of this EOR method, it is necessary to understand the microscopic phenomenon. With this objective, it was performed the numerical modeling of the flow of a droplet immerse in continuous phase through a constricted capillary using the Lattice-Boltzmann method. This method was chosen due to its simplicity on being applied to complex rock geometries and multiphasic flow and due to its good commitment in representing microscopic phenomena. Focusing on understanding the influence of each parameter on flow behavior, several simulation studies were performed altering flow velocity, viscosity ratio, ratio between droplet s and pipe s diameter and interfacial tension. A reduction in mobility is observed as the droplet passes through the restriction caused by the conjunction of viscous and capillary effects. A negative correlation of mobility reduction factor 𝑓 in relation to the viscosity ratio and to droplet size was noticed, as it had been observed before by Roca-Reyes (2011) in a numerical implementation of level-set method. Weak correlation between 𝑓 and capillary number was determined, as in previous experimental essay performed by Robles-Castillo (2011). In this study, it was verified the importance of determining the appropriate set of system parameters, in order to achieve success in the implementation of emulsion injection.

Extremely small pore size, low porosity, and ultra-low permeability are among the characteristics of shale rocks. In tight shale reservoirs, the nano-confinement effects that include large gas-oil capillary pressure and critical property shifts could alter the phase behaviors, thereby affecting the oil or gas production. In this research, two in-house simulation models, i.e., a compositionally extended black-oil model and a fully composition model are developed to examine the nano-pore confinement effects on convective and diffusive multicomponent multiphase fluid transport. Meanwhile, the effect of nano-confinement and rock intrinsic properties (porosity and tortuosity factor) on predicting effective diffusion coefficient are investigated. First, a previously developed compositionally extended black-oil simulation approach is modified, and extended, to include the effect of large gas-oil capillary pressure for modeling first contact miscible (FCM), and immiscible gas injection. The simulation methodology is applied to gas flooding in both high and very low permeability reservoirs. For a high permeability conventional reservoir, simulations use a five-spot pattern with different reservoir pressures to mimic both FCM and immiscible displacements. For a tight oil-rich reservoir, primary depletion and huff-n-puff gas injection are simulated including the effect of large gas-oil capillary pressure in flow and in flash calculation on recovery estimations. A dynamic gas-oil relative permeability correlation that accounts for the compositional changes owing to the produced gas injection is introduced and applied to correct for changes in interfacial tension (IFT), and its effect on oil recovery is examined. The results show that the simple modified black-oil approach can model well both immiscible and miscible floods, as long as the minimum miscibility pressure (MMP) is matched. It provides a fast and robust alternative for large-scale reservoir simulation with the purpose of flaring/venting reduction through reinjecting the produced gas into the reservoir for EOR. Molecular diffusion plays an important role in oil and gas migration in tight shale formations. However, there are insufficient reference data in the literature to specify the diffusion coefficients within porous media. Another objective of this research is to estimate the diffusion coefficients of shale gas, shale condensate, and shale oil at reservoir conditions with CO2 injection for EOR/EGR. The large nano-confinement effects including large gas-oil capillary pressure and critical property shifts could alter the phase behaviors. This study estimates the diffusivities of shale fluids in nanometer-scale shale rock from two perspectives: 1) examining the shift of diffusivity caused by nanopore confinement effects from phase change (phase composition and fluid property) perspective, and 2) calculating the effective diffusion coefficient in porous media by incorporating rock intrinsic properties (porosity and tortuosity factor). The tortuosity is obtained by using tortuosity-porosity relations as well as the measured tortuosity of shale from 3D imaging techniques. The results indicated that nano-confinement effects could affect the diffusion coefficient through altering the phase properties, such as phase compositions and densities. Compared to bulk phase diffusivity, the effective diffusion coefficient in porous shale rock is reduced by 102 to 104 times as porosity decreases from 0.1 to 0.03. Finally, a fully compositional model is developed, which enables us to process multi-component multi-phase fluid flow in shale nano-porous media. The validation results for primary depletion, water injection, and gas injection show a good match with the results of a commercial software (CMG, GEM). The nano-confinement effects (capillary pressure effect and critical property shifts) are incorporated in the flash calculation and flow equations, and their effects on Bakken oil production and Marcellus shale gas production are examined. The results show that including oil-gas capillary pressure effect could increase the oil production but decrease the gas production. Inclusion of critical property shift could increase the oil production but decrease the gas production very slightly. The effect of molecular diffusion on Bakken oil and Marcellus shale gas production are also examined. The effect of diffusion coefficient calculated by using Sigmund correlation is negligible on the production from both Bakken oil and Marcellus shale gas huff-n-puff. Noticeable increase in oil and gas production happens only after the diffusion coefficient is multiplied by 10 or 100 times. / Doctor of Philosophy / Shale reservoir is one type of unconventional reservoir and it has extremely small pore size, low porosity, and ultra-low permeability. In tight shale reservoirs, the pore size is in nanometer scale and the oil-gas capillary pressure reaches hundreds of psi. In addition, the critical properties (such as critical pressure and critical temperature) of hydrocarbon components will be altered in those nano-sized pores. In this research, two in-house reservoir simulation models, i.e., a compositionally extended black-oil model and a fully composition model are developed to examine the nano-pore confinement effects on convective and diffusive multicomponent multiphase fluid transport. The large nano-confinement effects (large gas-oil capillary pressure and critical property shifts) on oil or gas production behaviors will be investigated. Meanwhile, the nano-confinement effects and rock intrinsic properties (porosity and tortuosity factor) on predicting effective diffusion coefficient are also studied.

Cette thèse est consacrée à l'étude de méthodes numériques pour la simulation des écoulements polyphasiques en milieu poreux, en vue de leur application à des problèmes d'ingénierie pétrolière ou environnementale. Nous présentons une formulation générique du modèle d'écoulements à nombre quelconque de composants présents dans un nombre quelconque de phases. Dans notre approche l'approximation des flux diffusifs (issus, par exemple, de la loi de Darcy) s'appuie sur de nouveaux schémas, appelés schémas gradient, qui ont plusieurs avantages sur les schémas industriels standard : ces derniers, qui sont des schémas volumes finis multi-points centrés aux mailles, ne sont généralement pas symétriques et convergent difficilement sur des cas à forts rapports d'anisotropie. Nous montrons en revanche que les schémas gradient conduisent naturellement à des approximations symétriques et convergentes. Parmi cette classe de schémas, nous étudions plus particulièrement le schéma "VAG" qui fait intervenir des inconnues au centre des mailles et aux sommets du maillage. Ce schéma conduit à la définition de flux entre le centre d'une maille et ses sommets, qui sont utilisés pour généraliser la méthode "VAG" au contexte polyphasique. Des tests numériques montrent alors que ce schéma est robuste, et conduit à un très bon compromis précision/coût, ce qui en fait un candidat idoine pour les applications industrielles. Nous présentons notamment un cas test, basé sur des observations de terrains, d'injection et de dissolution de CO2 dans la région proche d'un puits foré dans un aquifère salin. Nous montrons alors que le schéma numérique permet de simuler l'assèchement et la précipitation de minéral observée en pratique. Un chapitre de la thèse est enfin consacré à l'étude pratique et théorique d'une méthode numérique générique pour contrôler l'effet d'axe lors de l'utilisation de schémas industriels / This thesis is focused on numerical methods dedicated to the simulation of multiphase flow in porous media, involved in petroleum or environmental engineering. We present a generic formulation of the flow model which is able to take into account any number of components within any number of phases. In our approach the approximation of the diffusive fluxes (mainly resulting from Darcy's law) is based on new schemes, called gradient schemes, which show several advantages over the standard industrial numerical schemes : these schemes, which belong to the class of the cell-centred MultiPoint Flux Approximation finite volume schemes, are not symmetric and may lead to difficulties of convergence in the case of high anisotropy ratios. We indeed show that gradient schemes are naturally providing symmetric and convergent approximations. We particularly study the Vertex Approximate Gradient scheme (called the VAG scheme), whose discrete unknowns are the values at the cell centres and at the vertices of the mesh. This scheme implies the definition of fluxes between the centre of a given cell and its vertices, used for the extension of the scheme to multiphase flow. Numerical tests show the robustness and the accuracy of the method for a low computational cost, which enables the use of the VAG scheme in an industrial framework. A test case, based on experimental data of injection and dissolution of CO2 in the near-well region within a saline aquifer, shows the aptitude of the scheme for reproducing drying and salt precipitation, which are practically observed. Finally, a chapter of the thesis is devoted to the theoretical and practical study of a general numerical method for controlling Grid Orientation Effect in industrial simulators

Neste trabalho é apresentada uma nova modelagem matemática para a descrição do escoamento de um líquido incompressível através de um meio poroso rígido homogêneo e isotrópico, a partir do ponto de vista da Teoria Contínua de Misturas. O fenômeno é tratado como o movimento de uma mistura composta por três constituintes contínuos: o primeiro representando a matriz porosa, o segundo representando o líquido e o terceiro representando um gás de baixíssima densidade. O modelo proposto possibilita uma descrição matemática realista do fenômeno de transição insaturado/saturado a partir de uma combinação entre um sistema de equações diferenciais parciais e uma desigualdade. A desigualdade representa uma limitação geométrica oriunda da incompressibilidade do líquido e da rigidez do meio poroso. Alguns casos particulares são simulados e os resultados comparados com resultados clássicos, mostrando as consequências de não levar em conta as restrições inerentes ao problema. / This work is concerned with a new mathematical modelling for describing the flow of an incompressible fluid (a liquid) through a rigid, homogeneous and isotropic porous medium, from a Continuum Mixtures point of view. The phenomenon is regarded as the motion of a mixture composed by three overlaping continuous constituents: the first one representing the porous matrix, the second one representing the liquid and the third one representing a (very) low density gas. The proposed mathematical modelling allows a realistic mathematical description for the unsaturated/saturated transition process by means of a combination between a system of partial differential equations and an inequality. This inequality represents a geometrical constraint arising from the liquid incompressibility merged with the porous matrix rigidity. The simulation of some interesting particular cases is carried out presenting a comparison between the obtained results and the classical ones, showing the consequences of disregarding the constraints associated to the phenomenon.

Neste trabalho é apresentada uma nova modelagem matemática para a descrição do escoamento de um líquido incompressível através de um meio poroso rígido homogêneo e isotrópico, a partir do ponto de vista da Teoria Contínua de Misturas. O fenômeno é tratado como o movimento de uma mistura composta por três constituintes contínuos: o primeiro representando a matriz porosa, o segundo representando o líquido e o terceiro representando um gás de baixíssima densidade. O modelo proposto possibilita uma descrição matemática realista do fenômeno de transição insaturado/saturado a partir de uma combinação entre um sistema de equações diferenciais parciais e uma desigualdade. A desigualdade representa uma limitação geométrica oriunda da incompressibilidade do líquido e da rigidez do meio poroso. Alguns casos particulares são simulados e os resultados comparados com resultados clássicos, mostrando as consequências de não levar em conta as restrições inerentes ao problema. / This work is concerned with a new mathematical modelling for describing the flow of an incompressible fluid (a liquid) through a rigid, homogeneous and isotropic porous medium, from a Continuum Mixtures point of view. The phenomenon is regarded as the motion of a mixture composed by three overlaping continuous constituents: the first one representing the porous matrix, the second one representing the liquid and the third one representing a (very) low density gas. The proposed mathematical modelling allows a realistic mathematical description for the unsaturated/saturated transition process by means of a combination between a system of partial differential equations and an inequality. This inequality represents a geometrical constraint arising from the liquid incompressibility merged with the porous matrix rigidity. The simulation of some interesting particular cases is carried out presenting a comparison between the obtained results and the classical ones, showing the consequences of disregarding the constraints associated to the phenomenon.