A class of linear parabolic equations is considered. We derive a framework for the a posteriori error analysis of time discretisations by Richardson extrapolation of arbitrary order combined with finite element discretisations in space. We use the idea of elliptic reconstructions and certain bounds for the Green's function of the parabolic operator. The crucial point in the analysis is the design of suitable polynomial reconstructions in time from approximations that are given only in the mesh points.