A model considering gas outflows due to supernova explosions is developed for the star formation and chemical enrichment for the globular clusters (GCs) in the Milky Way galaxy. Through Monte Carlo simulations, the observed global properties of GCs can be well reproduced, including the metallicity distribution, no-correlation between cluster masses and galactocentric distances, etc. The predicted mass function of the parent clouds for the observed GCs at present day can be well described as a power law with the index of -1.8, which is consistent with the current observations for the molecular clouds.