Asymptotic behavior of the unique solution for a fractional Kirchhoff problem with singularity

In this paper, we consider the following fractional Kirchhoff problem with singularity $ \left \{\begin{array}{lcl} \Big(1+ b\int_{\mathbb{R}^3}\int_{\mathbb{R}^3} \frac{|u(x)-u(y)|^2}{|x-y|^{3+2s}}\mathrm{d}x \mathrm{d}y \Big)(-\Delta)^s u+V(x)u = f(x)u^{-\gamma}, &&\quad x\in\mathbb{R}^3,\\ u>0,&&\quad x\in\mathbb{R}^3, \end{array}\right. $ where $ (-\Delta)^s $ is the fractional Laplacian with $ 0