Apresentamos um modelo teórico para explicar as tran sições de valéncia do Sm nas ligas Smi _ x M xS, onde M é um metal de transição (por ex. Y ou La), como uma função da concentração x. O sistema é descrito por um modelo de duas bandas: uma banda de largura nula e de energia E0 (o nível 4f do Sm) e uma banda de condução, a qual é tratada na aproximação do potencial coerente (CPA). Os elétrons 4f interagem entre si via uma repulsão coulombiana finita U, que separa a banda f em duas sub-bandas de energias E0 e Eo+U. Consi dera-se que as bandas f e de condução se hibridizam, sendo V o parâmetro de hibridização. Com este mo delo é possível obter uma transição continua de um estado funda mental não-magnético a um estado não magnético ou magnético, con forme se varia a razão U/V. A variação de valência como uma função da pressão nos monocalcogenetos de Sm (x = O) é também calculada. O modelo usa do é uma extensão do descrito acima, ao qual nOs incluimos uma repulsão coulombiana G entre elétrons f e de condução, mantendo finito o valor de U. Para diferentes valores de G/V obtemos tran sições de primeira ou de segunda ordem de um estado semicondutor a um estado metálico. A resistividade elétrica em ambas as fases é também calculada como uma função da temperatura e os re sultados teoricos obtidos estão em boa concordância com os experimentais. / We present a theoretical model to explain the valence transitions of Sm in Sm1-xMxS alloys, where M is a transition metal (e.g. Y or La) as a function of the concentration x. The system is described by a two-band model: a zero width band of energy E0 (the 4f-level of Sm) and a conduction band, which is treated in the coherent potencial approximation (CPA). The 4felectrons interact between them via a finite intra-atomic coulomb repulsion U, which splits the f-band in two sub-bands at energies E0 and Eo+U. The f -and conduction bands hybridize, being V the hybridization parameter. Within this model it is possible to obtain a continuous valence transition from a non- -magnetic to a non-magnetic or a magnetic ground-state, varying the ratio U/V. The change of valence as a function of pressure in the Sm monochalcogenides (x =0) is also computed. The model used is an extension of the one described above, in which we include a Coulomb repulsion G between f - and conduction electrons, preserving the finite value of U. For different values of G/V we obtain first or second order transitions from a semiconducting to a metallic state. The electrical resistivity in both fases as a function of temperature is also computed and the theoretical results obtained are in good agreement with the experimental ones.

Neste trabalho estudamos o modelo da rede de Kondo em uma rede kagome, buscando uma maior compreensão dos efeitos da frustração geométrica em sistemas de férmions pesados. Para tanto, fizemos uma aproximação de campo médio no hamiltoniano do sistema que serve para todas as fases do sistema. Analisamos inicialmente o caso não magnético. Obtemos neste limite as energias eletrônicas e as funções de Green necessárias ao cálculo numérico autoconsistente das ocupações e do parâmetro de Kondo. Os resultados encontrados estão em concordância qualitativa com trabalhos publicados em outras geometrias. A seguir analisamos o caso magnético, onde introduzimos uma aproximação suplementar, a qual é compatível com a de campo médio já considerada e, em princípio, existente apenas em sistemas com frustração geométrica. Realizamos cálculos autoconsistentes através de somas sobre as frequências de Matsubara. Os resultados mostram que não há coexistência entre ordem magnética e efeito Kondo, além de haver a supressão do antiferromagnetismo com o aumento de temperatura e variações no preenchimento de bandas. / In this work we study the Kondo Lattice model for the kagome lattice, in order to understand better the effects of geometrical frustration in heavy-fermion systems. In this context, we consider a mean field scheme valid for all the system’s phases. Firstly, we analyzed the nonmagnetic case. In this approximation the electron energies and spectral functions are reachable, then we use the density of states to calculate the occupations selfconsistently. Our results are qualitatively compared with previous works in other geometries. In the second part we introduce an approximation for magnestism, which takes into account the mean field scheme considered and the presence of geometrical frustration. Self-consistent calculations are done through the frequencies summation method. Our results show that the magnetism is supressed when the temperature is increased or the band filling deviates from half-filling. Besides, the coexistence of magnetic order and Kondo effect is not observable.

