Bonferroni mean (BM) is a very useful aggregation operator, which can consider the correlations between the aggregated arguments and the multi-valued neutrosophic set can be much more convenient to denote the incomplete, indeterminate and inconsistent information, in this paper, we applied the Bonferroni mean to the multi-valued neutrosophic set, and proposed some Bonferroni mean operators of multi-valued neutrosophic numbers (MVNNs). First, we gave some operational laws and a comparison method of MVNNs, then we presented the weighted Bonferroni mean (WBM) operator and weighted geometric Bonferroni mean (WGBM) operator. Further, we proposed the multi-valued neutrosophic weighted Bonferroni mean (MVNWBM) operator and the multi-valued neutrosophic weighted geometric Bonferroni mean (MVNWGBM) operator and some properties of them are also investigated. Finally, the decision making methods are developed based on MVNWGBM operator and MVNWBM operator, and an example about investment selection is given to illustrate the applications of the developed methods and the influence of different parameter values on the decision-making results.