We apply the Object Kinetic Monte Carlo method to simulate the hydrogen (H) clustering behaviour with large vacancy (V) clusters in tungsten. The main advantage of this method is to consider the temperature effects, large V clusters, and long-time evolution. It is thus a very useful tool to study the microstructure evolution of defects in irradiated materials covering from atomic scale to mesoscale. Simulation results show that the number of trapped H atoms by V and V clusters decreases with the increase of temperature. This is expected, because the de-trapping ability of H atoms increases when the temperature is increased. With the correction of the zero-point-energy in calculating binding energies, the present results show that the number of H atoms trapped by monovacancy decreases for all studied temperatures. We also found that the trapping ability of V clusters for H decreases with the increase of the sizes of V clusters, independent on temperature. The H concentration can largely affect the clustering behaviour of H m V n. When the H saturation concentration is reached, the clustering behaviour of H m V n is mainly determined by the sizes of V clusters.