Helium ionization plays an important role in the energy balance of the upper chromosphere and transition region. Helium spectral lines are also often used as diagnostics of these regions. We carry out one-dimensional radiation-hydrodynamics simulations of the solar atmosphere and find that the helium ionization is set mostly by photoionization and direct collisional ionization, counteracted by radiative recombination cascades. By introducing an additional recombination rate mimicking the recombination cascades, we construct a simplified three-level helium model atom consisting of only the ground states. This model atom is suitable for modeling nonequilibrium helium ionization in three-dimensional numerical models. We perform a brief investigation of the formation of the He I 10830 and He II 304 spectral lines. Both lines show nonequilibrium features that are not recovered with statistical equilibrium models, and caution should therefore be exercised when such models are used as a basis for interpretating observations. © 2014 The American Astronomical Society