  We obtain predictions for $B_{(s)}\to e^+e^-\gamma$ and $B_{(s)}\to
\mu^+\mu^-\gamma$ decays. All the contributions containing long-distance QCD
effects are calculated in the framework of relativistic quark model. The
contributions of the light vector-meson resonances related to the virtual
photon emission from valence quarks of the $B$-meson are included. The highest
branching ratios for the radiative leptonic B-decays are ${\cal B}(\bar
B^0_s\to e^+e^-\gamma)=18.8\times10^{-9}$ and ${\cal B}(B^0_s\to
\mu^+\mu^-\gamma)=12.2\times10^{-9}$. We also give the distribution of the
foward-backward asymmetry.
