^{1}

^{2}

^{*}

^{2}

^{1}

On the basis of similar structure of solutions of ordinary differential equation (ODE) boundary value problem, the similar construction method was put forward by solving problems of fluid flow in porous media through the homogeneous reservoir. It is indicate that the pressure distribution of dimensionless reservoir and bottom hole in Laplace space, which take on the radial flow, also shows similar structure, and the internal relationship between the above solutions were illustrated in detail.

Due to the permeability of porous media, reservoir engineers could simulate fluid flow by using media such as ground rock, filters and catalyst beds as well, it is useful to those researches who are interested in the behavior of porous media in different engineering applications. Mechanics of porous media flow plays an important role in many branches of engineering, including geotechnical engineering, material science, biomechanics and petroleum industry.

In 2004, S. C. Li [1,2] proposed a important conjecture that the solution formula of some differential equations under different conditions have the similarity, which is similar to that real numbers can be expressed as continued fraction and geometric graphics have certain similarity, we call it the similar structure of solutions. Over the past six years, for some second-order homogeneous ODE [1-14], and some second-order homogeneous linear partial differential equations (PDE) in Laplace space [15-21] as well as some mathematical models of fluids flow in porous media [22-34], the right smart evolution had been obtained on the study of the similar structure of their solutions.

The previous studies are separate, since the similar structure of solutions merely deduced from one certain model of the differential equation boundary value problem, and that was apparently gone against the analysis of the internal relationship between solutions of different models. This paper is intended to reveal the internal relationships between different mathematical models. First of all, in section two, we provide the theoretical background materials of the similar structure of solutions for solving the modified Bessel equation boundary value problem. Secondly, in section tree and four, we devote to solve the problem of fluid flow in porous media through the homogeneous reservoir. At last, in the fifth section, we reveal the inherent laws between above solutions.

In [

where a, b, , R and Q are real constants, and, its solution can be expressed as the following similar structure form

where is called the similar kernel function and defined as

where

,

denotes the first kind modified Bessel function with order, and denotes the second kind modified Bessel function with order. Furthermore, we have

When in the modified Bessel Equation (1), if we setting (where is real constant, and), then y satisfies the ODE

Similarly as Equations (2) and (3) the boundary value conditions are

Quite similar to solve the ODE boundary value problem, by Equations (9) and (11) have the following similar structure form of solution

where is also called the similar kernel function and defined as

where

In particular, we also have

According to [

where p_{D}, r_{D}, t_{D}, C_{D}, R_{D}, and q_{D} are dimensionless variables respectively represent the pressure of any point in reservoir, the distance of reservoir from any point to the wellbore, time, the wellbore storage coefficient, the outer boundary radius of circular reservoir and the liquid flow output of bottom hole (which may depend on time t), and S is the skin factor (well bottom pollution factor). The outer boundary conditions

,

and

denote that the outer boundary of circular reservoir are infinite, constant pressure and closed respectively.

Taking the Laplace transform of with respect to t_{D}, we obtain the ODE with parameter z (where z is variable in Laplace space)

and with the boundary conditions

If let

then the boundary value problem, Equations (23)-(26) becomes the boundary value problem of Equations (9)- (11). Consequently, according to Equation (12), we know that the similar structure formula of dimensionless reservoir pressure distribution at any point for homogeneous reservoir radial flow model in Laplace space is

where is called the similar kernel function and defined as

Similarly, according to Equation (17), we see that the similar structure formula of dimensionless well bottom hole pressure distribution for homogeneous reservoir radial flow model in Laplace space is

In particular, when the liquid flow output of bottom hole is constant, set, then Equations (28) and (32) still hold.

According to [

where, and are given as follows, respectively

, , (37)

The outer boundary conditions

,

and

denote that the outer boundary of circular reservoir are infinite, constant pressure and closed, respectively.

Taking the Laplace transform of with respect to, we obtain the ODE with parameter z (where z is variable in Laplace space)

and with the boundary value conditions

If let

then the boundary value problem Equations (38)-(40) becomes the boundary value problem Equations (9)-(11). Consequently, according to Equation (12), we know that the similar structure formula of dimensionless reservoir pressure distribution at any point for homogeneous reservoir radial flow model in Laplace space is

where is called the similar kernel function and defined as

Similarly, according to Equation (17), we know that the similar structure formula of dimensionless bottom hole pressure distribution for homogeneous reservoir radial flow model in Laplace space is

In particular, when the liquid flow output of bottom hole is constant, set, then Equations (42) and (46) still hold.

Form above discussion, the conclusion can be reached that 1) In section three and four, we use the property of similar structure of ODE boundary value problem, and in the process of solving problem of radial flow through homogeneous reservoir in Laplace space, we can get the expecting outcome by means of algebra constructive theory in Laplace space, and avoid the procedure of solving Bessel equation and the complicated partial derivative operation, but making simple change of variables (such as Equation (27) or Equation (41)) only.

We find this is a simple, convenient and effective method, and an innovation idea for solving problem of fluid flow in porous media.

2) Different similar kernel functions (see Equations (5)-(7), (13)-(15), (29)-(31) and (43)-(45)) corresponding to different right (outer) boundary value conditions (see Equations (3), (11), (22) and (36)), respectively. Actually, the similar kernel functions are the basic solutions which corresponding to the especial left (inner) boundary value conditions (see, , and) respectively.

3) The similar structure formula of solutions (see Equations (4) or (8), (12) or (17), (28) or (32) and (42) or (46)) depend on the left (inner) boundary conditions (see Equations (2), (10), (24), (25) and (39)), instead of the determine equations (see Equations (1), (9), (23) and (38)) and the right (outer) boundary value conditions (see Equations (3), (11), (22) and (36)). From the unified expression in Equation (28) and (32), it is easy to see how the wellbore storage effects, skin factors, as well as the outer boundaries influence the reservoir pressure and bottom hole pressure. This expression brought great convenience for programming well test analysis software, and simplified the program algorithm.

4) If take the Laplace transform to the measured pressure data, well test analysis (dynamic analysis of pressure) can be conducted directly in Laplace space, this also can reflect the advantage of the similar structure of solutions in this paper. In general, using Stehfest numerical inversion formula, the solutions we got in this paper can be transformed into numerical solution in real space, which fully meat the need of the application of well test analysis.

5) This arouses our interest in studying the similar structure of solutions. On the one hand, it is clearly that we can get direct solutions satisfying more especial left (inner) boundary conditions (such as the real constants and in Equation (2) satisfying restricted conditions) from the similar structure expression. On the other hand, using the similar structure expression and the basic solutions (i.e. the similar kernel functions), we can construct some more complicated practical problems. Therefore, we put forward the construction method for solving the boundary value problem of ODE, and it is called the similar construction method.

This work was supported by the Natural Science Key Projects of the Sichuan Education Bureau of China (No. 12ZA164), the Research Fund of Key Disciplinary of Application Mathematics of Xihua University (No. XZD0910-09-1), and the Graduate Innovation Found of Xihua University (No. YCJJ201110).