In general, borehole stability analysis and borehole response analysis are carried out considering that both hydraulic and mechanical parameters of the rock mass are deterministic. However, it is a well known fact, that rock masses and in particular sedimentary rock masses may show a considerable degree of heterogeneity, in micro, meso and macro scale. These heterogeneities produce spatial variability in mechanical and hydraulic properties of the rock medium. This variability can be very pronounced. The present work proposes the development of numerical analysis procedures, using finite elements, in order to analyze single and two phases flow, coupled fluid mechanical processes that take into account the spatial variability of hydraulic and mechanical properties and the variability of the initial stresses and pore pressures. In this study, the developed numerical procedures are used in two distinct phases. In the borehole stochastic stability analysis and in the borehole stochastic response analysis during the production, considering fluid mechanical coupling and two phase flow.