Neste trabalho estudamos modelos teóricos que descrevem sistemas eletrônicos fortemente correlacionados, em especial o modelo t-J, e suas aplicações a compostos de óxidos de cobre, notadamente os compostos que apresentam supercondutividade de alta temperatura crítica e o composto Sr2CuO2Cl2. No primeiro capítulo do trabalho, fazemos uma exposição de três modelos que envolvem o tratamento das interações elétron-elétron, que são os modelos de Hubbard de uma banda, o modelo de Heisenberg e o modelo t-J. Na dedução deste último fazemos uma expansão canônica do hamiltoniano de Hubbard, no limite de acoplamento forte, levando-nos a obter um novo hamiltoniano que pode ser utilizado para descrever um sistema antiferromagnético bidimensional na presen- ça de lacunas, que é exatamente o que caracteriza os compostos supercondutores de alta temperatura crítica na sua fase de baixa dopagem.Após termos obtido o hamiltoniano que descreve o modelo t-J, aplicamos à este uma descrição de polarons de spin, numa representação de holons, que são férmions sem spin, e spinons, que são bósons que carregam somente os graus de liberdade de spin. Utilizando uma função de Green para descrever a propagação do polaron pela rede, obtemos uma equação para a sua autoenergia somando uma série de diagramas de Feynman, sendo que para este cálculo utilizamos a aproxima ção de Born autoconsistente[1]. Do ponto de vista numérico demonstramos que a equação integral de Dyson resultante do tratamento anterior não requer um procedimento iterativo para sua solução, e com isto conseguimos trabalhar com sistemas com grande número de partículas. Os resultados mostram, como um aspecto novo, que o tempo de vida média do holon tem um valor bastante grande no ponto (π,0 ) da rede recíproca, perto da singularidade de Van Hove mencionada na literatura[2]. Este aspecto, e suas implicações, é amplamente discutido neste capítulo. No capítulo 3 estudamos o modelo estendido t-t'-J, com tunelamento à segundos vizinhos e a incorporação dos termos de três sítios[3]. Fazemos a mesma formulação do capítulo anterior, e discutimos as aplicações dos nossos resultados ao óxido mencionado anteriormente. Finalmente, no último capítulo apresentamos uma aplicação original do modelo t-J à uma rede retangular, levemente distorcida, e demonstramos que os resultados do capítulo 3 são reproduzidos sem necessidade de introduzir termos de tunelamento adicionais no hamiltoniano. Esta aplicação pode se tornar relevante para o estudo das fases de tiras encontradas recentemente nesses materiais de óxidos de cobre.

