Statistical calibration of probabilistic numerical methods for parabolic PDE-constrained inverse problems