To elucidate the regularities inherent in the kinetics of ultrafast charge recombination ensuing photoinduced charge separation in donor-acceptor dyads in solutions, the simulations of the kinetics have been performed within the stochastic multichannel point-transition model. Increasing the solvent relaxation time scales has been shown to strongly vary the dependence of the charge recombination rate constant on the free energy gap. In slow relaxing solvents the non-equilibrium charge recombination occurring in parallel with solvent relaxation is very effective so that the charge recombination terminates at the non-equilibrium stage. This results in crucial difference between the free energy gap laws for the ultrafast charge recombination and the thermal charge transfer. For the thermal reactions the well-known Marcus bell-shaped dependence of the rate constant on the free energy gap is realized while for the ultrafast charge recombination only descending branch is predicted in the whole area of free energy gap exceeding 0.1 eV. From the available experimental data on population kinetics of the second and first excited states for a series of the Zn-porphyrin-imide dyads in toluene and tetrahydrofuran solutions, effective rate constants of the charge recombination into the first excited state has been calculated. The obtained rate constants being very high is nearly invariable in the area of the charge recombination free energy gap from 0.1 to 0.4 eV that supports the theoretical prediction.