Two-photoactive-center model applied to photorefractive self-focusing in biased LiNbO3