From Gerald Edgar's 1991/1992 explanation of Pi in the Mandelbrot set, we learn that iterating a function of the type $x \mapsto x^2+x+\epsilon\;$ will take approximately $\frac{\pi}{\sqrt{\epsilon}}$ iterations to escape, where we start from the critical point at x=-0.5. From Gerald Edgar's Pi and the Mandelbrot set, "So our equation now reads $y'(n) = y^2 + \epsilon$. This has the solution $a\tan(an+c)$ where $a = \sqrt{\epsilon}$"
It turns out we can we put the equation iterating the Op's tetration expression into a similar form, for $$b>\exp\left(\frac{1}{e}\right);\;\;\; x \mapsto b^x$$
First we observe that if $\epsilon=\ln(\ln(b))+1$, and $y=x\cdot\ln(b)+(-1+\epsilon)$, then iterating
$$y \mapsto \exp(y)-1+\epsilon\;\;\;\text{is exactly congruent to}\;\;\; x \mapsto b^x$$
Also notice that $\epsilon$ approaches zero as b approaches $\exp(1/e)$, and $\exp(y)-1=y+\frac{y^2}{2}+\frac{y^3}{6}...$ so this is also close to the desired form except it has $\frac{y^2}{2}$ instead of $y^2$
So instead, we need $$z=\frac{y}{2}= \frac{x\cdot\ln(b)+(-1+\epsilon)}{2}$$ and then we iterate
$$z \mapsto \frac{\exp(2z)-1+\epsilon}{2}\;\;\;\text{is exactly congruent to}\;\;\; x \mapsto b^x$$
And $\frac{\exp(2z)-1}{2}=z+z^2+\frac{2z^3}{3}+\frac{z^4}{3}...$ has exactly the desired form to approximate $\frac{\pi}{\sqrt{\epsilon}}$ iterations by the tangent approximation; with the Op using x as the base: find n such that $x \uparrow \uparrow n = 10^{100}$, where $10^{100}$ is a good approximation for escapes to infinity. With a little algebra, we need an equation for $\epsilon$ in terms of x from the Op's equation. Here; I suggest using the shorthand $\eta=\exp(1/e)$, and then if $x>\eta$, but only a little bigger than $\eta$.
$$\epsilon = \ln\left(\ln\left(x + \eta-\eta \right)\right)+1$$
$$\epsilon = \frac{e\cdot \left(x-\eta\right)}{\eta} + \mathcal{O}
\left(x-\eta\right)^2$$
We are iterating $z\mapsto f(z)+\epsilon/2\;\;$ so we expect the approximation for the total escape time to be
$$n \approx \frac{\pi}{\sqrt{x-\eta}} \cdot \sqrt{\frac{2\eta}{e}}$$
$$n \approx \frac{\pi}{\sqrt{x-\eta}}\sqrt{\frac{2\eta}{e}};\;\;\;x\uparrow \uparrow n >\approx 10^{100}$$
For $x=\eta+10^{-10}\;\;\;n\approx 323893.054$ so I would use this approximation. I give both the equation for n in terms of x, and for x in terms of n.
$$n \approx -0.94 + \frac{\pi}{\sqrt{x-\eta}}\sqrt{\frac{2\eta}{e}};\;\;\;x\approx \eta+\frac{2\eta\pi^2}{e}\left(\frac{1}{n+0.94}\right)^2;\;\;\;x\uparrow \uparrow n \approx 10^{100}$$
This approximation seems to work remarkably well for values $x=\eta+\epsilon$ where $\epsilon<2.5\cdot 10^{-3}$. This approximation for the number of iterations required is far more accurate than one would naively expect. For example, it is accurate with an error of 0.11 or less iterations when $60 \leq n \leq 10^8$; For example, when used with the correct value of x for n=10^7, (found by iteration), the equation for n predicts 10000000.00141 iterations. The error 10^6,10^7, and 10^8 are almost the same.