I'll take a statistical field theory approach with which I am more familiar with. But I think it's the same as yours. I'll omit any source term $J$ for simplicity.
The partition function is given by
$$ Z = \mathrm{Tr} \left [ e^{-\beta \hat{H}} \right ] = \sum_\lambda \langle\lambda| e^{-\beta \hat{H}}|\lambda\rangle \\ = \int\mathcal{D}\phi\,\mathrm{e}^{-S_\mathrm{E}[\phi]},$$
where the "Euclidean action" $S_\mathrm{E}$ is related to the "Euclidean Lagrangian density" $\mathcal{L}_{\mathrm{E}}[\phi(x)]$ by $S_\mathrm{E} = \int_0^\beta \mathrm{d}\tau \int\mathrm{d}^3x \, \mathcal{L}_{\mathrm{E}}[\phi(x)]$.
Regarding the limits of integration, especially the imaginary time $\tau$:
If you were to employ periodic boundary conditions, you can set $\phi(\tau=0) = \phi(\beta)$ for bosons, but you would have to impose antiperiodic boundary conditions for fermions $\psi(\tau=0, \mathbf{x}) = -\psi(\tau=\beta, \mathbf{x}) $.
The result of a path integral is a number... and numbers "normally" commute. So, unknowingly and unintentionnaly, the above path integral is only valid for bosons.
In order to describe fermions with a path integral, we need numbers that anticommute. That is Grassmann numbers:
$\eta \xi = -\xi \eta$,
$\eta^2 = 0.$
Then, as shown for example here, you get a similar expression for the partition function:
$$ Z = \int\mathcal{D}\psi\,\mathcal{D}\bar{\psi}\,\mathrm{e}^{-S_\mathrm{E}[\psi, \bar{\psi}]},$$
where you have to use the antiperiodic boundary conditions outlined above.
I don't know how to "show" the existence of the Fermi surface... but I'd say it should be automatic? If two numbers (particles) do not commute, then they cannot live in the same state and they need to "pile up". The level of this "pile" is the Fermi surface. I would assume this translates to the field perspective as well but I wouldn't know how to show it...