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.