A CFD-PBM (Computational Fluid Dynamic-Population Balance Model) coupled model has been developed to investigate the effects of different number and position of bottom tuyeres and gas flow rate on the bubbly plume flow, inclusion removal and mixing phenomena in gas-stirred ladle. It is found that the dual blowing gives a shorter mixing time and higher inclusion removal ratio in comparison with the center blowing or eccentric blowing with one tuyere. With the increasing of separation angle of two tuyeres, the inclusion removal ratio increases, while mixing time decreases first and then increases. With the increasing of radial position of two tuyeres, the inclusion removal first increases and then decreases, and the mixing time decreases until the radial position exceeds 0.7R from the bottom center, where R is the bottom radius of ladle. It is recommended to use the two tuyeres placed at radial position of 0.6R and the angle of 135 deg in ladle to improve the joint efficiency both the inclusion removal and mixing. With the gas flow rate increasing, the efficiency both mixing and inclusion removal with the optimized tuyeres arrangement increases, however when the gas flow rate exceeds 300 NL/min in 150 ton ladle, the removal ratio and mixing time change little.