Numerical study of the thermal characteristics of a shallow tunnel section with a two-phase closed thermosyphon group in a permafrost region under climate warming