You must log in to edit PetroWiki. Help with editing
Content of PetroWiki is intended for personal use only and to supplement, not replace, engineering judgment. SPE disclaims any and all liability for your use of such content. More information
Solving unsteady flow problems with Green's and source functions: Difference between revisions
No edit summary |
|||
Line 1: | Line 1: | ||
As discussed in [[Source function solutions of the diffusion equation]], the conventional development of the source function solutions uses the instantaneous point-source solution as the building block with the appropriate integration (superposition) in space and time. In 1973, Gringarten and Ramey<ref name="r1" /> introduced the use of the source and [[ | As discussed in [[Source_function_solutions_of_the_diffusion_equation|Source function solutions of the diffusion equation]], the conventional development of the source function solutions uses the instantaneous point-source solution as the building block with the appropriate integration (superposition) in space and time. In 1973, Gringarten and Ramey<ref name="r1">Gringarten, A.C. and Ramey Jr., H.J. 1973. The Use of Source and Green's Functions in Solving Unsteady-Flow Problems in Reservoirs. SPE J. 13 (5): 285-296. SPE-3818-PA. http://dx.doi.org/10.2118/3818-PA</ref> introduced the use of the source and [[Green’s_function_for_solving_transient_flow_problems|Green’s function]] method to the petroleum engineering literature with a more efficient method of developing the source solutions. Specifically, they suggested the use of infinite-plane sources as the building block with Newman’s product method.<ref name="r2">Newman, A.B. 1936. Heating and Cooling Rectangular and Cylindrical Solids. Ind. Eng. Chem. 28 (5): 545–548. http://dx.doi.org/10.1021/ie50317a010</ref> In this page we discuss the use of Green’s functions and source functions in solving unsteady-flow problems in reservoirs. | ||
==Green's functions and source functions in solving unsteady flow problems== | == Green's functions and source functions in solving unsteady flow problems == | ||
[[File:Vol1 page 0110 eq 002.png]]....................(1) | Green’s function for transient flow in a porous medium is defined as the pressure at ''M'' (''x'', ''y'', ''z'') at time ''t'' because of an instantaneous point source of unit strength generated at point ''M′''(''x′'', ''y′'', ''z′'') at time ''τ'' < ''t'' with the porous medium initially at zero pressure and the boundary of the medium kept at zero pressure or impermeable to flow.<ref name="r1">Gringarten, A.C. and Ramey Jr., H.J. 1973. The Use of Source and Green's Functions in Solving Unsteady-Flow Problems in Reservoirs. SPE J. 13 (5): 285-296. SPE-3818-PA. http://dx.doi.org/10.2118/3818-PA</ref><ref name="r3">Carslaw, H.S. and Jaeger, J.C. 1959. Conduction of Heat in Solids, second edition, 353–386. Oxford, UK: Oxford University Press.</ref> If we let ''G''(''M'', ''M′'', ''t'' − ''τ'') denote the Green’s function, then it should satisfy the diffusion equation; that is, | ||
[[File:Vol1 page 0110 eq 002.png|RTENOTITLE]]....................(1) | |||
Because ''G'' is a function of ''t'' − ''τ'', it should also satisfy the adjoint diffusion equation, | Because ''G'' is a function of ''t'' − ''τ'', it should also satisfy the adjoint diffusion equation, | ||
[[File:Vol1 page 0110 eq 003.png]]....................(2) | [[File:Vol1 page 0110 eq 003.png|RTENOTITLE]]....................(2) | ||
Green’s function also has the following properties: <ref name="r1" /><ref name="r3" /> | Green’s function also has the following properties: <ref name="r1">Gringarten, A.C. and Ramey Jr., H.J. 1973. The Use of Source and Green's Functions in Solving Unsteady-Flow Problems in Reservoirs. SPE J. 13 (5): 285-296. SPE-3818-PA. http://dx.doi.org/10.2118/3818-PA</ref><ref name="r3">Carslaw, H.S. and Jaeger, J.C. 1959. Conduction of Heat in Solids, second edition, 353–386. Oxford, UK: Oxford University Press.</ref> | ||
1. ''G'' is symmetrical in the two points ''M'' and ''M′''; that is, Green’s function is invariant as the source and the observation points are interchanged. | 1. ''G'' is symmetrical in the two points ''M'' and ''M′''; that is, Green’s function is invariant as the source and the observation points are interchanged. 2. As ''t'' → ''τ'', ''G'' vanishes at all points in the porous medium; that is, [[File:Vol1 page 0110 inline 001.png|RTENOTITLE]], except at the source location, ''M'' = ''M′'', where it becomes infinite, so that ''G'' satisfies the delta function property, | ||
2. As ''t'' → ''τ'', ''G'' vanishes at all points in the porous medium; that is, [[File:Vol1 page 0110 inline 001.png]], except at the source location, ''M'' = ''M′'', where it becomes infinite, so that ''G'' satisfies the delta function property, | |||
[[File:Vol1 page 0110 eq 004.png]]....................(3) | [[File:Vol1 page 0110 eq 004.png|RTENOTITLE]]....................(3) | ||
where ''D'' indicates the domain of the porous medium, and ''φ''(''M'') is any continuous function. | where ''D'' indicates the domain of the porous medium, and ''φ''(''M'') is any continuous function. 3. Because ''G'' corresponds to the pressure because of an instantaneous point source of unit strength, it satisfies | ||
3. Because ''G'' corresponds to the pressure because of an instantaneous point source of unit strength, it satisfies | |||
[[File:Vol1 page 0110 eq 005.png]]....................(4) | [[File:Vol1 page 0110 eq 005.png|RTENOTITLE]]....................(4) | ||
4. ''G'' or its normal derivative, ''∂G''/''∂n'', vanishes at the boundary, Γ, of the porous medium. If the porous medium is infinite, then ''G'' vanishes as ''M'' or ''M′''→∞. | 4. ''G'' or its normal derivative, ''∂G''/''∂n'', vanishes at the boundary, Γ, of the porous medium. If the porous medium is infinite, then ''G'' vanishes as ''M'' or ''M′''→∞. | ||
Line 26: | Line 25: | ||
Let ''p''(''M′'' , ''τ'') be the pressure in the porous medium and ''G''(''M'', ''M′'', ''t'' - ''τ'') be the Green’s function. Let ''D'' denote the domain of the porous medium. Then, ''p'' and ''G'' satisfy the following differential equations: | Let ''p''(''M′'' , ''τ'') be the pressure in the porous medium and ''G''(''M'', ''M′'', ''t'' - ''τ'') be the Green’s function. Let ''D'' denote the domain of the porous medium. Then, ''p'' and ''G'' satisfy the following differential equations: | ||
[[File:Vol1 page 0111 eq 001.png]]....................(5) | [[File:Vol1 page 0111 eq 001.png|RTENOTITLE]]....................(5) | ||
and | and | ||
[[File:Vol1 page 0111 eq 002.png]]....................(6) | [[File:Vol1 page 0111 eq 002.png|RTENOTITLE]]....................(6) | ||
Then, we can write | Then, we can write | ||
[[File:Vol1 page 0111 eq 003.png]]....................(7) | [[File:Vol1 page 0111 eq 003.png|RTENOTITLE]]....................(7) | ||
or | or | ||
[[File:Vol1 page 0111 eq 004.png]]....................(8) | [[File:Vol1 page 0111 eq 004.png|RTENOTITLE]]....................(8) | ||
where ''ε'' is a small positive number. Changing the order of integration and applying the Green’s theorem, | where ''ε'' is a small positive number. Changing the order of integration and applying the Green’s theorem, | ||
[[File:Vol1 page 0111 eq 005.png]]....................(9) | [[File:Vol1 page 0111 eq 005.png|RTENOTITLE]]....................(9) | ||
where ''D'' and Γ indicate the volume and boundary surface of the domain, respectively; ''S'' denotes the points on the boundary; and ''∂''/''∂n'' indicates differentiation in the normal direction of the surface Γ; we obtain | where ''D'' and Γ indicate the volume and boundary surface of the domain, respectively; ''S'' denotes the points on the boundary; and ''∂''/''∂n'' indicates differentiation in the normal direction of the surface Γ; we obtain | ||
[[File:Vol1 page 0111 eq 006.png]]....................(10) | [[File:Vol1 page 0111 eq 006.png|RTENOTITLE]]....................(10) | ||
Taking the limit as ''ε''→0 and noting the delta-function property of the Green’s function ('''Eq. 3'''), '''Eq. 10''' yields | Taking the limit as ''ε''→0 and noting the delta-function property of the Green’s function ('''Eq. 3'''), '''Eq. 10''' yields | ||
[[File:Vol1 page 0111 eq 007.png]]....................(11) | [[File:Vol1 page 0111 eq 007.png|RTENOTITLE]]....................(11) | ||
where ''p''<sub>''i''</sub>(''M'') = ''p''(''M'', ''t'' = 0) is the initial pressure at point ''M''. | where ''p''<sub>''i''</sub>(''M'') = ''p''(''M'', ''t'' = 0) is the initial pressure at point ''M''. | ||
Line 56: | Line 55: | ||
In '''Eq. 11''', the boundary of the porous medium consists of two surfaces: the inner boundary that corresponds to the surface of the wellbore, Γ<sub>''w''</sub>, and the outer boundary of the reservoir, Γ<sub>''e''</sub>. '''Eq. 11''' may be written as | In '''Eq. 11''', the boundary of the porous medium consists of two surfaces: the inner boundary that corresponds to the surface of the wellbore, Γ<sub>''w''</sub>, and the outer boundary of the reservoir, Γ<sub>''e''</sub>. '''Eq. 11''' may be written as | ||
[[File:Vol1 page 0112 eq 001.png]]....................(12) | [[File:Vol1 page 0112 eq 001.png|RTENOTITLE]]....................(12) | ||
As the fourth property of Green’s function noted previously requires, if the outer boundary of the reservoir is impermeable, [[File:Vol1 page 0112 inline 001.png]] or at infinity, then ''G'' vanishes at the outer boundary; that is, ''G''(Γ<sub>''e''</sub>) = 0. Thus, '''Eq. 12''' becomes | As the fourth property of Green’s function noted previously requires, if the outer boundary of the reservoir is impermeable, [[File:Vol1 page 0112 inline 001.png|RTENOTITLE]] or at infinity, then ''G'' vanishes at the outer boundary; that is, ''G''(Γ<sub>''e''</sub>) = 0. Thus, '''Eq. 12''' becomes | ||
[[File:Vol1 page 0112 eq 002.png]]....................(13) | [[File:Vol1 page 0112 eq 002.png|RTENOTITLE]]....................(13) | ||
Similarly, if the flux, [[File:Vol1 page 0112 inline 002.png]], is specified at the inner boundary, then the normal derivative of Green’s function, [[File:Vol1 page 0112 inline 003.png]], vanishes at that boundary. This yields | Similarly, if the flux, [[File:Vol1 page 0112 inline 002.png|RTENOTITLE]], is specified at the inner boundary, then the normal derivative of Green’s function, [[File:Vol1 page 0112 inline 003.png|RTENOTITLE]], vanishes at that boundary. This yields | ||
[[File:Vol1 page 0112 eq 003.png]]....................(14) | [[File:Vol1 page 0112 eq 003.png|RTENOTITLE]]....................(14) | ||
If the initial pressure, ''p''<sub>''i''</sub>, is uniform over the entire domain (porous medium), ''D'', then, by the third property of Green’s function ('''Eq. 4'''), we should have | If the initial pressure, ''p''<sub>''i''</sub>, is uniform over the entire domain (porous medium), ''D'', then, by the third property of Green’s function ('''Eq. 4'''), we should have | ||
[[File:Vol1 page 0112 eq 004.png]]....................(15) | [[File:Vol1 page 0112 eq 004.png|RTENOTITLE]]....................(15) | ||
Also, the flux law for porous medium (Darcy’s law) requires that the volume of fluid passing through the point, ''M′''<sub>''w''</sub>, on the inner-boundary surface, Γ<sub>''w''</sub>, at time ''t'' be equal to | Also, the flux law for porous medium (Darcy’s law) requires that the volume of fluid passing through the point, ''M′''<sub>''w''</sub>, on the inner-boundary surface, Γ<sub>''w''</sub>, at time ''t'' be equal to | ||
[[File:Vol1 page 0112 eq 005.png]]....................(16) | [[File:Vol1 page 0112 eq 005.png|RTENOTITLE]]....................(16) | ||
The substitution of '''Eqs. 15''' and '''16''' into '''Eq. 14''' yields | The substitution of '''Eqs. 15''' and '''16''' into '''Eq. 14''' yields | ||
[[File:Vol1 page 0112 eq 006.png]]....................(17) | [[File:Vol1 page 0112 eq 006.png|RTENOTITLE]]....................(17) | ||
Not surprisingly, '''Eq. 17''' is the same as '''Eq. 18''' because ''G'' in '''Eq. 17''' is the instantaneous point-source solution of unit strength denoted by ''S'' in '''Eq. 18'''. | Not surprisingly, '''Eq. 17''' is the same as '''Eq. 18''' because ''G'' in '''Eq. 17''' is the instantaneous point-source solution of unit strength denoted by ''S'' in '''Eq. 18'''. | ||
[[File:Vol1 page 0104 eq 005.png]]....................(18) | [[File:Vol1 page 0104 eq 005.png|RTENOTITLE]]....................(18) | ||
The expression given in '''Eq. 17''' may be simplified further by assuming that the flux, [[File:Vol1 page 0112 inline 004.png]], is uniformly distributed on the inner-boundary surface (wellbore), Γ<sub>''w''</sub>. This yields | The expression given in '''Eq. 17''' may be simplified further by assuming that the flux, [[File:Vol1 page 0112 inline 004.png|RTENOTITLE]], is uniformly distributed on the inner-boundary surface (wellbore), Γ<sub>''w''</sub>. This yields | ||
[[File:Vol1 page 0113 eq 001.png]]....................(19) | [[File:Vol1 page 0113 eq 001.png|RTENOTITLE]]....................(19) | ||
where, | where, | ||
[[File:Vol1 page 0113 eq 002.png]]....................(20) | [[File:Vol1 page 0113 eq 002.png|RTENOTITLE]]....................(20) | ||
The integration in the right side of '''Eq. 20''' represents the distribution of instantaneous point sources over the length, area, or volume of the source (well), and ''S'' denotes the corresponding instantaneous source function. [[Source function solutions of the diffusion equation]] discusses the conventional derivation of the source functions starting from the basic instantaneous point-source solution. Here, we discuss the use of infinite-plane sources as the building block with Newman’s product method.<ref name="r2" /> | The integration in the right side of '''Eq. 20''' represents the distribution of instantaneous point sources over the length, area, or volume of the source (well), and ''S'' denotes the corresponding instantaneous source function. [[Source_function_solutions_of_the_diffusion_equation|Source function solutions of the diffusion equation]] discusses the conventional derivation of the source functions starting from the basic instantaneous point-source solution. Here, we discuss the use of infinite-plane sources as the building block with Newman’s product method.<ref name="r2">Newman, A.B. 1936. Heating and Cooling Rectangular and Cylindrical Solids. Ind. Eng. Chem. 28 (5): 545–548. http://dx.doi.org/10.1021/ie50317a010</ref> | ||
Newman’s product method<ref name="r2" /> may be stated for transient-flow problems in porous media as follows: <ref name="r1" /> if a well/reservoir system can be visualized as the intersection of 1D or 2D systems, then the instantaneous source or Green’s function for this well/reservoir system can be constructed by the product of the source or [[ | Newman’s product method<ref name="r2">Newman, A.B. 1936. Heating and Cooling Rectangular and Cylindrical Solids. Ind. Eng. Chem. 28 (5): 545–548. http://dx.doi.org/10.1021/ie50317a010</ref> may be stated for transient-flow problems in porous media as follows: <ref name="r1">Gringarten, A.C. and Ramey Jr., H.J. 1973. The Use of Source and Green's Functions in Solving Unsteady-Flow Problems in Reservoirs. SPE J. 13 (5): 285-296. SPE-3818-PA. http://dx.doi.org/10.2118/3818-PA</ref> if a well/reservoir system can be visualized as the intersection of 1D or 2D systems, then the instantaneous source or Green’s function for this well/reservoir system can be constructed by the product of the source or [[Green’s_function_for_solving_transient_flow_problems|Green’s functions]] for the 1D and/or 2D systems. For example, an infinite line-source well at ''x'' = ''x′'', ''y'' = ''y′'', and −∞ ≤ ''z′'' ≤ +∞ in an infinite reservoir may be visualized as the intersection of two infinite, 1D plane sources; one at ''x'' = ''x′'', −∞ ≤ ''y′'' ≤ +∞, and −∞ ≤ ''z′'' ≤ +∞, and the other at −∞ ≤ ''x′'' ≤ +∞, ''y'' = ''y′'', and −∞ ≤ ''z′'' ≤ +∞. Then, the instantaneous source function for this infinite line-source well, ''S''(''x'', ''x′'', ''y'', ''y′'', ''t'' − ''τ''), may be obtained as the product of two infinite, 1D plane sources, given by | ||
[[File:Vol1 page 0113 eq 003.png]]....................(21) | [[File:Vol1 page 0113 eq 003.png|RTENOTITLE]]....................(21) | ||
as follows | as follows | ||
[[File:Vol1 page 0113 eq 004.png]]....................(22) | [[File:Vol1 page 0113 eq 004.png|RTENOTITLE]]....................(22) | ||
As expected, this solution is the same as '''Eq. 23''', which was obtained by integration of an instantaneous point source in an infinite reservoir. For a radially isotropic reservoir (''η''<sub>''x''</sub> = ''η''<sub>''y''</sub> = ''η''<sub>''z''</sub>), '''Eq. 22''' yields | As expected, this solution is the same as '''Eq. 23''', which was obtained by integration of an instantaneous point source in an infinite reservoir. For a radially isotropic reservoir (''η''<sub>''x''</sub> = ''η''<sub>''y''</sub> = ''η''<sub>''z''</sub>), '''Eq. 22''' yields | ||
[[File:Vol1 page 0104 eq 002.png]]....................(23) | [[File:Vol1 page 0104 eq 002.png|RTENOTITLE]]....................(23) | ||
[[File:Vol1 page 0113 eq 005.png]]....................(24) | [[File:Vol1 page 0113 eq 005.png|RTENOTITLE]]....................(24) | ||
where ''d'' is the distance between the line source and the observation line in the ''x''-''y'' plane (see '''Fig. 1''') and is given by | where ''d'' is the distance between the line source and the observation line in the ''x''-''y'' plane (see '''Fig. 1''') and is given by | ||
[[File:Vol1 page 0113 eq 006.png]]....................(25) | [[File:Vol1 page 0113 eq 006.png|RTENOTITLE]]....................(25) | ||
Similarly, intersecting three infinite instantaneous plane sources (or a line source and a plane source), we can obtain the instantaneous point-source solution in an infinite reservoir as | Similarly, intersecting three infinite instantaneous plane sources (or a line source and a plane source), we can obtain the instantaneous point-source solution in an infinite reservoir as | ||
[[File:Vol1 page 0114 eq 001.png]]....................(26) | [[File:Vol1 page 0114 eq 001.png|RTENOTITLE]]....................(26) | ||
Instantaneous plane sources in slab reservoirs can be generated with the plane sources in infinite reservoirs and the method of images as discussed in [[Source function solutions of the diffusion equation]]. Similarly, the instantaneous slab sources can be obtained by integrating plane sources over the thickness of the slab source. '''Table 1''', compiled from the work of Gringarten and Ramey,<ref name="r1" /> presents the basic instantaneous source functions in infinite reservoirs, and '''Table 2''' shows the corresponding geometries of the source-reservoir systems. The basic instantaneous source functions given in '''Table 2''' may be used to construct the source functions that represent the appropriate well geometry by Newman’s product method. | Instantaneous plane sources in slab reservoirs can be generated with the plane sources in infinite reservoirs and the method of images as discussed in [[Source_function_solutions_of_the_diffusion_equation|Source function solutions of the diffusion equation]]. Similarly, the instantaneous slab sources can be obtained by integrating plane sources over the thickness of the slab source. '''Table 1''', compiled from the work of Gringarten and Ramey,<ref name="r1">Gringarten, A.C. and Ramey Jr., H.J. 1973. The Use of Source and Green's Functions in Solving Unsteady-Flow Problems in Reservoirs. SPE J. 13 (5): 285-296. SPE-3818-PA. http://dx.doi.org/10.2118/3818-PA</ref> presents the basic instantaneous source functions in infinite reservoirs, and '''Table 2''' shows the corresponding geometries of the source-reservoir systems. The basic instantaneous source functions given in '''Table 2''' may be used to construct the source functions that represent the appropriate well geometry by Newman’s product method. | ||
<gallery widths="300px" heights="200px"> | <gallery widths="300px" heights="200px"> | ||
Line 124: | Line 123: | ||
</gallery> | </gallery> | ||
Gringarten and Ramey<ref name="r1" /> have also presented the approximating forms of the instantaneous linear sources and the time ranges for these approximations to be valid. The approximate solutions are very useful in obtaining expressions for pressure distributions at early and late times and identifying the flow regimes during the corresponding time periods. '''Table 3''' presents the short- and long-time approximating forms for instantaneous linear sources and their time ranges. Examples 1 and 2 present some applications of the product-solution method and the derivation of the approximate solutions for pressure distributions. | Gringarten and Ramey<ref name="r1">Gringarten, A.C. and Ramey Jr., H.J. 1973. The Use of Source and Green's Functions in Solving Unsteady-Flow Problems in Reservoirs. SPE J. 13 (5): 285-296. SPE-3818-PA. http://dx.doi.org/10.2118/3818-PA</ref> have also presented the approximating forms of the instantaneous linear sources and the time ranges for these approximations to be valid. The approximate solutions are very useful in obtaining expressions for pressure distributions at early and late times and identifying the flow regimes during the corresponding time periods. '''Table 3''' presents the short- and long-time approximating forms for instantaneous linear sources and their time ranges. Examples 1 and 2 present some applications of the product-solution method and the derivation of the approximate solutions for pressure distributions. | ||
<gallery widths=300px heights=200px> | <gallery widths="300px" heights="200px"> | ||
File:Vol1 Page 117 Image 0001.png|'''Table 3''' | File:Vol1 Page 117 Image 0001.png|'''Table 3''' | ||
</gallery> | </gallery> | ||
==Example 1 - Ppartially penetrating vertical fracture in an infinite homogeneous slab reservoir== | == Example 1 - Ppartially penetrating vertical fracture in an infinite homogeneous slab reservoir == | ||
Considering transient flow toward a partially penetrating vertical fracture of vertical penetration ''h''<sub>''f''</sub> and horizontal penetration 2''x''<sub>''f''</sub> in an infinite, homogeneous, slab reservoir of uniform thickness, ''h'', and initial pressure, ''p''<sub>''i''</sub>, with impermeable top and bottom boundaries. | Considering transient flow toward a partially penetrating vertical fracture of vertical penetration ''h''<sub>''f''</sub> and horizontal penetration 2''x''<sub>''f''</sub> in an infinite, homogeneous, slab reservoir of uniform thickness, ''h'', and initial pressure, ''p''<sub>''i''</sub>, with impermeable top and bottom boundaries. | ||
''Solution.'' | ''Solution.'' '''Fig. 2''' shows the geometry of the well reservoir system of interest. Approximate the fracture by a vertical plane of height ''h''<sub>''f''</sub> and length 2''x''<sub>''f''</sub>. The corresponding source geometry may be visualized as the intersection of an infinite plane source at ''y'' = ''y′'' in an infinite reservoir (Source I in '''Tables 1''' and '''2'''), an infinite-slab source of thickness 2''x''<sub>''f''</sub> at ''x'' = ''x′'' in an infinite reservoir (Source IV), and an infinite-slab source of thickness ''h''<sub>''p''</sub> = ''h''<sub>''f''</sub> at ''z'' = ''z''<sub>''w''</sub> in an infinite-slab reservoir of thickness ''h'' (Source VIII). Then, by Newman’s product method, the appropriate source function is given by | ||
[[File:Vol1 page 0115 eq 001.png]] | [[File:Vol1 page 0115 eq 001.png|RTENOTITLE]] | ||
[[File:Vol1 page 0116 eq 001.png]]....................(27) | [[File:Vol1 page 0116 eq 001.png|RTENOTITLE]]....................(27) | ||
Assuming that the production is at a constant rate, [[File:Vol1 page 0116 inline 001.png|RTENOTITLE]] and using '''Eq. 27''' for the source function, ''S'', in '''Eq. 19''', we obtain | |||
[[File:Vol1 page 0116 eq 002.png|RTENOTITLE]] | |||
[[File:Vol1 page 0117 eq 001.png|RTENOTITLE]]....................(28) | |||
[[File:Vol1 page 0117 eq 001.png]]....................(28) | |||
If the fracture penetrates the entire thickness of the reservoir (i.e., ''h''<sub>''f''</sub> = ''h'') as shown in '''Fig. 3''', then '''Eq. 28''' yields | If the fracture penetrates the entire thickness of the reservoir (i.e., ''h''<sub>''f''</sub> = ''h'') as shown in '''Fig. 3''', then '''Eq. 28''' yields | ||
[[File:Vol1 page 0117 eq 002.png]]....................(29) | [[File:Vol1 page 0117 eq 002.png|RTENOTITLE]]....................(29) | ||
The fully penetrating fracture solution given in '''Eq. 29''' also could be obtained by constructing the source function as the product of an infinite plane source at ''y'' = ''y′'' in an infinite reservoir (Source I in '''Tables 1''' and '''2''') and an infinite-slab source of thickness 2''x''<sub>''f''</sub> at ''x'' = ''x′'' in an infinite reservoir (Source IV). This source function then could be used in '''Eq. 19'''. | The fully penetrating fracture solution given in '''Eq. 29''' also could be obtained by constructing the source function as the product of an infinite plane source at ''y'' = ''y′'' in an infinite reservoir (Source I in '''Tables 1''' and '''2''') and an infinite-slab source of thickness 2''x''<sub>''f''</sub> at ''x'' = ''x′'' in an infinite reservoir (Source IV). This source function then could be used in '''Eq. 19'''. | ||
Line 160: | Line 159: | ||
'''Fig. 4''' presents an example of transient-pressure responses computed from '''Eq. 29'''. To obtain the results shown in '''Fig. 4''', numerical integration has been used to evaluate the right side of '''Eq. 29'''. It is also of interest to obtain an early-time approximation for the solution given in '''Eq. 29'''. Substituting the early-time approximating forms for the slab sources in an infinite reservoir (approximations given in '''Table 3''' for Source Functions IV and VIII), we obtain | '''Fig. 4''' presents an example of transient-pressure responses computed from '''Eq. 29'''. To obtain the results shown in '''Fig. 4''', numerical integration has been used to evaluate the right side of '''Eq. 29'''. It is also of interest to obtain an early-time approximation for the solution given in '''Eq. 29'''. Substituting the early-time approximating forms for the slab sources in an infinite reservoir (approximations given in '''Table 3''' for Source Functions IV and VIII), we obtain | ||
[[File:Vol1 page 0118 eq 001.png]]....................(30) | [[File:Vol1 page 0118 eq 001.png|RTENOTITLE]]....................(30) | ||
where | where | ||
[[File:Vol1 page 0119 eq 001.png]]....................(31) | [[File:Vol1 page 0119 eq 001.png|RTENOTITLE]]....................(31) | ||
and | and | ||
[[File:Vol1 page 0119 eq 002.png]]....................(32) | [[File:Vol1 page 0119 eq 002.png|RTENOTITLE]]....................(32) | ||
<gallery widths="300px" heights="200px"> | <gallery widths="300px" heights="200px"> | ||
Line 174: | Line 173: | ||
</gallery> | </gallery> | ||
Assuming a constant production rate, [[File:Vol1 page 0116 inline 001.png]], and substituting the source function given by '''Eq. 30''' in '''Eq. 19''', we obtain | Assuming a constant production rate, [[File:Vol1 page 0116 inline 001.png|RTENOTITLE]], and substituting the source function given by '''Eq. 30''' in '''Eq. 19''', we obtain | ||
[[File:Vol1 page 0119 eq 003.png]]....................(33) | [[File:Vol1 page 0119 eq 003.png|RTENOTITLE]]....................(33) | ||
where erfc (''z'') is the complementary error function defined by | where erfc (''z'') is the complementary error function defined by | ||
[[File:Vol1 page 0119 eq 004.png]]....................(34) | [[File:Vol1 page 0119 eq 004.png|RTENOTITLE]]....................(34) | ||
== Example 2 - Uniform-flux horizontal well in a closed, homogeneous, rectangular parallelepiped reservoir == | |||
Consider transient flow toward a uniform-flux horizontal well of length ''L''<sub>''h''</sub> located at (''x′'', ''y′'', ''z''<sub>''w''</sub>) in a closed, homogeneous, rectangular parallelepiped of dimensions 0 ≤ ''x'' ≤ ''x''<sub>''e''</sub>, 0 ≤ ''y'' ≤ ''y''<sub>''e''</sub>, 0 ≤ ''z'' ≤ ''h'' and of initial pressure, ''p''<sub>''i''</sub>. | Consider transient flow toward a uniform-flux horizontal well of length ''L''<sub>''h''</sub> located at (''x′'', ''y′'', ''z''<sub>''w''</sub>) in a closed, homogeneous, rectangular parallelepiped of dimensions 0 ≤ ''x'' ≤ ''x''<sub>''e''</sub>, 0 ≤ ''y'' ≤ ''y''<sub>''e''</sub>, 0 ≤ ''z'' ≤ ''h'' and of initial pressure, ''p''<sub>''i''</sub>. | ||
''Solution.'' | ''Solution.'' '''Fig. 5''' shows the sketch of the horizontal well/reservoir system considered in this example. If we approximate the horizontal well by a horizontal line source of length ''L''<sub>''h''</sub>, then the resulting source/reservoir system may be visualized as the intersection of three sources: an infinite plane source at ''y'' = ''y′'' in an infinite-slab reservoir of thickness ye with impermeable boundaries (Source V in '''Tables 1''' and '''2'''), an infinite-slab source of thickness ''L''<sub>''h''</sub> at ''x'' = ''x′'' in an infinite-slab reservoir of thickness ''x''<sub>''e''</sub> with impermeable boundaries (Source VIII), and an infinite-plane source at ''z'' = ''z''<sub>''w''</sub> in an infinite-slab reservoir of thickness ''h'' with impermeable boundaries (Source V). Then, by Newman’s product method, the appropriate source function can be obtained as | ||
[[File:Vol1 page 0120 eq 001.png]]....................(35) | [[File:Vol1 page 0120 eq 001.png|RTENOTITLE]]....................(35) | ||
Assuming that the production is at a constant rate, [[File:Vol1 page 0120 inline 001.png]], and using '''Eq. 35''' for the source function, ''S'', in '''Eq. 19''', we obtain | Assuming that the production is at a constant rate, [[File:Vol1 page 0120 inline 001.png|RTENOTITLE]], and using '''Eq. 35''' for the source function, ''S'', in '''Eq. 19''', we obtain | ||
[[File:Vol1 page 0120 eq 002.png]]....................(36) | [[File:Vol1 page 0120 eq 002.png|RTENOTITLE]]....................(36) | ||
'''Table 4''' presents the pressure responses for a uniform-flux horizontal well in a closed square computed from '''Eq. 36'''. We may obtain a short-time approximation for '''Eq. 3.216''' with the early-time approximations given in '''Table 3''' for Source Functions V and VIII. This yields | '''Table 4''' presents the pressure responses for a uniform-flux horizontal well in a closed square computed from '''Eq. 36'''. We may obtain a short-time approximation for '''Eq. 3.216''' with the early-time approximations given in '''Table 3''' for Source Functions V and VIII. This yields | ||
[[File:Vol1 page 0120 eq 003.png]]....................(37) | [[File:Vol1 page 0120 eq 003.png|RTENOTITLE]]....................(37) | ||
where Ei(−''x'') is the exponential integral function defined by '''Eq. 38'''. '''Eq. 37''' indicates that the early-time flow is radial in the ''y''-''z'' plane around the axis of the horizontal well. This solution corresponds to the time period during which none of the reservoir boundaries influence the pressure response. | where Ei(−''x'') is the exponential integral function defined by '''Eq. 38'''. '''Eq. 37''' indicates that the early-time flow is radial in the ''y''-''z'' plane around the axis of the horizontal well. This solution corresponds to the time period during which none of the reservoir boundaries influence the pressure response. | ||
[[File:Vol1 page 0092 eq 005.png]]....................(38) | [[File:Vol1 page 0092 eq 005.png|RTENOTITLE]]....................(38) | ||
<gallery widths="300px" heights="200px"> | <gallery widths="300px" heights="200px"> | ||
Line 209: | Line 209: | ||
It is also possible to obtain another approximation for '''Eq. 36''' that covers the intermediate time-flow behavior. If we approximate the source function in the ''x'' direction (Source Function VIII) by its early and intermediate time approximation and the source function in the ''y'' direction (Source Function V) by its early time approximation given in '''Table 3''', we obtain | It is also possible to obtain another approximation for '''Eq. 36''' that covers the intermediate time-flow behavior. If we approximate the source function in the ''x'' direction (Source Function VIII) by its early and intermediate time approximation and the source function in the ''y'' direction (Source Function V) by its early time approximation given in '''Table 3''', we obtain | ||
[[File:Vol1 page 0121 eq 001.png]]....................(39) | [[File:Vol1 page 0121 eq 001.png|RTENOTITLE]]....................(39) | ||
This approximation should correspond to the time period during which the influence of the top and/or bottom boundaries may be evident but the lateral boundaries in the ''x'' and ''y'' directions do not have an influence on the pressure response. This solution also could have been obtained by assuming a laterally infinite reservoir. In this case, the source function would have been constructed as the product of three source functions: an infinite-plane source at ''y'' = ''y′'' in an infinite reservoir (Source I in '''Tables 1''' and '''2'''), an infinite-slab source of thickness ''L''<sub>''h''</sub> at ''x'' = ''x′'' in an infinite reservoir (Source IV), and an infinite-plane source at ''z'' = ''z''<sub>''w''</sub> in an infinite-slab reservoir of thickness ''h'' with impermeable boundaries (Source V). | This approximation should correspond to the time period during which the influence of the top and/or bottom boundaries may be evident but the lateral boundaries in the ''x'' and ''y'' directions do not have an influence on the pressure response. This solution also could have been obtained by assuming a laterally infinite reservoir. In this case, the source function would have been constructed as the product of three source functions: an infinite-plane source at ''y'' = ''y′'' in an infinite reservoir (Source I in '''Tables 1''' and '''2'''), an infinite-slab source of thickness ''L''<sub>''h''</sub> at ''x'' = ''x′'' in an infinite reservoir (Source IV), and an infinite-plane source at ''z'' = ''z''<sub>''w''</sub> in an infinite-slab reservoir of thickness ''h'' with impermeable boundaries (Source V). | ||
==Nomenclature== | == Nomenclature == | ||
{| | {| | ||
|- | |- | ||
|'' | | ''c'' | ||
|= | | = | ||
| | | fluid compressibility, atm<sup>−1</sup> | ||
|- | |- | ||
|'' | | ''d'' | ||
|= | | = | ||
| | | distance to a linear boundary, cm | ||
|- | |- | ||
|'' | | ''D'' | ||
|= | | = | ||
| | | domain | ||
|- | |- | ||
|'' | | ''Ei''(''x'') | ||
|= | | = | ||
| | | exponential integral function | ||
|- | |- | ||
|'' | | ''G'' | ||
|= | | = | ||
| | | Green’s function | ||
|- | |- | ||
|''h'' | | ''h'' | ||
|= | | = | ||
| | | formation thickness, cm | ||
|- | |- | ||
|'' | | ''h''<sub>''f''</sub> | ||
|= | | = | ||
| | | fracture height (vertical penetration), cm | ||
|- | |- | ||
|'' | | ''k'' | ||
|= | | = | ||
| | | isotropic permeability, md | ||
|- | |- | ||
|'' | | ''L''<sub>''h''</sub> | ||
|= | | = | ||
| | | horizontal-well length, cm | ||
|- | |- | ||
|'' | | ''M'' | ||
|= | | = | ||
| | | point in space | ||
|- | |- | ||
|'' | | ''M′'' | ||
|= | | = | ||
|point in | | source point in space | ||
|- | |- | ||
|'' | | ''M''<sub>''w''</sub> | ||
|= | | = | ||
| | | point in Γ<sub>''w''</sub> | ||
|- | |- | ||
|'' | | ''M′''<sub>''w''</sub> | ||
|= | | = | ||
| | | source point in Γ<sub>''w''</sub> | ||
|- | |- | ||
|''p'' | | ''p'' | ||
|= | | = | ||
| | | pressure, atm | ||
|- | |- | ||
|'' | | ''p''<sub>''i''</sub> | ||
|= | | = | ||
| | | initial pressure, atm | ||
|- | |- | ||
| | | ''q'' | ||
|= | | = | ||
| | | production rate, cm<sup>3</sup>/s | ||
|- | |- | ||
| | | [[File:Vol1 page 0102 inline 001.png|RTENOTITLE]] | ||
|= | | = | ||
| | | instantaneous production rate for a point source, cm<sup>3</sup>/s | ||
|- | |- | ||
|'' | | ''r'' | ||
|= | | = | ||
| | | radial coordinate and distance, cm | ||
|- | |- | ||
|'' | | ''r′'' | ||
|= | | = | ||
|source | | source coordinate in ''r''-direction, cm | ||
|- | |- | ||
|'' | | ''S'' | ||
|= | | = | ||
| | | source function | ||
|- | |- | ||
|'' | | ''t'' | ||
|= | | = | ||
| | | time, s | ||
|- | |- | ||
|'' | | ''x'' | ||
|= | | = | ||
| | | distance in ''x''-direction, cm | ||
|- | |- | ||
|'' | | ''x′'' | ||
|= | | = | ||
| | | source coordinate in ''x''-direction, cm | ||
|- | |- | ||
|''x''<sub>'' | | ''x''<sub>''e''</sub> | ||
|= | | = | ||
| | | distance to the external boundary in ''x''-direction, cm | ||
|- | |- | ||
|'' | | ''x''<sub>''f''</sub> | ||
|= | | = | ||
| | | fracture half-length, cm | ||
|- | |- | ||
|'' | | ''y'' | ||
|= | | = | ||
| | | distance in ''y''-direction, cm | ||
|- | |- | ||
|'' | | ''y′'' | ||
|= | | = | ||
| | | source coordinate in ''y''-direction, cm | ||
|- | |- | ||
|'' | | ''y''<sub>''e''</sub> | ||
|= | | = | ||
|distance in '' | | distance to the external boundary in ''y''-direction, cm | ||
|- | |- | ||
|'' | | ''z'' | ||
|= | | = | ||
| | | distance in ''z''-direction, cm | ||
|- | |- | ||
|'' | | ''z′'' | ||
|= | | = | ||
| | | source coordinate in ''z''-direction, cm | ||
|- | |- | ||
|'' | | ''z''<sub>''w''</sub> | ||
|= | | = | ||
| | | well coordinate in ''z''-direction, cm | ||
|- | |- | ||
|'' | | ''α'' | ||
|= | | = | ||
|permeability direction, | | permeability direction, | ||
|- | |- | ||
| | | ''β'' | ||
|= | | = | ||
| | | permeability direction, | ||
|- | |- | ||
|Γ | | Γ | ||
|= | | = | ||
| | | boundary surface, cm<sup>2</sup> | ||
|- | |- | ||
|Γ<sub>'' | | Γ<sub>''e''</sub> | ||
|= | | = | ||
| | | external boundary surface | ||
|- | |- | ||
| | | Γ<sub>''w''</sub> | ||
|= | | = | ||
| | | length, surface, or volume of the source | ||
|- | |- | ||
| | | Δ | ||
|= | | = | ||
| | | difference operator | ||
|- | |- | ||
|''η'' | | ''η'' | ||
|= | | = | ||
|diffusivity constant | | diffusivity constant | ||
|- | |- | ||
|'' | | ''η''<sub>''i''</sub> | ||
|= | | = | ||
| | | diffusivity constant in i direction, ''i'' = ''x'', ''y'', ''z'', or ''r'' | ||
|- | |- | ||
|'' | | ''θ'' | ||
|= | | = | ||
| | | angle from positive x-direction, degrees radian | ||
|- | |- | ||
|'' | | ''θ′'' | ||
|= | | = | ||
| | | source coordinate in ''θ''-direction, degrees radian | ||
|- | |- | ||
|'' | | ''μ'' | ||
|= | | = | ||
| | | viscosity, cp | ||
|- | |- | ||
|'' | | ''τ'' | ||
|= | | = | ||
| | | time, s | ||
|- | |- | ||
|''φ''(''M'') | | ''Φ'' | ||
|= | | = | ||
|any continuous function | | porosity, fraction | ||
|- | |||
| ''φ''(''M'') | |||
| = | |||
| any continuous function | |||
|} | |} | ||
==References== | == References == | ||
< | <references /> | ||
== Noteworthy papers in OnePetro == | |||
Use this section to list papers in OnePetro that a reader who wants to learn more should definitely read | Use this section to list papers in OnePetro that a reader who wants to learn more should definitely read | ||
==External links== | == External links == | ||
Use this section to provide links to relevant material on websites other than PetroWiki and OnePetro | Use this section to provide links to relevant material on websites other than PetroWiki and OnePetro | ||
==See also== | == See also == | ||
[[Transient analysis mathematics]] | |||
[[Transient_analysis_mathematics|Transient analysis mathematics]] | |||
[[Green’s_function_for_solving_transient_flow_problems|Green’s function for solving transient flow problems]] | |||
[[ | [[Source_function_solutions_of_the_diffusion_equation|Source function solutions of the diffusion equation]] | ||
[[ | [[Solving_unsteady_flow_problems_with_Laplace_transform_and_source_functions|Solving unsteady flow problems with Laplace transform and source functions]] | ||
[[ | [[Mathematics_of_fluid_flow|Mathematics of fluid flow]] | ||
[[ | [[Differential_calculus_refresher|Differential calculus refresher]] | ||
[[ | [[PEH:Mathematics_of_Transient_Analysis]] | ||
[[ | [[Category:5.3.1 Flow in porous media]] |
Latest revision as of 10:29, 8 June 2015
As discussed in Source function solutions of the diffusion equation, the conventional development of the source function solutions uses the instantaneous point-source solution as the building block with the appropriate integration (superposition) in space and time. In 1973, Gringarten and Ramey[1] introduced the use of the source and Green’s function method to the petroleum engineering literature with a more efficient method of developing the source solutions. Specifically, they suggested the use of infinite-plane sources as the building block with Newman’s product method.[2] In this page we discuss the use of Green’s functions and source functions in solving unsteady-flow problems in reservoirs.
Green's functions and source functions in solving unsteady flow problems
Green’s function for transient flow in a porous medium is defined as the pressure at M (x, y, z) at time t because of an instantaneous point source of unit strength generated at point M′(x′, y′, z′) at time τ < t with the porous medium initially at zero pressure and the boundary of the medium kept at zero pressure or impermeable to flow.[1][3] If we let G(M, M′, t − τ) denote the Green’s function, then it should satisfy the diffusion equation; that is,
Because G is a function of t − τ, it should also satisfy the adjoint diffusion equation,
Green’s function also has the following properties: [1][3]
1. G is symmetrical in the two points M and M′; that is, Green’s function is invariant as the source and the observation points are interchanged. 2. As t → τ, G vanishes at all points in the porous medium; that is, , except at the source location, M = M′, where it becomes infinite, so that G satisfies the delta function property,
where D indicates the domain of the porous medium, and φ(M) is any continuous function. 3. Because G corresponds to the pressure because of an instantaneous point source of unit strength, it satisfies
4. G or its normal derivative, ∂G/∂n, vanishes at the boundary, Γ, of the porous medium. If the porous medium is infinite, then G vanishes as M or M′→∞.
Let p(M′ , τ) be the pressure in the porous medium and G(M, M′, t - τ) be the Green’s function. Let D denote the domain of the porous medium. Then, p and G satisfy the following differential equations:
and
Then, we can write
or
where ε is a small positive number. Changing the order of integration and applying the Green’s theorem,
where D and Γ indicate the volume and boundary surface of the domain, respectively; S denotes the points on the boundary; and ∂/∂n indicates differentiation in the normal direction of the surface Γ; we obtain
Taking the limit as ε→0 and noting the delta-function property of the Green’s function (Eq. 3), Eq. 10 yields
where pi(M) = p(M, t = 0) is the initial pressure at point M.
In Eq. 11, the boundary of the porous medium consists of two surfaces: the inner boundary that corresponds to the surface of the wellbore, Γw, and the outer boundary of the reservoir, Γe. Eq. 11 may be written as
As the fourth property of Green’s function noted previously requires, if the outer boundary of the reservoir is impermeable, or at infinity, then G vanishes at the outer boundary; that is, G(Γe) = 0. Thus, Eq. 12 becomes
Similarly, if the flux, , is specified at the inner boundary, then the normal derivative of Green’s function, , vanishes at that boundary. This yields
If the initial pressure, pi, is uniform over the entire domain (porous medium), D, then, by the third property of Green’s function (Eq. 4), we should have
Also, the flux law for porous medium (Darcy’s law) requires that the volume of fluid passing through the point, M′w, on the inner-boundary surface, Γw, at time t be equal to
The substitution of Eqs. 15 and 16 into Eq. 14 yields
Not surprisingly, Eq. 17 is the same as Eq. 18 because G in Eq. 17 is the instantaneous point-source solution of unit strength denoted by S in Eq. 18.
The expression given in Eq. 17 may be simplified further by assuming that the flux, , is uniformly distributed on the inner-boundary surface (wellbore), Γw. This yields
where,
The integration in the right side of Eq. 20 represents the distribution of instantaneous point sources over the length, area, or volume of the source (well), and S denotes the corresponding instantaneous source function. Source function solutions of the diffusion equation discusses the conventional derivation of the source functions starting from the basic instantaneous point-source solution. Here, we discuss the use of infinite-plane sources as the building block with Newman’s product method.[2]
Newman’s product method[2] may be stated for transient-flow problems in porous media as follows: [1] if a well/reservoir system can be visualized as the intersection of 1D or 2D systems, then the instantaneous source or Green’s function for this well/reservoir system can be constructed by the product of the source or Green’s functions for the 1D and/or 2D systems. For example, an infinite line-source well at x = x′, y = y′, and −∞ ≤ z′ ≤ +∞ in an infinite reservoir may be visualized as the intersection of two infinite, 1D plane sources; one at x = x′, −∞ ≤ y′ ≤ +∞, and −∞ ≤ z′ ≤ +∞, and the other at −∞ ≤ x′ ≤ +∞, y = y′, and −∞ ≤ z′ ≤ +∞. Then, the instantaneous source function for this infinite line-source well, S(x, x′, y, y′, t − τ), may be obtained as the product of two infinite, 1D plane sources, given by
as follows
As expected, this solution is the same as Eq. 23, which was obtained by integration of an instantaneous point source in an infinite reservoir. For a radially isotropic reservoir (ηx = ηy = ηz), Eq. 22 yields
where d is the distance between the line source and the observation line in the x-y plane (see Fig. 1) and is given by
Similarly, intersecting three infinite instantaneous plane sources (or a line source and a plane source), we can obtain the instantaneous point-source solution in an infinite reservoir as
Instantaneous plane sources in slab reservoirs can be generated with the plane sources in infinite reservoirs and the method of images as discussed in Source function solutions of the diffusion equation. Similarly, the instantaneous slab sources can be obtained by integrating plane sources over the thickness of the slab source. Table 1, compiled from the work of Gringarten and Ramey,[1] presents the basic instantaneous source functions in infinite reservoirs, and Table 2 shows the corresponding geometries of the source-reservoir systems. The basic instantaneous source functions given in Table 2 may be used to construct the source functions that represent the appropriate well geometry by Newman’s product method.
Gringarten and Ramey[1] have also presented the approximating forms of the instantaneous linear sources and the time ranges for these approximations to be valid. The approximate solutions are very useful in obtaining expressions for pressure distributions at early and late times and identifying the flow regimes during the corresponding time periods. Table 3 presents the short- and long-time approximating forms for instantaneous linear sources and their time ranges. Examples 1 and 2 present some applications of the product-solution method and the derivation of the approximate solutions for pressure distributions.
Example 1 - Ppartially penetrating vertical fracture in an infinite homogeneous slab reservoir
Considering transient flow toward a partially penetrating vertical fracture of vertical penetration hf and horizontal penetration 2xf in an infinite, homogeneous, slab reservoir of uniform thickness, h, and initial pressure, pi, with impermeable top and bottom boundaries.
Solution. Fig. 2 shows the geometry of the well reservoir system of interest. Approximate the fracture by a vertical plane of height hf and length 2xf. The corresponding source geometry may be visualized as the intersection of an infinite plane source at y = y′ in an infinite reservoir (Source I in Tables 1 and 2), an infinite-slab source of thickness 2xf at x = x′ in an infinite reservoir (Source IV), and an infinite-slab source of thickness hp = hf at z = zw in an infinite-slab reservoir of thickness h (Source VIII). Then, by Newman’s product method, the appropriate source function is given by
Assuming that the production is at a constant rate, and using Eq. 27 for the source function, S, in Eq. 19, we obtain
If the fracture penetrates the entire thickness of the reservoir (i.e., hf = h) as shown in Fig. 3, then Eq. 28 yields
The fully penetrating fracture solution given in Eq. 29 also could be obtained by constructing the source function as the product of an infinite plane source at y = y′ in an infinite reservoir (Source I in Tables 1 and 2) and an infinite-slab source of thickness 2xf at x = x′ in an infinite reservoir (Source IV). This source function then could be used in Eq. 19.
Fig. 4 presents an example of transient-pressure responses computed from Eq. 29. To obtain the results shown in Fig. 4, numerical integration has been used to evaluate the right side of Eq. 29. It is also of interest to obtain an early-time approximation for the solution given in Eq. 29. Substituting the early-time approximating forms for the slab sources in an infinite reservoir (approximations given in Table 3 for Source Functions IV and VIII), we obtain
where
and
Assuming a constant production rate, , and substituting the source function given by Eq. 30 in Eq. 19, we obtain
where erfc (z) is the complementary error function defined by
Example 2 - Uniform-flux horizontal well in a closed, homogeneous, rectangular parallelepiped reservoir
Consider transient flow toward a uniform-flux horizontal well of length Lh located at (x′, y′, zw) in a closed, homogeneous, rectangular parallelepiped of dimensions 0 ≤ x ≤ xe, 0 ≤ y ≤ ye, 0 ≤ z ≤ h and of initial pressure, pi.
Solution. Fig. 5 shows the sketch of the horizontal well/reservoir system considered in this example. If we approximate the horizontal well by a horizontal line source of length Lh, then the resulting source/reservoir system may be visualized as the intersection of three sources: an infinite plane source at y = y′ in an infinite-slab reservoir of thickness ye with impermeable boundaries (Source V in Tables 1 and 2), an infinite-slab source of thickness Lh at x = x′ in an infinite-slab reservoir of thickness xe with impermeable boundaries (Source VIII), and an infinite-plane source at z = zw in an infinite-slab reservoir of thickness h with impermeable boundaries (Source V). Then, by Newman’s product method, the appropriate source function can be obtained as
Assuming that the production is at a constant rate, , and using Eq. 35 for the source function, S, in Eq. 19, we obtain
Table 4 presents the pressure responses for a uniform-flux horizontal well in a closed square computed from Eq. 36. We may obtain a short-time approximation for Eq. 3.216 with the early-time approximations given in Table 3 for Source Functions V and VIII. This yields
where Ei(−x) is the exponential integral function defined by Eq. 38. Eq. 37 indicates that the early-time flow is radial in the y-z plane around the axis of the horizontal well. This solution corresponds to the time period during which none of the reservoir boundaries influence the pressure response.
It is also possible to obtain another approximation for Eq. 36 that covers the intermediate time-flow behavior. If we approximate the source function in the x direction (Source Function VIII) by its early and intermediate time approximation and the source function in the y direction (Source Function V) by its early time approximation given in Table 3, we obtain
This approximation should correspond to the time period during which the influence of the top and/or bottom boundaries may be evident but the lateral boundaries in the x and y directions do not have an influence on the pressure response. This solution also could have been obtained by assuming a laterally infinite reservoir. In this case, the source function would have been constructed as the product of three source functions: an infinite-plane source at y = y′ in an infinite reservoir (Source I in Tables 1 and 2), an infinite-slab source of thickness Lh at x = x′ in an infinite reservoir (Source IV), and an infinite-plane source at z = zw in an infinite-slab reservoir of thickness h with impermeable boundaries (Source V).
Nomenclature
References
- ↑ 1.0 1.1 1.2 1.3 1.4 1.5 Gringarten, A.C. and Ramey Jr., H.J. 1973. The Use of Source and Green's Functions in Solving Unsteady-Flow Problems in Reservoirs. SPE J. 13 (5): 285-296. SPE-3818-PA. http://dx.doi.org/10.2118/3818-PA
- ↑ 2.0 2.1 2.2 Newman, A.B. 1936. Heating and Cooling Rectangular and Cylindrical Solids. Ind. Eng. Chem. 28 (5): 545–548. http://dx.doi.org/10.1021/ie50317a010
- ↑ 3.0 3.1 Carslaw, H.S. and Jaeger, J.C. 1959. Conduction of Heat in Solids, second edition, 353–386. Oxford, UK: Oxford University Press.
Noteworthy papers in OnePetro
Use this section to list papers in OnePetro that a reader who wants to learn more should definitely read
External links
Use this section to provide links to relevant material on websites other than PetroWiki and OnePetro
See also
Transient analysis mathematics
Green’s function for solving transient flow problems
Source function solutions of the diffusion equation
Solving unsteady flow problems with Laplace transform and source functions