Inclusion sizing, especially for large inclusions greater than 30μm provides important reference for metallurgical process control and fatigue life assessment of steel. Ultrasonic non-destructive testing (NDT) shows great advantages in detecting infrequently occurred large inclusions than eddy current, magnetic particle, microscopic or macroscopic examination procedures. In this paper, the performance of inclusion sizing by immersion ultrasonic C-scan imaging is studied numerically. A two-dimensional model that consists of spherically focused transducer, water couplant and steel with embedded inclusion is established and solved numerically by the finite element method. The signal intensity distributions of inclusion with different sizes are acquired and the effects of inclusion type, shape, orientation on signal intensity distribution are analysed. The results show that the 6dB-drop threshold has the smallest relative error compared with the 12dB-drop threshold and the full-drop threshold, which is better for determining inclusion size larger than 100μm. Experiment is also performed to validate the simulated results.