The integral doesn't appear to admit a closed form solution. But, with b=1 (which is the lower bound for b),we get:
$$\phi(x;\alpha_0,\sigma)=\frac{1}{x^2 }\sum _{i=0}^\infty \frac{1}{i!}L (\alpha_0-1)^i e^{\frac{1}{2} i (i-1) \sigma ^2} \frac{(i+\log (\frac{L}{x})}{\log (\frac{L}{x})^{1-i} }$$
This result is obtained by expanding $\alpha$ around its lower bound $b$ (which we simplified to $b=1$) and integrating each summand.
We know that as $x \to \infty$ the density $\phi$ collapses to power laws with a tail $\alpha=1$, i.e., where K is a norming constant, $\phi(x)= K x^{-2}$.
$\textbf{Expectation: }$ Even in the absence of a fully explicit density, we can extract the explicit expectation as follows. By integrating first with respect to $x$:
$$ \int_b^\infty \int_L^\infty x \phi(x;\alpha)\,\mathrm{d}x \mathrm{d}\alpha= \int_L^\infty \frac{\alpha L \exp \left(-\frac{\left(2 \log (\alpha -b)-2 \log (\text{$\alpha $0}-b)+\sigma ^2\right)^2}{8 \sigma ^2}\right)}{\sqrt{2 \pi } (\alpha -1) \sigma (\alpha -b)}\mathrm{d}\alpha$$
With, again, b=1:
$$ \mathbb{E}(X)= \frac{L \left(\alpha_0+e^{\sigma ^2}-1\right)}{\alpha_0-1}$$
Note that
\begin{align}
&F(x) = \int\limits_0^xf(\xi)\,\mathrm d\xi = \frac{\alpha ^{\frac\alpha2}}{B\left(\frac12,\frac{\alpha }{2}\right)}\int\limits_0^x \xi^{-\frac12}(\alpha +\xi)^{\frac12 (-\alpha -1)} \,\mathrm d\xi,\\[4pt]
&t=\dfrac{\xi}{\alpha + \xi},\quad \xi=\dfrac{\alpha t}{1-t},\quad \alpha+\xi=\dfrac\alpha{1-t},\quad \mathrm d\xi = \dfrac\alpha{(1-t)^2}\mathrm dt,\\[4pt]
&F(x) = \frac{\alpha ^{\frac\alpha2}}{B\left(\frac12,\frac\alpha2\right)} \int\limits_0^{\frac x{\alpha+x}} \left(\dfrac{\alpha t}{1-t}\right)^{-\frac12}\left(\dfrac\alpha{1-t}\right)^{\frac12 (-\alpha -1)}\dfrac\alpha{(1-t)^2}\mathrm dt\\
&\qquad = \frac1{B\left(\frac12,\frac\alpha2\right)} \int\limits_0^\frac x{\alpha+x} t^{-\frac12}(1-t)^{\frac\alpha2-1}\,\mathrm dt,\\[4pt]
&F(x) = \frac{B_{\frac x{\alpha+x}}\left(\frac12,\frac\alpha2\right)}{B\left(\frac12,\frac\alpha2\right)} = I_{\frac x{\alpha+x}}\left(\frac12,\frac\alpha2\right),\tag1\\[4pt]
\end{align}
where $B_t(a,b)$ is incomplete beta function and $I_t(a,b)$ is regularized beta function.
If $\mathbf{n=1},$ then solution is trivial:
$$f_1(y) = f(y).\tag2$$
If $\mathbf{n=2},$ then
\begin{align}
&f_2(y) = \int\limits_0^yf(\xi)f(y-\xi)\,\mathrm d\xi = f(y)\ast f_1(y),\tag3\\
&f_2(y) = \dfrac{\alpha^\alpha}{B^2\left(\frac12,\frac\alpha2\right)} \int\limits_0^y\dfrac{\left((\alpha+\xi)(\alpha+y-\xi)\right)^{-\frac{\alpha+1}2}}{\sqrt{\xi(y-\xi)}}\,\mathrm d\xi,\\[4pt]
&\xi=\frac y2(1-\cos t),\quad d\xi=\dfrac y2\sin t\,\mathrm dt,\\[4pt]
&f_2(y) = \dfrac{\alpha^\alpha}{B^2\left(\frac12,\frac\alpha2\right)} \int\limits_0^\pi\dfrac{\left(\left(\alpha+\frac y2\right)^2 - \left(\frac y2\sin t\right)^2\right)^{-\frac{\alpha+1}2}}{\frac y2\sin t}\dfrac y2\sin t\,\mathrm dt\\
&\qquad = \dfrac{\alpha^\alpha}{2B^2\left(\frac12,\frac\alpha2\right)} \int\limits_0^{2\pi}\left(\left(\alpha+\frac y2\right)^2 - \left(\frac y2\sin t\right)^2\right)^{-\frac{\alpha+1}2}\,\mathrm dt,\\
&z=e^{it},\quad dt = -i\dfrac{dz}z,\\
&f_2(y) = -i\dfrac{\alpha^\alpha}{B^2\left(\frac12,\frac\alpha2\right)} \oint\limits_{|z|=1}\left(\left(\alpha+\frac y2\right)^2 +\frac{y^2}4 \left(z-\frac 1z\right)^2\right)^{-\frac{\alpha+1}2}\,\dfrac{\mathrm dz}z,\\
\end{align}
wherein
\begin{align}
&\left(\left(\alpha+\frac y2\right)^2 +\frac{y^2}4 \left(z-\frac 1z\right)^2\right)^{-\frac{\alpha+1}2} = \left(\alpha^2+\alpha y - \frac {y^2}4 +\frac{y^2}4 \left(z^2+\frac 1{z^2}\right)\right)^{-\frac{\alpha+1}2}\\[4pt]
&\qquad = \left(\alpha^2+\alpha y - \frac {y^2}4\right)^{-\frac{\alpha+1}2}\left(1+\dfrac{y^2}{4\alpha^2+4\alpha y - y^2}\left(z^2+\frac 1{z^2}\right)\right)^{-\frac{\alpha+1}2}\\[4pt]
& = \left(\alpha^2+\alpha y - \frac {y^2}4\right)^{-\frac{\alpha+1}2}\Bigg(1-\frac{\alpha+1}2\dfrac{y^2}{4\alpha^2+4\alpha y - y^2}\left(z^2+\frac 1{z^2}\right)\\[4pt]
&+\dfrac1{2!}\frac{\alpha+1}2 \frac{\alpha+3}2\left(\dfrac{y^2}{4\alpha^2+4\alpha y - y^2}\right)^2\left(z^2+\frac 1{z^2}\right)^2\\[4pt]
&-\dfrac1{3!}\frac{\alpha+1}2\frac{\alpha+3}2\frac{\alpha+5}2\left(\dfrac{y^2}{4\alpha^2+4\alpha y - y^2}\right)^3\left(z^2+\frac 1{z^2}\right)^3\\[4pt]
&+\dfrac1{4!}\frac{\alpha+1}2\frac{\alpha+3}2\frac{\alpha+7}2\frac{\alpha+9}2\left(\dfrac{y^2}{4\alpha^2+4\alpha y - y^2}\right)^4\left(z^2+\frac 1{z^2}\right)^4 + \dots\Bigg)\\
&\qquad = c_0(y) + c_2(y)\left(z^2+\dfrac1{z^2}\right) + c_4(y)\left(z^4+\dfrac1{z^4}\right) + \dots,\\
&c_0(y) = \left(\alpha^2+\alpha y - \frac {y^2}4\right)^{-\frac{\alpha+1}2}\Bigg(1
+{\dfrac1{2!}\frac{\Gamma\left(\frac{\alpha+5}2\right)}{\Gamma\left(\frac{\alpha+1}2\right)}\binom21\left(\dfrac{y^2}{4\alpha^2+4\alpha y - y^2}\right)^2}\\[4pt]
&+\dfrac1{4!}\frac{\Gamma\left(\frac{\alpha+9}2\right)}{\Gamma\left(\frac{\alpha+1}2\right)}\binom42\left(\dfrac{y^2}{4\alpha^2+4\alpha y - y^2}\right)^4
+\dfrac1{6!}\frac{\Gamma\left(\frac{\alpha+13}2\right)}{\Gamma\left(\frac{\alpha+1}2\right)}\binom63\left(\dfrac{y^2}{4\alpha^2+4\alpha y - y^2}\right)^6+\dots\Bigg),\\[4pt]
&c_0(y) = \left(\alpha^2+\alpha y - \frac {y^2}4\right)^{-\frac{\alpha+1}2}\frac1{\Gamma\left(\frac{\alpha+1}2\right)}\sum\limits_{n=0}^\infty \dfrac{\Gamma\left(\frac{\alpha+4n+1}2\right)}{(n!)^2}\left(\dfrac{y^2}{4\alpha^2+4\alpha y - y^2}\right)^{2n},\tag4\\[4pt]
&f_2(y) = 2\pi c_0(y)\dfrac{\alpha^\alpha}{B^2\left(\frac12,\frac\alpha2\right)}.\tag5\\[4pt]
\end{align}
If $\mathbf{n>2},$ then solution can be obtained as convolution in the forms of
\begin{align}
&f_n(y) = \int\limits_0^yf(y-\xi)f_{n-1}(\xi)\,\mathrm d\xi = f(y)\ast f_{n-1}(y),\tag6\\
\end{align}
or
\begin{align}
&f(y)=\int\limits_0^y\int\limits_0^{y-\xi_1}\dots\int\limits_0^{y-\xi_1\dots-\xi_{n-2}}f(\xi_1)f(\xi_2)\dots f(\xi_{n-1})f(y-\xi_{n-1})\,\mathrm d\xi_1\,\mathrm d\xi_2\dots\,\mathrm d\xi_{n-1}.\tag7\\
\end{align}
Best Answer
The first integral is
$$\frac1{\sqrt{2 \pi} \sigma x} \int_0^{\infty} d\alpha \, \alpha \, e^{\alpha \log{(L/x)}} e^{-(\alpha-\mu)^2/(2 \sigma^2)} $$
This basically has the form
$$C \int_0^{\infty} d\alpha \, \alpha \, e^{-A (\alpha+B)^2} = C \int_B^{\infty} d\alpha \, (\alpha-B) e^{-A \alpha^2}$$
which is an error function plus an exponential term.
The second integral is
$$\frac{b^{-a/b}}{x \Gamma \left ( \frac{a}{b} \right )}\int_0^{\infty} d\alpha \, \alpha^{a/b} \, e^{- \alpha (1/b + \log{(x/L)})} = \frac{a}{b} \frac{b^{-a/b}}{x \left ( \frac1{b} + \log{\left (\frac{L}{x} \right )} \right )^{a/b+1} }$$
because the integral is a Gamma function.