Aquí está una manera más directa.
La primera cosa a tener en cuenta es la función de $f(z)$ es continua para todos los $z > - 1$, incluyendo en el punto de $z = 1$.
De
$$f(z) = \int_0^\infty \ln (1 + z \operatorname{sech} x) \, dx,$$
sobre la diferenciación con respecto a $z$ hemos
\begin{align}
f'(z) &= \int_0^\infty \frac{\operatorname{sech} x}{1 + z \operatorname{sech}} \, dx\\
&= \int_0^\infty \frac{dx}{\cosh x + z}\\
&= 2 \int_0^\infty \frac{e^x}{e^{2x} + 2z e^x + 1} \, dx.
\end{align}
En la última línea que hemos utilizado el hecho de que $\cosh x = (e^{x} + e^{-x})/2$.
Una sustitución de $u = e^x, du = e^x \, dx$ conduce a
$$f'(z) = 2 \int_1^\infty \frac{du}{(u + z)^2 + 1 - z^2}. \qquad (*)$$
Caso 1: $-1 < z < 1$
Para este caso, como el término lo $(1 - z^2)$ será positivo podemos escribir ($*$) como:
\begin{align}
f'(z) &= 2 \int_1^\infty \frac{du}{(u + z)^2 + (\sqrt{1 - z^2})^2}\\
&= \frac{2}{\sqrt{1 - z^2}} \left [\tan^{-1} \left (\frac{u + z}{\sqrt{1 - z^2}} \right ) \right ]_1^\infty\\
&= \frac{\pi}{\sqrt{1 - z^2}} - \frac{2}{\sqrt{1 - z^2}} \tan^{-1} \left (\frac{1 + z}{\sqrt{1 - z^2}} \right )\\
&= \frac{\pi}{\sqrt{1 - z^2}} - \frac{2}{\sqrt{1 - z^2}} \left [\frac{\pi}{2} - \tan^{-1} \left (\frac{\sqrt{1 - z^2}}{1 + z} \right ) \right ]\\
&= \frac{2}{\sqrt{1 - z^2}} \tan^{-1} \left (\frac{\sqrt{1 - z^2}}{1 + z} \right ), \qquad (**)
\end{align}
donde hemos hecho uso de los resultados $\tan^{-1} (x) = \pi/2 - \tan^{-1} (1/x)$ para $x > 0$.
Tomando nota de que
$$\tan \left (\frac{1}{2} \cos^{-1} z \right ) = \frac{\sqrt{1 - z^2}}{1 + z},$$
el resultado en ($**$) puede escribirse como
$$f'(z) = \frac{1}{\sqrt{1 - z^2}} \cos^{-1} z.$$
La integración con el respeto a $z$ es fácil ver que
$$f(z) = \int \frac{\cos^{-1} z}{\sqrt{1 - z^2}} \, dz = -\frac{1}{2} (\cos^{-1} z )^2 + C.$$
Para encontrar la constante de $C$, la elección de $z= 0$ vemos que $f(0) = 0$. Por lo tanto $C = \pi^2/8$, llevando a
$$f (z) = \frac{\pi^2}{8} - \frac{1}{2} (\cos^{-1} z)^2, \, -1 < z < 1. \qquad (\dagger)$$
Caso 2: $z > 1$
Para este caso, como el término lo $(1 - z^2)$ va a ser negativo, podemos escribir ($*$) como:
\begin{align}
f'(z) &= 2 \int_1^\infty \frac{du}{(u + z)^2 - (\sqrt{z^2 - 1})^2}\\
&= -\frac{2}{\sqrt{z^2 - 1}} \left [\coth^{-1} \left (\frac{u + z}{\sqrt{z^2 - 1}} \right ) \right ]_1^\infty\\
&= \frac{2}{\sqrt{z^2 - 1}} \coth^{-1} \left (\frac{1 + z}{\sqrt{z^2 - 1}} \right )\\
&= \frac{2}{\sqrt{z^2 - 1}} \tanh^{-1} \left (\frac{\sqrt{z^2 - 1}}{z+ 1} \right ), \qquad (***)
\end{align}
donde hemos hecho uso de los resultados $\coth^{-1} (x) = \tanh^{-1} (1/x), x \neq 0$.
Tomando nota de que
$$\tanh \left (\frac{1}{2} \cosh^{-1} z \right ) = \frac{\sqrt{z^2 - 1}}{z + 1},$$
el resultado en ($***$) puede ser escrita como
$$f'(z) = \frac{1}{\sqrt{z^2 - 1}} \cosh^{-1} z.$$
La integración con el respeto a $z$ es fácil ver que
$$f(z) = \int \frac{\cosh^{-1} z}{\sqrt{z^2 - 1}} \, dz = \frac{1}{2} (\cosh^{-1} z)^2 + C.$$
Para encontrar lo desconocido constante $C$, como la función de $f(z)$ es continua para todos los $z > - 1$, $f$ ser continua en $z = 1$ ($\dagger$) requerimos $f(1) = \pi^2/8$. Por lo tanto $C = \pi^2/8$, llevando a
$$f(z) = \frac{\pi^2}{8} + \frac{1}{2} (\cosh^{-1} z)^2, \qquad z > 1.$$