The technique of resonance enhanced multiphoton ionization (REMPI) spectroscopy has already become an important method for studying the energy level structure of atoms and molecules. Analytic expression of 1+2+1 double REMPI probability of four-level system has been deduced with the theory of rate equation. Based on the expression, the present paper simulated the variations of ionization probability versus laser intensity, laser pulse duration and collision relaxation rate. It was found that in the 1+2+1 ionization mechanism the ionization probability increases with laser intensity, until saturation phenomenon appears in one and two excitation steps. If laser intensity increases further, the ionization probability will oscillate around the saturation value, and the oscillation amplitude will increase with laser intensity. With regard to the influence of laser pulse duration and collision relaxation rate on the ionization probability, the results show that the ionization probability increases from zero to the saturation value 1 with the increase in laser pulse duration, while it decreases linearly with the increase in collision relaxation rate.