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.

Submitted by Johnny Rodrigues (johnnyrodrigues@ufcg.edu.br) on 2018-07-24T14:45:11Z No. of bitstreams: 1 LUCIANO MARTINS BARROS - DISSERTAÇÃO PPGMAT 2010..pdf: 2620845 bytes, checksum: dc95731aa66d4dab123e738a1dc6b49c (MD5) / Made available in DSpace on 2018-07-24T14:45:11Z (GMT). No. of bitstreams: 1 LUCIANO MARTINS BARROS - DISSERTAÇÃO PPGMAT 2010..pdf: 2620845 bytes, checksum: dc95731aa66d4dab123e738a1dc6b49c (MD5) Previous issue date: 2010-03 / Neste trabalho obtivemos uma solução do problema de Riemann associado a um sistema de duas leis de conservação proveniente da modelagem matemática de um escoamento trifásico num meio poroso. Consideramos o caso de um reservatório petrolífero contendo inicialmente uma mistura arbitrária do tipo gás/óleo à ser deslocada pela injeção de uma mistura do tipo água/gás, também arbitrária. Usando uma combinação de métodos analíticos e computacionais determinamos a geometria das chamadas curvas de onda sob a condição de entropia de viscosidade, com matriz de viscosidade sendo a identidade. Determinamos todas as possíveis sequências de ondas que descrevem o escoamento para cada par de misturas de injeção e de produção representandoosdadosdeRiemann. Mostramosqueparadadosdeproduçãorepresentando uma mistura próxima de óleo puro, ou de gás puro, apenas duas ondas estão presentes no escoamento, independentemente da mistura de injeção. No entanto, para dados de produção representando uma melhor proporção gás/óleo mostramos a existência de uma faixa de dados de injeção para a qual três ondas estão presentes no escoamento, uma delas sendo uma onda de choque transicional. / In this work we describe a Riemann solution for a system of two conservation laws modeling a three-phase flow in a porous media. We consider the case where a petroleum reservoir is initially filled with an arbitrary gas/oil mixture to be displaced by the injection of a gas/water mixture, also arbitrary. By using a combination of analytical and computational methods we obtain the geometry of the so called wave curves under the viscous profile entropy condition, with the viscosity matrix as the identity. We determine all wave sequences describing the flow for each pair of injection andproductionmixtures,representingtheRiemanndata. Weshowthatforproduction mixture data close to pure oil, or pure gas, only two waves are present in the flow independentlyontheinjectionmixture. Nevertheless, forproductiondatarepresenting a more proportional gas/oil mixture we show the existence of a injection data range for which three waves are present in the flow, one of them being a transitional shock wave.

