Leading logarithmic QCD corrections to the Bs → γγ decays in the two-Higgs doublet model