Conventional models for fluid flow in well tests have not been consistent with material balance. According to the slightly compressible fluid assumption, the quadratic gradient term in the nonlinear partial differential equation has been usually neglected. This approach is questionable for live oil and low permeability reservoirs. We have already known that linearization by neglecting quadratic gradient terms may lead to errors for large values of well-test time. In this paper, a method that is consistent with material balance was proposed on the spherical flow system. All terms in the nonlinear partial eqiation were retained. Exact solution for the resulting nonlinear partial differential equation in an infinite reservoir was obtained by using the Laplace transform considering wellbore storage. Analytical solution for nonlinear partial differential equation are resulted by using orthogonal transforms under both closed and constant-pressure outer boundary conditions. The law of pressure changes for a fluid compressibility α and a storage coefficient CD were discussed.