In this paper, we are interested in the modelling and the finite element approximation of a petroleum reservoir, in axisymmetric form. The flow in the porous medium is governed by the Darcy-Forchheimer equation coupled with a rather exhaustive energy equation. The semi-discretized problem is put under a mixed variational formulation, whose approximation is achieved by means of conservative Raviart-Thomas elements for the fluxes and of piecewise constant elements for the pressure and the temperature. The discrete problem thus obtained is well-posed and a posteriori error estimates are also established. Numerical tests are presented validating the developed code.