Elétrons em bandas estreitas de energia são fortemente afetados pela interação coulombiana como também por desordem na rede. Ambos os efeitos podem levar à localização, mas de naturezas diferentes: estado isolante de Mott, induzido por correlação, e localização de Anderson, induzida por desordem. A existência da fase de Mott ´e também significativamente dependente do preenchimento da banda. Abordagens teóricas para lidar com esse tipo de sistema são usualmente baseadas no hamiltoniano de Hubbard ou modelos relacionados, incluindo desordem como uma distribuição de energias locais. Neste trabalho, utilizando a Teoria de Campo Médio Dinâmico (DMFT), estudamos o modelo de Anderson-Falicov-Kimball e sua versão de três bandas, obtida como uma simplificação do modelo de Hubbard de três bandas associado aos planos de CuO2 dos cupratos supercondutores de alta temperatura crítica, realizando nossa análise para os casos magnético e não magnético. A densidade de estados de uma partícula é obtida por medias aritmética e geométrica sobre a desordem, já que somente a última pode detectar a localização na ausência de um gap de energia. Variando as intensidades de interação coulombiana e desordem, construímos diagramas de fases para esse modelo, onde identificamos transições metal-isolante mediadas por correlação e desordem, bem como a inter-relação entre esses efeitos. Isso é feito para vários preenchimentos de banda, já que nosso principal interesse aqui é estudar como a variação da densidade de elétrons (dopagem) afeta os diagramas de fases previamente obtidos na ausência de dopagem. Para o modelo de uma banda no caso paramagnético, as informações reveladas pela densidade de estados são confirmadas pela análise das condutividades estática e dinâmica, incluindo efeitos de temperatura. Quando consideramos a solução magnética, observamos o comportamento da temperatura de Néel e podemos apresentar um diagrama de fases mais completo. Além de uma análise bastante extensa do modelo de uma banda, fazemos um estudo inicial do modelo de três bandas, focalizando comportamentos que possam vir a ser comparados ao que se observa nos óxidos supercondutores. / Electrons in narrow-band solids are strongly affected by the Coulomb interaction as well as lattice disorder. Both effects can lead to localization, but of different nature: correlation-induced Mott insulating state, and disorder-induced Anderson localization. The existence of the Mott phase is also significantly dependent on the band filling. Theoretical approaches to deal with this kind of system are usually based on the Hubbard Hamiltonian or related models, including disorder as a distribution of the on-site energies. In this work, utilizing Dynamic Mean Field Theory (DMFT), we study the Anderson- Falicov-Kimball and its three-band version, obtained as a simplification of the three-band Hubbard model associated to the CuO2 planes of high-critical-temperature cuprate superconductors, performing our analysis for the magnetic and non-magnetic cases. The one-particle density of states is obtained by both arithmetic and geometrical averages over disorder, since only the latter can detect localization in the absence of an energy gap. Varying the strengths of Coulomb interaction and disorder, we construct phase diagrams for these models, where we identify metal-insulator transitions driven by correlation and disorder, as well as the interplay between these effects. This is done for various band fillings, since our main interest here is to study how the variation of the electron density affects the phase diagrams previously obtained in the absence of doping. For the one-band model in the paramagnetic case, the picture revealed by the density of states is further checked by evaluating the static and dynamic conductivities, including temperature effects. When we consider the magnetic solution, we observe the N´eel temperature behavior, and we are able to present a more complete phase diagram. Besides a quite extensive analysis of the one-band model, we develop an initial study of the threeband model, focusing on behaviors that might be linked with what is observed on the superconducting oxides.