Le but de cette thèse est de faire l'étude de méthodes de volumes finis pour des équations aux dérivées partielles déterministes et stochastiques; nous effectuons des simulations numériques et démontrons des résultats de convergence d'algorithmes.Au Chapitre 1, nous appliquons un schéma semi-implicite en temps combiné avec la méthode de volumes finis généralisés SUSHI pour la simulation d'écoulements à densité variable en milieu poreux; il vient à résoudre une équation de convection-diffusion parabolique pour la concentration couplée à une équation elliptique en pression. Nous présentons ensuite une méthode de simulation numérique pour un problème d'écoulements à densité variable couplé à un transfert de chaleur.Au Chapitre 2, nous effectuons une étude numérique de l'équation de Burgers non visqueuse en dimension un d'espace, avec des conditions aux limites périodiques, un terme source stochastique de moyenne spatiale nulle et une condition initiale déterministe. Nous utilisons un schéma de volumes finis combinant une intégration en temps de type Euler-Maruyama avec le flux numérique de Godunov. Nous effectuons des simulations par la méthode de Monte-Carlo et analysons les résultats pour différentes régularités du terme source. Il apparaît que la moyenne empirique des réalisations converge vers la moyenne en espace de la condition initiale déterministe quand t → ∞. Par ailleurs, la variance empirique converge elle aussi en temps long, vers une valeur qui dépend de la régularité et de l'amplitude du terme stochastique.Au Chapitre 3, nous démontrons la convergence d'une méthode de volumes finis pour une loi de conservation du premier ordre avec une fonction de flux monotone et un terme source multiplicatif faisant intervenir un processus Q-Wiener. Le terme de convection est discrétisé à l'aide d'un schéma amont. Nous présentons des estimations a priori pour la solution discrète dont en particulier une estimation de type BV faible. A l'aide d'une interpolation en temps, nous démontrons deux inégalité entropiques vérifiées par la solution discrète, ce qui nous permet de prouver que la solution discrète converge selon une sous-suite vers une solution stochastique faible entropique à valeurs mesures de la loi de conservation.Au Chapitre 4, nous obtenons des résultats similaires à ceux du Chapitre 3 dans le cas où la fonction flux n'est pas monotone; le terme de convection est discrétisé à l'aide d'un schéma monotone. / This thesis bears on numerical methods for deterministic and stochastic partial differential equations; we perform numerical simulations by means of finite volume methods and prove convergence results.In Chapter 1, we apply a semi-implicit time scheme together with the generalized finite volume method SUSHI for the numerical simulation of density driven flows in porous media; it amounts to solve a nonlinear convection-diffusion parabolic equation for the concentration coupled with an elliptic equation for the pressure. We then propose a numerical scheme to simulate density driven flows in porous media coupled to heat transfer. We use adaptive meshes, based upon square or cubic volume elements.In Chapter 2, We perform Monte-Carlo simulations in the one-dimensional torus for the first order Burgers equation forced by a stochastic source term with zero spatial integral. We suppose that this source term is a white noise in time, and consider various regularities in space. We apply a finite volume scheme combining the Godunov numerical flux with the Euler-Maruyama integrator in time. It turns out that the empirical mean converges to the space-average of the deterministic initial condition as t → ∞. The empirical variance also stabilizes for large time, towards a limit which depends on the space regularity and on the intensity of the noise.In Chapter 3, we study a time explicit finite volume method with an upwind scheme for a first order conservation law with a monotone flux function and a multiplicative source term involving a Q-Wiener process. We present some a priori estimates including a weak BV estimate. After performing a time interpolation, we prove two entropy inequalities for the discrete solution and show that it converges up to a subsequence to a stochastic measure-valued entropy solution of the conservation law in the sense of Young measures.In Chapter 4, we obtain similar results as in Chapter 3, in the case that the flux function is non-monotone, and that the convection term is discretized by means of a monotone scheme.

[pt] Tradicionalmente, os testes de formação em poços de petróleo buscam caracterizar o campo de permeabilidades a partir da interpretação dos transientes de pressão (PTA) nos períodos de fluxo e estática baseados em modelos isotérmicos de escoamento em meios porosos. Com o avanço da instrumentação dos testes, registros mais precisos de temperatura passaram a estar disponíveis e fomentaram a pesquisa baseada em modelos não isotérmicos que possibilitaram a análise a partir dos transientes de temperatura (TTA). Além da caracterização de parâmetros do reservatório como permeabilidade e porosidade com a interpretação dos transientes de temperatura, os dados de pressão obtidos a partir de um modelo não isotérmico representa de forma mais fidedigna o fenômeno físico sobretudo quando os testes são submetidos a maiores diferenciais de pressão. Este trabalho consiste no desenvolvimento de um simulador para teste de formação que considera a modelagem não isotérmica de reservatório unidimensional radial acoplado a um poço produtor e na utilização deste simulador, associado a métodos de otimização multivariável, para resolução do problema inverso da caracterização de parâmetros do reservatório. Alguns métodos de otimização foram testados e o algoritmo do Simplex de Nelder-Mead apresentou melhor eficácia. Foram estabelecidos três tipos de problemas e utilizados em três casos hipotéticos considerando inclusive a imposição artificial de ruídos nos sinais de pressão e temperatura utilizados para resolução do problema inverso. / [en] Traditionally, oil well formation tests aim to characterize the reservoir permeability field from pressure transient analysis (PTA) of drawdown and build up based on isothermal flow models in porous media. With the advancement of well test instrumentation, more accurate temperature records became available and have encouraged researches based on non-isothermal models that made possible the temperature transient analysis (TTA). In addition to the characterization of reservoir parameters such as permeability and porosity by TTA, the pressure data obtained from a non-isothermal model represent better the physical phenomenon, especially when the tests are subjected to greater drawdowns. This work consists in the development of a simulator for formation test that considers non-isothermal modeling of a unidimensional radial reservoir coupled to a production well and in the use of this simulator, associated with multivariable optimization methods, to solve the inverse problem of reservoir parameters characterization. Some optimization methods were tested and the Nelder-Mead Simplex algorithm presented better efficiency. Three types of problems were established and used in three hypothetical cases, including artificially imposed noise in pressure and temperature signals used to solve the inverse problem.

