# Similar Construction Method of Solution for Solving the Mathematical Model of Fractal Reservoir with Spherical Flow

Authors: Cui-Cui Sheng; Jin-Zhou Zhao; Yong-Ming Li; Shun-Chu Li; Hu Jia
Journal: Journal of Applied Mathematics
 (2013)

---

## Abstract

On the basis of similar structure of solution for a second-order linear differential equation's boundary value problem, and our analysis of the relationship between this similar structure and its kernel function, the differential equation, and the boundary conditions, we propose a new simple solution—similar constructive method of solution (SCMS)—and sum up its detailed steps. We set up a porous media model under three kinds of outer boundary conditions (infinite, constant pressure, and closed), in which the influences of fractal dimension, spherical flow, skin effect, and storage are taken into consideration. And then we use SCMS to solve it. The SCMS is a straightforward method for differential equation's boundary value problems with complex boundary conditions, especially for solving the reservoir models in petroleum engineering.


---

## Body

## 1. Introduction

In the early 21st century, Jia and Li [1] put forward an important hypothesis that the solutions of differential equations under different boundary conditions have similarities—similar structure of solutions, just as a real number can be expressed in continued fraction and geometric graphic have certain similarity. After that, many studies have been carried out on the similar structure of solutions for some boundary value problems of second-order linear homogeneous differential equations [2–8].In the past two decades, many researchers have studied the fractal nature of reservoir rocks and other porous media. Fractal theory can be applied to effectively describe the complexities of a real reservoir, especially of the naturally fractured reservoir or heterogeneous media. Chang and Yortsos [9] were the first to apply the fractal theory to reservoir model, and then set up the mathematical model of fractal reservoir. Beier [10] afterwards extended the fractal model to the study of hydraulically fractured well, and discovered a power-law behavior during the linear and radial flow period. Camacho-Velázquez et al. [11] successfully applied the fractal model to explore the naturally fractured reservoir and obtained a series of decline curves. Using the same model, Cossio [12] proposed a semianalytic solution for finite-conductive flow in vertical fractures, and from the expression of the solution, people can easily perceive the relationship between the fractal parameters {df,θ} and the pressure waves or streamlines. Li et al. [13–16] built mathematical models for fluid flow in porous media, such as homogeneous reservoir mathematical model, double porosity reservoir mathematical model, multilayer reservoir mathematical model, and composite reservoir mathematical model. They found that under three kinds of outer boundary conditions (infinite, constant pressure, and closed), their pressure can also be expressed in similar structures in the Laplace space.When a reservoir is mined with a small opening in the upper reservoir, or the thickness of reservoir stratum is larger than the radius of well bore, there still exists the vertical flow in the flow region, so it cannot be considered as cylindrical flow. Ershaghi et al. [17] analyzed pressure transient data in reservoirs with natural fractures and spherical flow. He found out that the cylindrical flow could not be applied to these reservoirs, or there would be large errors. Proett and Chin [18] provided a new exact spherical flow solution with storage and skin taken into consideration in early-time interpretation with applications to wire line formation and early-evaluation drill stem testing. In recent years, we [19–23] have built some mathematical models for reservoir with spherical flow, and through a complicated process of solution, we obtained their expressions of pressure in the Laplace space.Cheng et al. [24, 25] put forward the idea of SCMS, but they did not specify its detailed procedures. Based on their research, the SCMS is extended in our paper in detail. In this paper, the first two sections focus on solving the boundary value problem of a second-order linear differential equation. Its solution can be expressed in a continued fraction. On the basis of the relationship among similar structure of solution, its kernel function, differential equation, and boundary conditions, we sum up the detailed construction steps of SCMS. In this paper, we also built a new fractal mathematical model for the reservoirs with spherical flow. To apply SCMS to solve this reservoir model, in Section 3 we first study a boundary value problem of the modified Bessel equation, which frequently occurred in reservoir engineering. So we could conveniently construct formation pressure and wellbore pressure. Through the expressions of formation pressure and wellbore pressure and their kernel functions, we could analyze the influences of the variable flow rate, the wellbore storage, the skin effect, the fractal dimension, and the fractal index on them. In a word, SCMS is a straightforward method for boundary value problems of differential equations with complex boundary conditions, especially for solving the reservoir models in petroleum engineering.


