Thermal pairing in Richardson model