The similar structure of solutions in fractal multilayer reservoir including a quadratic gradient term

When the classical nonlinear partial differential equations are used to model the fractal reservoir, based on the assumption of low compressibility fluids, the effects of the quadratic gradient term are ignored, which would be questionable for mixed gas reservoirs and low permeability reservoirs. To consider the influence of the wellbore storage, the nonlinear mathematical flow model of the fractal multilayer reservoir is built in this paper, with three kinds of outer boundaries (infinite boundaries, constant pressure boundaries and closed boundaries). Using the Laplace transform method, the solutions for the dimensionless reservoir pressure and the bottom hole pressure in the Laplace space are obtained. An analysis shows that the solutions involve similar structures even for three different kinds of outer boundaries, and can be unified as a continuous fraction. The unified expression would make it more convenient to analyze the formation parameters, which greatly facilitates the development of the well test analysis software.