Scaling relations and mass bias in hydrodynamical f (R) gravity simulations of galaxy clusters