Numerical simulation of hydrate slurry flow characteristics in vertical pipes based on population balance theory