In this paper, an effective nonstationary phase boundary estimation scheme in electrical impedance tomography (EIT) is presented based on the interacting multiple model (IMM) algorithm. The inverse problem is treated as a stochastic nonlinear state estimation problem with the nonstationary phase boundary (state) being estimated online with the aid of the IMM algorithm. In the design of the IMM algorithm multiple models with different process noise covariances are incorporated to improve estimation performance in spite of the modelling uncertainty. Computer simulations are provided to illustrate the proposed algorithm.