[pt] A teoria da homogeneização de equações diferenciais tornou-se um campo aberto de pesquisa em diversas áreas das ciências exatas e mostrou-se ser uma poderosa ferramenta para a compreensão do comportamento global de materiais heterogêneos. Apesar de ser conhecido que a dedução da lei de Darcy através das equações de Navier-Stokes já é tema debatido há décadas muitas questões continuam em aberto, principalmente a respeito de condições de contorno mais complexas, casos envolvendo fluxos multifásicos e técnicas de homogeneização numérica. Sabe-se que a lei de Darcy se apresenta sob forma de uma relação linear apenas para um intervalo de gradiente hidráulico e que este intervalo se sobrepõe ao intervalo de fluxo laminar do fluido através dos vazios do solo.Se propõe neste trabalho, então, compreender a perda de linearidade na lei de Darcy, a partir da teoria da homogeneização, modificando e explorando os resultados obtidos anteriormente na literatura. / [en] The theory of homogenization of differential equations has become an open field of research in several areas of the exact sciences and has proved to be a powerful tool for understanding the global behavior of heterogeneous materials. Despite knowing that the deduction of Darcy s law through the Navier-Stokes equations has been debated for decades, many questions remain open, mainly regarding more complex boundary cenditions, cases involving multiphase flows and the numerical homogenization techniques. It is known that Darcy s law is presented in the form of a linear relationship only for a range of hydraulic gradient that overlaps the range of laminar flow of fluids through soil voids. Therefore, it is proposed in this work to understand the loss of linearity in Darcy s law, based on the theory of homogenization, modifying and exploring the limit results obtained in the literature.

Cette thèse porte sur la modélisation de l’écoulement et du transport en milieu poreux ;nous effectuons des simulations numériques et démontrons des résultats de convergence d’algorithmes.Au Chapitre 1, nous appliquons des méthodes de volumes finis pour la simulation d’écoulements à densité variable en milieu poreux ; il vient à résoudre une équation de convection diffusion parabolique pour la concentration couplée à une équation elliptique en pression.Nous nous appuyons sur la méthode des volumes finis standard pour le calcul des solutions de deux problèmes spécifiques : une interface en rotation entre eau salée et eau douce et le problème de Henry. Nous appliquons ensuite la méthode de volumes finis généralisés SUSHI pour la simulation des mêmes problèmes ainsi que celle d’un problème de bassin salé en dimension trois d’espace. Nous nous appuyons sur des maillages adaptatifs, basés sur des éléments de volume carrés ou cubiques.Au Chapitre 2, nous nous appuyons de nouveau sur la méthode de volumes finis généralisés SUSHI pour la discrétisation de l’équation de Richards, une équation elliptique parabolique pour le calcul d’écoulements en milieu poreux. Le terme de diffusion peut être anisotrope et hétérogène. Cette classe de méthodes localement conservatrices s’applique àune grande variété de mailles polyédriques non structurées qui peuvent ne pas se raccorder.La discrétisation en temps est totalement implicite. Nous obtenons un résultat de convergence basé sur des estimations a priori et sur l’application du théorème de compacité de Fréchet-Kolmogorov. Nous présentons aussi des tests numériques.Au Chapitre 3, nous discrétisons le problème de Signorini par un schéma de type gradient,qui s’écrit à l’aide d’une formulation variationnelle discrète et est basé sur des approximations indépendantes des fonctions et des gradients. On montre l’existence et l’unicité de la solution discrète ainsi que sa convergence vers la solution faible du problème continu. Nous présentons ensuite un schéma numérique basé sur la méthode SUSHI.Au Chapitre 4, nous appliquons un schéma semi-implicite en temps combiné avec la méthode SUSHI pour la résolution numérique d’un problème d’écoulements à densité variable ;il s’agit de résoudre des équations paraboliques de convection-diffusion pour la densité de soluté et le transport de la température ainsi que pour la pression. Nous simulons l’avance d’un front d’eau douce assez chaude et le transport de chaleur dans un aquifère captif qui est initialement chargé d’eau froide salée. Nous utilisons des maillages adaptatifs, basés sur des éléments de volume carrés. / This thesis bears on the modelling of groundwater flow and transport in porous media; we perform numerical simulations by means of finite volume methods and prove convergence results. In Chapter 1, we first apply a semi-implicit standard finite volume method and then the generalized finite volume method SUSHI for the numerical simulation of density driven flows in porous media; we solve a nonlinear convection-diffusion parabolic equation for the concentration coupled with an elliptic equation for the pressure. We apply the standard finite volume method to compute the solutions of a problem involving a rotating interface between salt and fresh water and of Henry's problem. We then apply the SUSHI scheme to the same problems as well as to a three dimensional saltpool problem. We use adaptive meshes, based upon square volume elements in space dimension two and cubic volume elements in space dimension three. In Chapter 2, we apply the generalized finite volume method SUSHI to the discretization of Richards equation, an elliptic-parabolic equation modeling groundwater flow, where the diffusion term can be anisotropic and heterogeneous. This class of locally conservative methods can be applied to a wide range of unstructured possibly non-matching polyhedral meshes in arbitrary space dimension. As is needed for Richards equation, the time discretization is fully implicit. We obtain a convergence result based upon a priori estimates and the application of the Fréchet-Kolmogorov compactness theorem. We implement the scheme and present numerical tests. In Chapter 3, we study a gradient scheme for the Signorini problem. Gradient schemes are nonconforming methods written in discrete variational formulation which are based on independent approximations of the functions and the gradients. We prove the existence and uniqueness of the discrete solution as well as its convergence to the weak solution of the Signorini problem. Finally we introduce a numerical scheme based upon the SUSHI discretization and present numerical results. In Chapter 4, we apply a semi-implicit scheme in time together with a generalized finite volume method for the numerical solution of density driven flows in porous media; it comes to solve nonlinear convection-diffusion parabolic equations for the solute and temperature transport as well as for the pressure. We compute the solutions for a specific problem which describes the advance of a warm fresh water front coupled to heat transfer in a confined aquifer which is initially charged with cold salt water. We use adaptive meshes, based upon square volume elements in space dimension two.