[pt] A modelagem de problemas relacionados a geomecânica do reservatório é tradicionalmente realizada por elementos finitos. Para utilizar esse método é preciso que o modelo englobe uma região consideravelmente superior a região em que o reservatório está inserido, além de necessitar imposição condições de contorno. Pensando em reduzir a necessidade de discretização de grandes regiões do maciço rochoso é proposto o método das funções de Green para análise geomecânica. Este método é baseado no uso de soluções analíticas clássicas (solução fundamental de Kelvin, solução fundamental de Melan, por exemplo) como soluções auxiliares para resolver problemas elasticamente heterogêneo e não lineares em meios saturados de fluidos. A não linearidade do material pode ser devido a deformações irreversíveis ou resposta de elasticidade não linear típica da análise 4D. O procedimento de solução geral depende de um método de colocação discreta e uma abordagem iterativa de ponto fixo para construir o campo de deslocamento. Esse método teve sua convergência verificada através de modelos simplificados que possuem solução analítica. Visando o avanço do desempenho computacional do método das funções de Green, foram feitas duas modificações independentes utilizando inteligência artificial. A primeira modificação é baseada na integração de dois conceitos principais: o teorema da reciprocidade e a capacidade de generalização das redes neurais artificiais. O teorema da reciprocidade é usado para formular a expressão matemática que rege o problema geomecânico, que é então discretizado no espaço em elementos inteligentes. O comportamento do campo de deformação dentro desses novos elementos é previsto usando uma rede neural artificial. Para fazer essas previsões, a rede neural usa condições de contorno de deslocamento, propriedades do material e a forma geométrica do elemento como dados de entrada. A segunda modificação consiste na utilização de soluções auxiliares que considerem a heterogeneidade de maciços estratificados. Essas soluções são obtidas através do treinamento de redes neurais artificiais que tem como dado de saída o deslocamento em um determinado ponto do maciço estratificado devido a aplicação de uma força pontual em um ponto no interior desse maciço. Para isso, as redes neurais de deslocamentos necessitam das propriedades elásticas e da espessura de cada camada do maciço bem como das coordenadas de aplicação da força pontual e do ponto onde será avaliado o deslocamento. Ao se utilizar essas soluções fundamentais baseadas em inteligência artificial é possível se obter todo o campo de deslocamentos de um problema heterogêneo e elástico de geomecânica do reservatório bastando apenas discretizar o reservatório. Cada uma das modificações do método da função de Green foi avaliada individualmente e observou-se um ganho de pelo menos 5 vezes no tempo de processo, utilizando o mesmo recurso computacional, quando se compara ao método clássico da função de Green. / [en] The analysis and simulation of problems associated with reservoir geomechanics are traditionally performed using the finite element method. However, to perform this analysis, it is necessary to consider a region much larger than the region in which the reservoir is inserted. This is done so that boundary conditions can be applied in an attempt to mimic the effect of the infinite media surrounding the reservoir. With the aim of reducing the need for discretization of large regions of the massif, a Green s functions approach was proposed for reservoir geomechanical analysis. This method is based on the use of classical analytical solutions (Kelvin s fundamental solution, Melan s fundamental solution, for example) as auxiliary solutions to solve elastically heterogeneous and nonlinear problems in fluid-saturated media. The non-linearity of the material can be due to irreversible deformations or non-linear elasticity response typical of 4D analysis. The general solution procedure relies on a discrete collocation method and an iterative fixed-point approach to build the displacement field. This method´s convergence was verified through simplified models that have analytical solutions. As the reduction in processing time is crucial for decision-makers to act in field applications, two improvements were proposed using artificial intelligence (AI) to reduce processing time of the Green s function approach. The first improvement is based on the generalization ability of artificial neural networks (ANN). Due to this characteristic, it was proposed to discretize the model with a coarse mesh of intelligent elements instead of refined mesh of traditional elements based on polynomials. The behavior of the strain field within these new elements is predicted using an ANN. To make these predictions, the neural network uses displacement boundary conditions, material properties and the geometric shape of the element as input data. The examples comparing the intelligent element approach and the traditional method were performed on a computer with 12 threads of 2,6GHz and 32GB RAM. This comparison showed reductions between five and ten times in CPU time, while maintaining the accuracy of the results. The second improvement consists in the use of auxiliary solutions that consider the heterogeneity of stratified massifs. These solutions are obtained through the training of artificial neural networks that have as output the displacement in a certain point of the stratified massif due to the application of a point load inside the massif. This ANN uses as input data elastic properties and the thickness of each layer of the massif, and of the semi-infinite media, as well as the coordinates of the point load and of the point where the displacement is to be evaluated. The use of the developed ANN-based Green’s function approach only demands the discretization of the reservoir itself, thus avoiding the discretization of other regions of the massif. Furthermore, it is possible to obtain the displacement at any point of the massif due to a pore pressure variation within the reservoir without having to solve for the other points in the massif. These two characteristics increase the efficient of the method in relation to traditional methods, such as the finite element method. A numerical example was performed on a computer with 12 threads of 2,6GHz and 32GB RAM to compare the ANN-based Green’s function approach with the traditional approach. The CPU time to obtain the solution using the ANN-based Green’s function approach was five times smaller than the that required by the traditional approach.

