Ing. Joachim Werther on the occasion of his 80th birthdayOptimization of ironmaking blast furnaces (BFs) must consider bottom and top operations. This is done by a recently developed integrated BF model covering the BF below the burden surface. Key operating parameters are focused, viz. coal rate, oxygen enrichment, blast temperature, and burden central opening extent. Numerical orthogonal experiments are thus conducted. Fitting the simulation data yields a set of empirical correlations predicting BF performance indicators. After the fitting accuracy test, the correlations are used for multi-objective optimization and operatable zone identification.