Consider the simplest QFT, namely the free scalar field. The equation of motion (in the Heisenberg picture) is
$$
(\partial_t^2-\nabla^2 +m^2)\phi(t,\mathbf{x})=0
\tag{1}
$$
and the equal-time commutation relations are
\begin{gather}
[\phi(t,\mathbf{x}),\,\dot\phi(t,\mathbf{y})]
=i\delta^3(\mathbf{x}-\mathbf{y})
\\
[\phi(t,\mathbf{x}),\,\phi(t,\mathbf{y})]=0
\hskip2cm
[\dot\phi(t,\mathbf{x}),\,\dot\phi(t,\mathbf{y})]=0.
\tag{2}
\end{gather}
We want to construct a unitary operator $U(\mathbf{x})$ that satisfies
$$
U(\mathbf{x})\phi(t,\mathbf{0})U^\dagger(\mathbf{x})=
\phi(t,\mathbf{x}),
\tag{3}
$$
which is the definition of the translation operator. The momentum operators $P_k$ are the hermitian generators of the translation group (by definition again), so
$$
U(\mathbf{x})=\exp(i\mathbf{x}\cdot \mathbf{P})
\hskip2cm
\mathbf{x}\cdot \mathbf{P}\equiv\sum_k x_k P_k.
\tag{4}
$$
in units where $\hbar=1$. Take the gradient of (3) with respect to $x_k$ to get
$$
[P_k\,\phi(t,\mathbf{x})]=-i\nabla_k\phi(t,\mathbf{x}).
\tag{5}
$$
This is the defining property of the operators $P_k$. What we want, though, is an explicit expression for $P_k$ in terms of the field operators. In general, we can use Noether's theorem to get an expression for $P_k$ in terms of the field operators. Or, instead of going through Noether's theorem, we can write down an ansatz and then prove that it works. (The second approach is easier when we already know the answer, and since I do already know the answer, I'll use the second approach here.) The commutation relations (2) imply that the operator
$$
P_k=\int d^3x\ \dot\phi(t,\mathbf{x})\nabla_k \phi(t,\mathbf{x})
\tag{6}
$$
is hermitian and satisfies the condition (5), so this ansatz works. Equation (6) expresses the momentum operators in terms of the field operators, and then equation (4) gives the translation operators. If the field operators are written in terms of the usual creation/annihilation operators, then (6) becomes
$$
P_k\propto \int d^3p\ p_k a^\dagger(\mathbf{p})a(\mathbf{p}).
\tag{7}
$$
The post
Derivation of total momentum operator QFT
writes out this last step a little more explicitly. Also see equation (2.21) in
which is one of the first sources I found in a quick search for "momentum operator in QFT."
The preceding equations explain how to express $U(\mathbf{x})$ in terms of the field operators. To address the original question about why the two-point correlation function is translation-invariant, we only need equations (3) and (4), together with the assumption that the ground state $|G\rangle$ is translation-invariant: $U(\mathbf{x})|G\rangle=|G\rangle$. This gives
\begin{align}
\langle G|\phi(\mathbf{x})\phi(\mathbf{x'})|G\rangle
&=
\langle G|U(\mathbf{x})\phi(\mathbf{0})U^\dagger(\mathbf{x})
U(\mathbf{x'})\phi(\mathbf{0})U^\dagger(\mathbf{x'})|G\rangle \\
&=
\langle G|\phi(\mathbf{0})
U(\mathbf{x'}-\mathbf{x})\phi(\mathbf{0})|G\rangle,
\end{align}
which shows that the correlation function depends only on
$\mathbf{x'}-\mathbf{x}$.