Numerical simulation of influence of thermokarst lake horizontal expansion rate on talik development beneath thermokarst lakes on Qinghai-Tibet Plateau