## 2. Similar Structure of Solution for the Boundary Value Problem of a Second-Order Linear Differential Equation

The boundary value problem of a second-order linear differential equation:(1)y′′+p(x)y′+q(x)y=0,(2)[ay+(1+ab)y′]x=α=c,(3)[ey+fy′]x=β=0,
where p(x),p′(x), and q(x) are continuous function within [α,β]. a,b,c,e,f are real numbers. When p(x)>0,q(x)≥0 and q(x) is not constantly zero, and a(1+ab)≤0, ef≥0, a2+(1+ab)2≠0, e2+f2≠0, 0<α<β, this boundary value problem has a unique solution [26]; this solution can be expressed in the following continued fraction:
(4)y(x)=c·1a+1/(b+ϕ(α))·1b+ϕ(α)·ϕ(x),
where ϕ(x) is a similar kernel function and defined as
(5)ϕ(x)=eφ0,0(x,β)+fφ0,1(x,β)eφ1,0(α,β)+fφ1,1(α,β).We assume thaty1(x) and y2(x) are the linearly independent solutions of the differential equation (1); then we define a binary function and calculate its partial derivative for x1, x2 respectively, listed as follows:
(6)φ0,0(x1,x2)=y1(x1)y2(x2)-y1(x2)y2(x1),(7)φ0,0(x,β)=y1(x)y2(β)-y1(β)y2(x),(8)φ0,1(α,β)=∂∂x2φ0,0(x1,x2)|(α,β)=y1(α)y2′(β)-y1′(β)y2(α),(9)φ1,0(α,β)=∂∂x1φ0,0(x1,x2)|(α,β)=y1′(α)y2(β)-y1(β)y2′(α),(10)φ1,1(α,β)=∂∂x1φ0,1(x1,x2)|(α,β)=∂∂x2φ1,0(x1,x2)|(α,β)=y1′(α)y2′(β)-y1′(β)y2′(α).


## 3. Similar Construction Method of Solution (SCMS)

On the basis of the relationship among the similar structure of solution (4), its similar kernel function (5), and boundary conditions (2)-(3), the SCMS steps are presented as below.Step 1.
Solve the linearly independent solutionsy1(x) and y2(x) of the differential equation (1).Step 2.
Construct a binary functionφ0,0(x1,x2) with y1(x) and y2(x) (e.g., (6)) and calculate its partial derivative for x1, x2, respectively (e.g., (7)–(10)).Step 3.
Construct a similar kernel functionϕ(x), (e.g., (5)) with the coefficients of the right boundary condition (3) and the binary function φ0,0(x1,x2) and its partial derivative for x1, x2.Step 4.
Construct similar structure of solution (4) with the similar kernel function ϕ(x) and the coefficients of the left boundary condition (2).


## 4. SCMS for Solving the Boundary Value Problem of the Modified Bessel Equation

