The nonlinear dual-porosity flow model, specifically considering the quadratic pressure gradient term, wellbore storage coefficient, well skin factor, and interporosity flow of matrix to natural fractures, was established for well production in a naturally fractured formation and then solved using a semianalytical method, including Laplace transform and a transformation of the pressure function. Analytical solution of the model in Laplace space was converted to numerical solution in real space using Stehfest numerical inversion. Nonlinear flow process for well production in a naturally fractured formation with different external boundaries was simulated and analyzed using standard pressure curves. Influence of the quadratic pressure gradient coefficient on pressure curves was studied qualitatively and quantitatively in conditions of a group of fixed model parameters. The research results show that the semianalytical modelling method is applicable in simulating the nonlinear dual-porosity flow behavior.