The interplay between symmetry and topological properties plays a very important role in modern physics. In the past decade, the concept of symmetry-enriched topological (SET) phases was proposed and their classifications have been systematically studied for bosonic systems. Very recently, the concept of SET phases has been generalized into fermionic systems and their corresponding classification schemes are also proposed. Nevertheless, how to realize all these fermionic SET (fSET) phases in lattice models remains to be a difficult open problem. In this paper, we first construct exactly solvable models for nonanomalous nonchiral fSET phases, namely, the symmetry-enriched fermionic string-net models, which are described by commuting-projector Hamiltonians whose ground states are the fixed-point wave functions of each fSET phase. Mathematically, we provide a partial definition to -graded super fusion category (when the total symmetry is a direct product of physical symmetry and fermion-parity symmetry ), which is the input data of a symmetry-enriched fermionic string-net model. Next, we construct exactly solvable models for nonchiral fSET phases with ’t Hooft anomaly, especially the fermionic ’t Hooft anomaly which is different from the well-known bosonic anomaly. In our construction, this fermionic ’t Hooft anomaly is characterized by a violation of fermion-parity conservation in some of the surface moves (a type of renormalization moves for the ground state wave functions of surface fSET phases), and also by a new fermionic obstruction in the surface pentagon equation. We demonstrate this construction in a concrete example that the surface topological order is a gauge theory embedded into a fermion system and the total symmetry . We further conjecture that the fermionic ’t Hooft anomaly is characterized by a possible violation of nonterminating -type (called -type) string in the fusion rules of surface models.