[pt] A Equação da Difusividade Hidráulica Não-Linear (EDHN) modela o escoamento monofásico de fluidos em meios porosos levando em conta a variação das propriedades da rocha e do fluido presente no interior de seus poros. Normalmente, a solução adimensional da linha-fonte pD(rD, tD) para escoamento de líquidos é encontrada por meio do uso da transformada de Laplace ou transformação de Boltzmann, o qual, o perfil transiente de pressões em coordenadas cartesianas é descrito pela função erro complementar erfc(xD, yD, tD) e, em coordenadas cilíndricas pela função integral exponencial Ei(rD, tD). Este trabalho propõe a solução analítica pelo método de expansão assíntotica de primeira ordem em séries, para solução de alguns problemas de escoamento de petróleo em meios porosos com permeabilidade dependente da pressão de poros e termo fonte. A solução geral será implementada no software Matlab (marca registrada) e a calibração do modelo matemático será realizada comparandose a solução obtida neste trabalho com a solução calculada por meio de um simulador de fluxo óleo em meios porosos denominado IMEX (marca registrada) , amplamente usado na indústria de petróleo e em pesquisas científicas e que usa o método de diferenças finitas. A solução geral da equação diferencial é dada pela soma da solução para escoamento de líquidos com permeabilidade constante e o termo de primeira ordem da expansão assintótica, composto pela não linearidade devido à variação de permeabilidade. O efeito da variação instantânea de permeabilidade em função da pressão de poros é claramente demonstrado nos gráficos diagnósticos e especializados apresentados. / [en] The Nonlinear Hydraulic Diffusivity Equation (NHDE) models the singlephase flow of fluids in porous media considering the variation in the properties of the rock and the fluid present inside its pores. Normally, the dimensionless linear solution for the flow of oil is performed using the Laplace and Fourier transform or Boltzmann transformation and provides the unsteady pressure profile in Cartesian coordinates given by complementary error function erfc(xD, yD, tD) and in cylindrical coordinates described by the exponential integral function Ei(rD, tD). This work develops a new analytical model based on an integro-differential solution to predict the formation mechanical damage caused by the permeability loss during the well-reservoir life-cycle for several oil flow problems. The appropriate Green s function (GF) to solve NHDE for each well-reservoir setting approached in this thesis is used. The general solution is implemented in the Matlab (trademark) and the mathematical model calibration will be carried out by comparing the solution obtained in this work to the porous media finite difference oil flow simulator named IMEX (trademark). The general solution of the NHDE is computed by the sum of the linear solution (constant permeability) and the first order term of the asymptotic series expansion, composed of the nonlinear effect of the permeability loss. The instantaneous permeability loss effect is clearly noticed in the diagnostic and specialized plots.