Leading logarithmic QCD corrections to the Bs-->gammagamma decays in the two-Higgs doublet model