In surface mining operations, the ease with which a muckpile may be dug is a dominant factor in mining efficiency. Within broad limits, this is largely dependent on muckpile shape and looseness within a muckpile._x000D_
Moreover, prediction of grade boundary movement and subsequent dilu- tion is of interest to most open cut mines. Such a predictive tool would enable blast design changes to be made to reduce mixing of certain grade regions. At the very least, a good knowledge of post-blast region bound- aries can reduce layout of the grade boundaries and decrease grade dis- crepancies between pit and plant. The model presented in this paper is able to consider various options for a production blast and predict consequent muckpile shapes and grade boundary movements in a three dimensional manner. The overburden geometry, hole patterns, hole inclination and charge distribution are all explicit model inputs. The model divides the rock mass into blocks and calculates a velocity for each block. Three controlling parameters are used in the model, which are related to the explosive energy per unit length of blast hole, energy attenuation of the rock mass and inter-row delay, respectively. The output of the model is a 3 dimensional muckpile shape (contour or any cross section) and a contour map of grade distribution within the muckpile. The model parameters are site specific, since their values depend on the rock condition and explosive-rock interaction which varies from site to site. In order to "tune" the model to a particular site, parameter values must be found which reproduce an observed muckpile profile shape from a given blast design. This calibration of the model requires at least one surveyed muckpile profile shape from a production blast. The model parameters which best reproduce that shape are found by a pattern search technique, and may then be used to predict the effect of changes in blast design on muckpile profile shape and grade boundary movement.