The dynamics of generalized Forchheimer equations for slightly compressible fluids are tudied by means of initial boundary value problem for the pressure. We prove that the solutions continuously depend on the boundary data and the Forchheimer polynomials. New bounds for the solutions are established in the Lα-norm for all α ≥ 1, and then are used to improve estimates for their spatial and time derivatives. New Poincaré–Sobolev inequalities and nonlinear Gronwall type estimates for nonlinear differential inequalities are utilized to achieve better asymptotic bounds. The methods developed can be applied to other degenerate parabolic equations. Bibliography: 25 titles.