[pt] O presente trabalho propõe um método estatístico computacionalmente eficiente (chamado Green-FOSM) para consideração de incertezas em problemas geomecânicos, com o objetivo de melhorar o processo de tomada de decisão ao analisar problemas associados com o processo de injeção ou depleção de fluídos. A novidade do método proposto está associada com a utilização do método das funções de Green (GFA), que, com o auxílio do método estatístico first-order second-moment (FOSM), é utilizado para propagar as inerentes incertezas associadas às propriedades mecânicas do material para o campo de deslocamento da formação geológica. Além disso, através dos conceitos de grid estocástico e função de autocorrelação, o método proposto permite a consideração da variabilidade espacial de variáveis aleatórias de entrada que representam essas propriedades mecânicas. O GFA utiliza as soluções fundamentais da mecânica clássica (solução fundamental de Kelvin, solução fundamental de Melan, entre outras) e o teorema da reciprocidade para determinar o campo de deslocamento de uma formação geológica com geometria irregular e diferentes tipos de materiais. A grande vantagem deste método em relação ao clássico método dos elementos finitos (MEF) é que ele não requer a imposição de condições de contorno e a análise do problema pode ser realizada considerando apenas o domínio do reservatório ou outras regiões de interesse. Esta estratégia de modelagem diminui os graus de liberdade do modelo e o tempo de processamento da análise. Desta forma, como o GFA requer menos esforço computacional, este método torna-se ideal para ser utilizado na propagação de incertezas em problemas geomecânicos. Inicialmente, baseado no método das funções de Green original proposto por Peres et al. (2021), foi proposto uma versão iterativa do método Green-FOSM, que apresenta resultados estatísticos semelhantes aos encontrados através da clássica simulação de Monte Carlo (SMC). Nesta versão original, o campo de deslocamento é PUC-Rio - Certificação Digital Nº 1912634/CA calculado usando um esquema numérico iterativo que diminui o desempenho computacional do método e pode gerar problemas de convergência. Tais limitações tem dificultado a aplicação do GFA original e do método Green-FOSM iterativo em problemas reais. Assim, o presente trabalho desenvolveu uma nova versão do GFA que utiliza um esquema numérico não-iterativo. Para os problemas de validação analisados, o método não-iterativo demonstra ser até 17.5 vezes mais rápido do que a versão original. Além disso, esta versão demonstra ser capaz de expandir a aplicabilidade do GFA, pois os problemas de convergência foram eliminados e os resultados obtidos por este método, ao analisar um perfil geológico representativo do pré-sal brasileiro, são semelhantes aos encontrados via MEF. Por fim, a partir do GFA não-iterativo foi proposta uma versão não-iterativa do método Green-FOSM. Esta versão não-iterativa é capaz de analisar probabilisticamente formações geológicas complexas, como é o caso das formações geológicas do présal brasileiro. Utilizando os mesmos recursos computacionais, o método GreenFOSM não-iterativo é no mínimo 200 vezes mais rápido que o método iterativo. De forma geral, os resultados encontrados nas análises realizadas (determinísticas e probabilísticas) são próximos dos resultados obtidos pelo método de referência (MEF e SMC, respectivamente). / [en] The present work proposes a computationally efficient stochastic statistical method (called Green-FOSM) that considers uncertainties in geomechanical problems, with the objective of improving the decision-making process related to problems associated with the process of fluid injection or depletion. The novelty of the method lies in the use of the Green s function approach (GFA), which, together, with the first-order second-moment statistical method (FOSM), is used to propagate uncertainties associated with the mechanical properties of material to the displacement field of the geological formation. Furthermore, using the concepts of stochastic grid and autocorrelation function, the proposed method allows the consideration of the spatial variability of random variables that represent these mechanical properties. The GFA uses the fundamental solutions of classical mechanics (Kelvin fundamental solution, Melan fundamental solution, among others) and the reciprocity theorem to calculate the displacement field of a geological formation with irregular geometry, and different types of materials. The great advantage of this method compared to the classical finite element method (FEM) is that it does not require the imposition of boundary conditions and the analysis of the problem can be performed considering only the reservoir or other regions of interest. This modeling strategy decreases the degrees of freedom of the model and the CPU time of the deterministic analysis. In this way, as the GFA requires less computational effort, this approach becomes ideal for propagating the uncertainties in geomechanical problems. Initially, an iterative version of the Green-FOSM method was proposed, which presents statistical results similar to those found through the classic Monte Carlo simulation (MCS). In this initial version, the displacement field is calculated using an iterative numerical scheme, which decreases the computational performance of the method and can generate convergence problems. Such limitations would restrict the application of the PUC-Rio - Certificação Digital number 1912634/CA original GFA and the iterative Green-FOSM method in real problems. Thus, the present work also developed a new version of the GFA, which uses a non-iterative numerical scheme. For the proposed validation problems, the non-iterative method proved to be up to 17.5 times faster than the original version. This version is able to expand the applicability of the GFA, since the convergence problems were eliminated and the results obtained by this method, when analyzing a representative geological profile of the Brazilian pre-salt, are similar to those found via FEM. Finally, based on the non-iterative GFA, a non-iterative version of the Green-FOSM method was proposed. This non-iterative version is capable of probabilistically analyzing complex geological formations, such as the Brazilian pre-salt geological formations. Using the same computational resources, the non-iterative GreenFOSM method is at least 200 times faster than the iterative Green-FOSM method. In general, the results found in the investigated analyzes (deterministic and probabilistic) are close to the results obtained by the reference method (FEM and MCS, respectively).

