This study proposes a generalized Darcy's law with considering phase lags in both the water flux and drawdown gradient to develop a lagging flow model for describing drawdown induced by constant-rate pumping (CRP) in a leaky confined aquifer. The present model has a mathematical formulation similar to the dual-porosity model. The Laplace-domain solution of the model with the effect of wellbore storage is derived by the Laplace transform method. The time-domain solution for the case of neglecting the wellbore storage and well radius is developed by the use of Laplace transform and Weber transform. The results of sensitivity analysis based on the solution indicate that the drawdown is very sensitive to the change in each of the transmissivity and storativity. Also, a study for the lagging effect on the drawdown indicates that its influence is significant associated with the lag times. The present solution is also employed to analyze a data set taken from a CRP test conducted in a fractured aquifer in South Dakota, USA. The results show the prediction of this new solution with considering the phase lags has very good fit to the field data, especially at early pumping time. In addition, the phase lags seem to have a scale effect as indicated in the results. In other words, the lagging behavior is positively correlated with the observed distance in the Madison aquifer.