Motivated by the reservoir engineering concept of the well productivity index (PI) we study a time dependent functional for general nonlinear Forchheimer equation. The PI of the well characterizes the well capacity with respect to drainage area of the well. Unlike the linear case for which this concept is well developed, there are only a few recent publications dedicated to the PI for nonlinear case. In this paper the PI is comprehensively studied both theoretically and numerically. The impact of the nonlinearity of the flow filtration on the value of the PI is analyzed. Exact formula for the so called "skin factor" in radial case is derived depending on the rate of the flow, the order of nonlinearity and the geometric parameters. Dynamics of the PI for the class of boundary conditions is studied and its convergence to the specific value of steady state PI was justified. Developed framework is applied to obtain nonlinear analog of Peaceman formula for the well-block pressure in unstructured grid. Numerical simulations sustain theoretical results.