Based on 3D Biot’s consolidation theory and nonlinear Duncan-Chang’s model, a 3D FEM (finite element method) program is developed considering the coupling of groundwater seepage and soil skeleton deformation during excavation. The comparison between the analysis result considering the variation of water head difference and that without considering it shows that the porewater pressure distribution of the former is distinctly different from that of the latter and that the foundation pit de- formations of the former are larger than those of the latter, so that the result without considering the variation of water head dif- ference is unreliable. The distribution rules of soil horizontal and vertical displacements around the pit and excess porewater pressure are analyzed in detail in time and space, which is very significant for guiding underground engineering construction and ensuring environment safety around the pit. Based on 3D Biot's consolidation theory and nonlinear Duncan-Chang's model, a 3D FEM (finite element method) program is developed considering the coupling of groundwater seepage and soil skeletal deformation during excavation. The comparison between the analysis result considering the variation of water head difference and that without considering it shows that the porewater pressure distribution of the former is distinctly different from from that of the latter and that the foundation pit de- formations of the former are larger than those of the latter, so that the result without considering the variation of The distribution rules of soil horizontal and vertical displacements around the pit and excess porewater pressure are analyzed in detail in time and space, which is very significant for guiding underground engineering construction and ensuring environment safety around the pit.