Heavy reduction (HR) is a novel technology that could effectively improve the internal porosities and other internal quality problems in continuously cast steel, during which a large reduction deformation is implemented at and after the strand solidification end. In the present paper, non-uniform solidification of the wide-thick slab was calculated with a two-dimensional (2D) heat transfer model. Based on the predicted temperature distribution at the solidification end of the casting strand, a three-dimensional (3D) thermal-mechanical coupled model was developed for investigating the deformation behavior of the internal porosities in wide-thick slab during HR. An Arrhenius-type constitutive model for the studied steel grade was derived based on the measured true stress-strain with single-pass thermosimulation compression experiments and applied to the 3D thermal-mechanical coupled model for improving the calculation accuracy. With the developed 3D thermal-mechanical coupled model, deformation behavior of the two artificial porosities located at the slab center of 1/2 width and 1/8 width during HR was investigated under different condition of HR deformation, HR start position, and HR reduction mode. Based on the calculated porosity closure degree (ηs) and the corresponding equivalent strain (εeq) under different HR conditions, a prediction model that describes the quantitative relationship between ηs and εeq was derived for directly and accurately evaluating the process effect of HR on improving the internal porosities in wide-thick slab.