Abstract Mining sequence is one of the important research contents for mining engineering. Different secondary stress field would be formed with different mining process and filling sequence. By means of reasonable arrangement and adjustment of the mining sequence of orebody, the stress distribution of the surrounding rock is able to be remarkably changed, the current excavated working face is able to be in the stress-reduced area and the phenomenon of stress concentration in course of mining process could be also avoided. Rational mining sequence is able to largely improve the stability of stope and the security of working face. According to multilayer orebody occurrence characteristics and mining technical conditions of Laochang 13-8^# ore groups, by using 3Dmine-MIDAS/GTS-FLAC^3D coupled numerical simulation method, numerical simulation study on the mining sequence of multilayer orebody is carried out and four kinds of optimization schemes of multilayer orebody mining sequence are put forward, namely, downward mining from upper to lower, overhead mining from upper to lower, downward mining from lower to upper, and overhead mining from lower to upper. By applying the multi-index comprehensive evaluating model based on index contentment degree and considering stress fields, displacement fields, and crack fields comprehensively, the four multi-index comprehensive contentment degrees of the optimization schemes for mining sequence are 0.684, 0.641, 0.653, and 0.683, respectively. The contentment degree of the overhead mining from lower to upper is 0.683. Therefore, the result indicates that the overhead mining from lower to upper is feasible in the theory.