Taking the fluidized pre-reduction process of iron ore powder bubbling fluidized bed as the background, for the problem of non-uniform structure in the bed of gas-solid fluidization process, the non-uniform fluidization characteristics of bicomponent particles are investigated in a cold two-dimensional bubbling fluidized bed by using a combination of physical experiments and mathematical simulations. Fluidization experiments were carried out under typical working conditions by using glass beads to study the effects of apparent gas velocity, mass ratio, and other factors on the non-uniform structure in the bed. Through the experimental observation of the bubble behavior, the effect of the cyclic change in bubble formation, rise and growth to rupture on the bed uniformity were analyzed. The experiments showed that the fluidized bed of two-component particles would be stratified, and the non-uniformity was strong in the upper part and weak in the lower part, and the apparent gas velocity and particle size were the main influencing factors. Based on the Euler-Lagrange reference frame modeling, the fluidization process of the two-dimensional bubble bed was simulated by the CFD-DEM method. The simulations of typical experimental conditions were carried out to further analyze the velocity distribution and the volume ratio of each phase in the bed from the gas-solid interaction level, revealing that the velocity distribution in the upper part of the bed is not uniform, and the gas flow is strongly perturbed, with intense bubble aggregation. The results reveal the reasons for the non-uniform phenomenon of gas-solid fluidization, which can provide a theoretical basis for the regulation of the non-uniform structure of the fluidization process.