Abstract: In order to explore the applicability and accuracy of different drag models in the numerical study of dense gas-solid two-phase flow in spouting beds, based on the developed coupled flow program of mesoscale Lattice Boltzmann discrete Element method（LBM-DEM）, three traction models of syamlal-O ’Brien, Gidaspow and EMMS were used to perform a numerical simulation of the micro-spout bed. The results under different drag models were compared and analyzed, and the accuracy of prediction results of each model were verified by experiments. The results show that the three drag models can completely show the whole process of bubble evolution and particle fluidization in a spouting bed. Among them, the EMMS model has the best overall applicability that the overall prediction results are most consistent with the experimental results, which can better show the details of the non-uniform fluidization process in the bed. Gidaspow model has more advantages in bubble evolution of homogeneous flow state. The drag force value calculated by Syamlal-O ’Brien model is always on the high side, so attention should be paid to avoid the possible side-wall effect. The above research results and methods can provide reference for the study of dense fluid-particle two-phase system in a wide range of fields.