Well-posedness and variational numerical scheme for an adaptive model in highly heterogeneous porous media