To apply SCMS to solve the reservoir model, we first introduce the boundary value problem of the modified Bessel equation, which frequently occurred in reservoir engineering:(11)x2y′′+Axy′+(E-Fxd′)y=0,(12)[a′y+(1+a′b′)y′]x=α′=c′,(13)[e′y+f′y′]x=β′=0,
where a′,b′,c′,d′,e′,f′,α′,β′,A,E, and F are real numbers. And they meet the conditions
(14)(1-A)2≥4E,F>0,a′(1+a′b′)≤0,e′f′≥0,d′≠0,a′2+(1+a′b′)2≠0,(e′)2+(f′)2≠0,0<α′<β′.According to SCMS in Section3, we can construct the similar structure and its similar kernel functions of the solutions of the boundary value problem (11)–(13).Step 1.
The linearly independent solutions of the modified Bessel equation (11) are as follows (with details in Appendix A):
(15)x(1-A)/2Iν(2Fd′xd′/2),x(1-A)/2Kν(2Fd′xd′/2),
where Ii(·) denotes the first type of the modified Bessel function with order i, Ki(·) denotes the second type of the modified Bessel function with order i, and ν=(1-A)2-4E/d′.Step 2.
Construct a binary function:(16)M0,0(x,τ)=x(1-A)/2Kν(2Fd′xd′/2)×τ(1-A)/2Iν(2Fd′τd′/2)-τ(1-A)/2Kν(2Fd′τd′/2)×x(1-A)/2Iν(2Fd′xd′/2)=(xτ)(1-A)/2ψν,ν(xd′/2,τd′/2,2Fd′),
where
(17)ψm,n(xε,yε,ξ)=Km(ξxε)In(ξyε)+(-1)m-n+1Im(ξxε)Kn(ξyε),
where m and n are real numbers.
Then, calculate the binary function’s partial derivatives forx,τ:
(18)M0,1(x,τ)=∂∂τM0,0(x,τ)=x(1-A)/2τ-(1+A)/2×[1-A+(1-A)2-4E2×ψν,ν(xd′/2,τd′/2,2Fd′)+Fτd′/2ψν,ν+1(xd′/2,τd′/2,2Fd′)1-A+(1-A)2-4E2],M1,0(x,τ)=∂∂xM0,0(x,τ)=x-(1+A)/2τ(1-A)/2×[1-A+(1-A)2-4E2×ψν,ν(xd′/2,τd′/2,2Fd′)-Fxd′/2ψν+1,ν(xd′/2,τd′/2,2Fd′)1-A+(1-A)2-4E2],M1,1(x,τ)=∂∂xM0,1(x,τ)=∂∂τM1,0(x,τ)=(xτ)-(1+A)/2×{[1-A+(1-A)2-4E]22×ψν,ν(xd′/2,τd′/2,2Fd′)+F[1-A+(1-A)2-4E]2×τd′/2ψν,ν+1(xd′/2,τd′/2,2Fd′)-F[1-A+(1-A)2-4E]2×xd′/2ψν+1,ν(xd′/2,τd′/2,2Fd′)-F(xτ)d′/2ψν+1,ν+1(xd′/2,τd′/2,2Fd′)[1-A+(1-A)2-4E]22}.Step 3.
Construct similar kernel functions with the coefficients of the right boundary condition (13), the binary function (16), and its partial derivatives (18):
(19)ϕ(x)=e′M0,0(x,β′)+f′M0,1(x,β′)e′M1,0(α′,β′)+f′M1,1(α′,β′).Step 4.
Construct the similar structure of solution with the similar kernel function (19) and coefficients of the left boundary condition (12):
(20)y(x)=c′·1a′+1/(b′+ϕ(α′))·1b′+ϕ(α′)·ϕ(x).Consequently, (20) is the solution of the boundary value problem of the modified Bessel equation (11)–(13), which is similar to the real numbers and can be expressed as continued fraction.


## 5. SCMS for Solving the Mathematical Model of Fractal Reservoir with Spherical Flow

