The present work focuses on the development of a novel statistical approach for the analysis of the behavior of a 2D fluidized bed operating under bubbling and slugging conditions. Experimental data were collected by means of a purposely built lab-scale 2D fluidized bed and an original digital image analysis technique (DIAT), allowing the acquisition of a number of data sufficient to perform a statistical analysis of bubble dynamics. In particular, measurements of bubble size distributions (BSD) along the bed were performed. The analysis allowed to characterize the shape of BSD for different particle systems and fluidization gas velocities; interestingly a bimodal shape of BSD is observed, likely due to the simultaneous presence of coalesce...