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.

[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 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.