Multi-cluster horizontal well fracturing is one of the key technologies to develop the unconventional reservoirs such as shales. However, the field data shows that some perforation clusters have little production contribution. In this study, a three-dimensional (3D) numerical model for simulating the multiple fracture propagation based on 3D block discrete element method was established, and this model considers the stress interference, perforation friction and fluid-mechanical coupling effect. In order to determine the most appropriate measures to improve the uniformity of multiple fracture propagation, the effect of the geologic and engineering parameters on the multiple fracture propagation in shale reservoirs is investigated. The modeling results show that the geometry of each fracture within a stage is different, and the outer fractures generally receive more fracturing fluid than the interior fractures. The vertical stress almost has no effect on the geometries of multiple fractures. However, higher horizontal stress difference, larger cluster spacing, smaller perforation number, higher injection rate, and smaller fracturing fluid viscosity are conducive to promote the uniform propagation of multiple fractures. The existence of bedding planes will increase the fluid filtration, resulting in a reduction in fracture length. The middle two fractures receive less fluid and the width of them is smaller. Through analyzing the numerical results, a large amount of fracturing fluid should be injected and the proppant with smaller size is suggested to be used to effectively prop the bedding planes. Cluster spacing and perforation number should be controlled in an appropriate range according to reservoir properties. Increasing the injection rate and reducing the viscosity of fracturing fluid are important means to improve the geometry of each fracture.