In AppendixB, the dimensionless mathematical model of fractal reservoir with spherical flow is presented. It takes dimensionless pressure as variable, and accords with Darcy flow. The details are as follows.Fundamental differential equation:(21)∂2pD∂rD2+ζrD∂pD∂rD=rDθe2S∂pD∂tD,rD>1,tD>0,where:(22)ζ=df-d-θ+2.Initial condition:(23)pD(rD,0)=0.Inner boundary condition:(24)dpwDdrD=-qD(tD)+CDdpwDdtD.Three kinds of outer boundary conditions:(25)pD(∞,tD)=0,pD(RD,tD)=0,∂pD∂rD|rD=RD=0.They denote that the outer boundaries of circular reservoir are infinite pressure, constant pressure, and closed, respectively.All the notations are explained in the nomenclature section.Taking the Laplace transform ofpD(rD,tD) with respect to tD and using the initial condition (23), we get the boundary value problem of the ordinary differential equation with parameter z (where z is Laplace space variable); that is,(26)d2p-DdrD2+ζrDdp-DdrD=rDθe2Szp-D,[-CDzp-wD(rD,z)+ddrDp-D(rD,z)]|rD=1=-q-D(z),p-D(∞,z)=0orp-D(RD,z)=0ordp-DdrD|rD=RD=0.With reference to SCMS in Section3 and its comparison with the boundary value problem of the modified Bessel equation, we can infer
(27)y=p-D,x=rD,A=ζ,E=0,F=ze2S,a′=-CDz,b′=0,c′=-q-D(z),α′=1,ν=ζ-12+θ
and setting e′=1, f′=0, β′→∞ or e′=1, f′=0, β′=RD or e′=0, f′=1, β′=RD represent the three kinds of outer boundary conditions, respectively. So we can construct the similar structure of solution of the dimensionless mathematical model of fractal reservoir with spherical flow (26):
(28)p-(rD,z)=-q-D(z)·1-CDz+1/ϕ(1,z)·1ϕ(1,z)·ϕ(rD,z),
where ϕ(rD,z) is also a similar kernel function. Its expression is the following.Forp-D(∞,z)=0,
(29)ϕ(rD,z)=-Kν((2/(2+θ))z/e2SrD(2+θ)/2)rDθ/2z/e2S[Kν+1((2/(2+θ))z/e2SrD(2+θ)/2)].In (29), we use the properties of limx→∞Ii(x)=∞ and limx→∞Ki(x)=0 to obtain the result above.Forp-D(RD,z)=0, (30)ϕ(rD,z)=-ψν,ν(rD(2+θ)/2,RD(2+θ)/2,(2/(2+θ))z/e2S)zrDθ/e2Sψν+1,ν(rD(2+θ)/2,RD(2+θ)/2,(2/(2+θ))z/e2S).Fordp-D/drD|rD=RD=0, (31)ϕ(rD,z)=-ψν,ν+1(rD(2+θ)/2,RD(2+θ)/2,(2/(2+θ))z/e2S)zrDθ/e2Sψν+1,ν+1(rD(2+θ)/2,RD(2+θ)/2,(2/(2+θ))z/e2S).Thus, (28) is the similar structure of formation pressure of the fractal reservoir with spherical flow in the Laplace space.Then, settingrD=1, we can obtain the similar structure of wellbore pressure of the fractal reservoir with spherical flow in the Laplace space:
(32)p-w(z)=-q-D(z)·1-CDz+1/ϕ(1,z).Substituting (29), (30), and (31) into (28) and (32), respectively, we can get the expressions of formation pressures and the wellbore pressures of thefractal reservoir with spherical flow under three kinds of boundary conditions in the Laplace space. Then, they can be inverted numerically from the Laplace space by using the Stehfest algorithm [27].Note that the solution described above is a general spherical flow model. Its special case are constant production rate(q(t)=q0), no skin effect (S=0), no wellbore storage (C=0), and homogeneous reservoir (df=0,θ=0).


## 6. Conclusions

Based on the reasoning mentioned above, the following conclusions can be drawn.(1)
When dealing with the boundary value problem of a general second-order linear differential equation, we find that the kernel function is related to right (outer) boundary condition and similar structure of solution is related to left (inner) boundary condition. So the SCMS is proposed to solve such complicated problems in reality.(2)
Many reservoir models can be adapted to the boundary value problem of the modified Bessel equation. So SCMS is a convenient, effective, and creative method to solve the problem of fluid flow in porous media since it can avoid complex procedures. SCMS optimizes the programming process of well test analysis software and also simplifies the arithmetic process.(3)
The mathematical model of fractal reservoir with spherical flow is a new model. It is particularly useful when the thick reservoir is mined in an imperfect model. On the basis of the continued fractions (see (28) and (32)), we can directly analyze the wellbore storage’s influence on the formation pressure and wellbore pressure. Based on (29)–(31), it is convenient to see how the fractal dimension, fractal exponent, skin effect, and the three outer boundary conditions influence the kernel functions. Then, substituting (29), (30), and (31) into (28) and (32), respectively, we can see the influence of these parameters on the formation pressure and wellbore pressure.(4)
After the pressure data is processed with the Laplace transformation, well test analysis can be conducted directly in the Laplace space, which can reflect the advantage of the similar structure of solutions. Finally, we can use the Stehfest algorithm to obtain the solution in real space, which meets the demand of well test analysis.



---

*Source: 219218-2013-05-27.xml*