A model for calculation of kinetic characteristics of electronic excitation decay in the impurity subsystem of doped crystals with taking into account the phenomenon of excitation trapping is suggested. Numerical calculations are carried out for LiNbO3:Yb3+ crystals: the probabilities of elementary acts of resonant and non-resonant nonradiative redistributions of electronic excitation energy in the impurity subsystem are calculated and the dependences of excitation lifetime and luminescence quantum yield on the concentration of pair centers in the matrix of LiNbO3 crystals are determined.