Let us consider the Schwarzschild metric, assuming as signature convention the negative sign on the time-time component of the metric.
$ds^2 = -g dt^2 + g^{-1} dr^2 + r^2 (d\theta^2 + \sin^2 \theta d\phi^2)$
where:
$G = 1$
$c = 1$
$g = (1 - 2M/r)$
We have three reference frames:
Schwarzschild with coordinates $(t, r, \theta, \phi)$
Stationary observer with coordinates $(\tau, r_{stat}, \theta, \phi)$$(\tau, r_{stat}, \theta_{stat}, \phi_{stat})$
Free falling observer with coordinates $(t', r', \theta, \phi)$
The coordinates $\theta$ and $\phi$ are in common.$(t', r', \theta', \phi')$
Stationary observer at constant $r$, $\theta$ and $\phi$
The time basis is built on the four-velocity $U_{stat} = \partial_\tau = dt / d\tau \partial_t$, with $dt / d\tau = g^{-1/2}$, giving $e_0 = \partial_\tau = g^{-1/2} \partial_t$. It is normalized with squared norm = $-1$. The radial basis is built as $\partial_r$ then normalized with squared norm = $+1$, giving $e_1 = \partial_{r_{stat}} = g^{1/2} \partial_r$. That is
$e_0 = (1 - 2M/r)^{-1/2} \partial_t$
$e_1 = (1 - 2M/r)^{1/2} \partial_r$
The $e_0$ and $e_1$ are orthonormal.
Free falling observer from rest at infinity
It is a radial path, that is at constant $\theta$ and $\phi$. You can relate the free falling frame to the stationary frame with the Lorentz transformation
$\tau = \gamma t' + \gamma v r'$
$r_{stat} = \gamma v t' + \gamma r' $
$v = -(2M/r)^{1/2}$
The velocity $v$ is got by comparing the energy of the free falling as measured by the stationary observer calculated both as $E = -p_\mu U_{stat}^\mu$ and as $E = \gamma m$ with $\gamma = (1 - v^2)^{-1/2}$.
The relation between the partial derivatives is
$\partial_{t'} = \partial \tau / \partial t' \partial_\tau + \partial r_{stat} / \partial t' \partial_{r_{stat}}$
$\partial_{r'} = \partial \tau / \partial r' \partial_\tau + \partial r_{stat} / \partial r' \partial_{r_{stat}}$
where:
$\partial \tau / \partial t' = \gamma$
$\partial r_{stat} / \partial t' = \gamma v$
$\partial \tau / \partial r' = \gamma v$
$\partial r_{stat} / \partial r' = \gamma$
As you have
$e'_0 = \partial_{t'}$
$e'_1 = \partial_{r'}$
you can write
$e'_0 = \gamma e_0 + \gamma v e_1$
$e'_1 = \gamma v e_0 + \gamma e_1$
expressing against $\partial_t$ and $\partial_r$
$e'_0 = \gamma g^{-1/2} \partial_t + \gamma v g^{1/2} \partial_r$
$e'_1 = \gamma v g^{-1/2} \partial_t + \gamma g^{1/2} \partial_r$
and remberingremembering that
$\gamma = g^{-1/2}$
$g = (1 - 2M/r)$
$v = -(2M/r)^{1/2}$
we have
$e'_0 = g^{-1} \partial_t + v \partial_r = (1 - 2M/r)^{-1} \partial_t - (2M/r)^{1/2} \partial_r$
$e'_1 = v g^{-1} \partial_t + \partial_r = -(2M/r)^{1/2} (1 - 2M/r)^{-1} \partial_t + \partial_r$
The $e'_0$ and $e'_1$ are orthonormal as well.
Note:
To complete the orthonormal basis, we have also
$e_2 = 1 / r \partial_\theta$$e_2 = e'_2 = 1 / r \partial_\theta$
$e_3 = 1 / (r \sin \theta) \partial_\phi$$e_3 = e'_3 = 1 / (r \sin \theta) \partial_\phi$