A variational approach to derive a piecewise constant conservative approximation of anisotropic diffusion equations is presented. A priori error estimates are derived assuming usual mesh regularity constraints and a posteriori error indicator is proposed and analyzed for the model problem