In the research field of ground water, hydraulic gradient is studied for decades. In the consolidation field, hydraulic gradient is yet to be investigated as an important hydraulic variable. So, the variation of hydraulic gradient in nonlinear finite strain consolidation was focused on in this work. Based on lab tests, the nonlinear compressibility and nonlinear permeability of Ningbo soft clay were obtained. Then, a strongly nonlinear governing equation was derived and it was solved with the finite element method.Afterwards, the numerical analysis was performed and it was verified with the existing experiment for Hong Kong marine clay. It can be found that the variation of hydraulic gradient is closely related to the magnitude of external load and the depth in soils. It is interesting that the absolute value of hydraulic gradient(AVHG) increases rapidly first and then decreases gradually after reaching the maximum at different depths of soils. Furthermore, the changing curves of AVHG can be roughly divided into five phases. This five-phase model can be employed to study the migration of pore water during consolidation.