[pt] Modelos analíticos que descrevam o comportamento da pressão são de extrema utilidade na área de avaliação de formações e caracterização de reservatório, pois eles fornecem estimativas sobre diversos parâmetros do reservatório. Este trabalho tem dois objetivos principais: primeiro, estender a solução existente para testes de injetividade e falloff em reservatórios com uma camada e poços horizontais de modo que ela possa ser aplicada também em formações multicamadas com poços horizontais multirramificados. Além disso, este trabalho aplica funções impulso para obter uma formulação alternativa para testes de injetividade em reservatórios estratificados com poços verticais e formações com uma camada e poços horizontais. / [en] Analytical models that describe the pressure behavior are extremely useful for pressure transient analysis and reservoir characterization as they provide estimates of reservoir parameters. This work has two main goals: first, to extend the existing solutions for injectivity/falloff tests in single-layer formations with horizontal wells so that they can be applied to multilayer stratified reservoirs with multilateral horizontal wells. Furthermore, this work applies impulse functions to obtain an alternative formulation for injectivity tests in multilayer commingled formations with vertical wells and single-layer reservoirs with horizontal wells.

[pt] Teste de injetividade é uma técnica convencional em engenharia de reservatórios, utilizada para a recuperação de óleo em reservatórios e avaliação de formações. Geralmente utiliza-se água como fluido injetado, que resulta em um deslocamento do óleo presente devido ao aumento da pressão nos poros. Durante o teste, a resposta de pressão medida fornece diversas informações sobre os parâmetros do reservatório, tal como dados de permeabilidade. Desta forma, pesquisadores têm se dedicado em encontrar equações matemáticas que modelam a resposta de pressão desses testes com objetivo de gerenciamento e manutenção preditiva do reservatório. Neste trabalho, apresentamos uma nova solução analítica para a análise de testes de injetividade, que combina a técnica de convolução pressão-pressão com um modelo radial composto de duas zonas. Essa solução permite avaliar o teste de injetividade mesmo na ausência de dados precisos de vazão, uma vez que a convolução pressão-pressão utiliza exclusivamente os dados de pressão adquiridos em diferentes posições do reservatório. O modelo considerado consiste em dois poços, um injetor, localizado na zona interna do reservatório, e um observador, na zona externa. A validação da solução proposta foi realizada por meio da comparação dos resultados analíticos com aqueles obtidos em um simulador comercial baseado em diferenças finitas. / [en] The injectivity test is a conventional technique in reservoir engineering used for oil recovery and formation evaluation. Typically, water is injected to displace the existing oil by increasing the pressure in the pores. In this test, the pressure response measurement provides valuable information about the reservoir parameters, including permeability data. Therefore, researchers aim to develop mathematical equations that could accurately model pressure response during these tests for reservoir management and maintenance prediction purposes. This work introduces a new analytical solution for injectivity test analysis. The solution combines the pressure-pressure convolution technique with a two-zone radial model. It allows the evaluation of the injectivity test without precise flow rate data, as the pressure-pressure convolution exclusively uses the pressure data acquired at different positions in the reservoir. The reservoir model comprises an injector well in the inner zone of the reservoir and an observation well in the outer zone for measuring pressure response. The proposed solution was validated by comparing the analytical results with those obtained from a finite differences-based commercial simulator.

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