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


Difference between revisions of "Fluid flow in horizontal wells"

PetroWiki
Jump to navigation Jump to search
 
 
Line 1: Line 1:
[[Estimating horizontal well productivity|Productivity estimates in horizontal wells]] are subject to more uncertainty than comparable estimates in vertical wells. Further, it is much more difficult to interpret well test data because of 3D flow geometry. The radial symmetry usually present in a vertical well does not exist. Several flow regimes can potentially occur and need to be considered in analyzing test data from horizontal wells. Wellbore storage effects can be much more significant and partial penetration and end effects commonly complicate interpretation.  
+
[[Estimating_horizontal_well_productivity|Productivity estimates in horizontal wells]] are subject to more uncertainty than comparable estimates in vertical wells. Further, it is much more difficult to interpret well test data because of 3D flow geometry. The radial symmetry usually present in a vertical well does not exist. Several flow regimes can potentially occur and need to be considered in analyzing test data from horizontal wells. Wellbore storage effects can be much more significant and partial penetration and end effects commonly complicate interpretation.
  
In vertical wells, variables such as average permeability, net vertical thickness, and skin are used. Horizontal wells need more detail. Not only is vertical thickness important, but the horizontal dimensions of the reservoir, relative to the horizontal wellbore, need to be known.  
+
In vertical wells, variables such as average permeability, net vertical thickness, and skin are used. Horizontal wells need more detail. Not only is vertical thickness important, but the horizontal dimensions of the reservoir, relative to the horizontal wellbore, need to be known.
  
 
== Steps in evaluating horizontal well-test data ==
 
== Steps in evaluating horizontal well-test data ==
 +
 
There are three basic steps in evaluating pressure-transient data from a horizontal well:
 
There are three basic steps in evaluating pressure-transient data from a horizontal well:
 +
 
#Identify the specific flow regimes in the test data
 
#Identify the specific flow regimes in the test data
 
#Apply the proper analytical and graphical procedures to the data
 
#Apply the proper analytical and graphical procedures to the data
#Evaluate the uniqueness and sensitivity of the results to properties derived from analysis or simply assumed  
+
#Evaluate the uniqueness and sensitivity of the results to properties derived from analysis or simply assumed
  
===Identify flow regimes===
+
=== Identify flow regimes ===
Evaluation of data from a vertical wellbore will generally center on a single flow regime, such as infinite-acting radial flow, known as the MTR. However, a pressure-transient test in a horizontal well can involve as many as five major and distinct regimes that need to be identified. These regimes may or may not occur in a given test and may or may not be obscured by wellbore storage effects.
 
  
===Apply the proper procedures===
+
Evaluation of data from a vertical wellbore will generally center on a single flow regime, such as infinite-acting radial flow, known as the MTR. However, a pressure-transient test in a horizontal well can involve as many as five major and distinct regimes that need to be identified. These regimes may or may not occur in a given test and may or may not be obscured by wellbore storage effects.
Each flow regime can be modeled by an equation that can be used to estimate important reservoir properties. At best, only groups of analytical parameters can be determined directly from equations. It is imperative that the proper analytical and graphical procedures be applied to the data. In many cases, when solving for specific parameters, the application of these analytical expressions may involve a complex iterative procedure.  
 
  
===Evaluate uniqueness and sensitivity===
+
=== Apply the proper procedures ===
Experience indicates that results of horizontal well test analysis are seldom unique, so it is important that the uniqueness and sensitivity of the results to assumed properties be evaluated. Simulation of the test using properties that have been determined from the test can confirm that at least the analysis is consistent with the test data. A simulator can also determine whether other sets of formation properties will also lead to a fit of the data.  
+
 
 +
Each flow regime can be modeled by an equation that can be used to estimate important reservoir properties. At best, only groups of analytical parameters can be determined directly from equations. It is imperative that the proper analytical and graphical procedures be applied to the data. In many cases, when solving for specific parameters, the application of these analytical expressions may involve a complex iterative procedure.
 +
 
 +
=== Evaluate uniqueness and sensitivity ===
 +
 
 +
Experience indicates that results of horizontal well test analysis are seldom unique, so it is important that the uniqueness and sensitivity of the results to assumed properties be evaluated. Simulation of the test using properties that have been determined from the test can confirm that at least the analysis is consistent with the test data. A simulator can also determine whether other sets of formation properties will also lead to a fit of the data.
  
 
== Horizontal well flow regimes ==
 
== Horizontal well flow regimes ==
Different formation properties can be calculated from the data in each of the five different flow regimes. Any flow regime may be absent from a [[Diagnostic plots|plot]] of test data because of geometry, wellbore storage, or other factors. Nor does the fact that they can appear mean that they do appear. The five different flow regimes that can occur are:
+
 
 +
Different formation properties can be calculated from the data in each of the five different flow regimes. Any flow regime may be absent from a [[Diagnostic_plots|plot]] of test data because of geometry, wellbore storage, or other factors. Nor does the fact that they can appear mean that they do appear. The five different flow regimes that can occur are:
 +
 
 
*Early radial
 
*Early radial
 
*Hemi-radial
 
*Hemi-radial
 
*Early linear
 
*Early linear
 
*Late pseudoradial
 
*Late pseudoradial
*Late linear  
+
*Late linear
  
'''Fig. 1''' shows a horizontal well with length, ''L''<sub>''w''</sub>, within a reservoir that is assumed to be a rectangular parallelepiped or a "box reservoir" drainage area. In this discussion, it is assumed that the axes of the coordinate system coincide with the direction of principal permeability and the well produces over its entire length, ''L''<sub>''w''</sub>.  
+
'''Fig. 1''' shows a horizontal well with length, ''L''<sub>''w''</sub>, within a reservoir that is assumed to be a rectangular parallelepiped or a "box reservoir" drainage area. In this discussion, it is assumed that the axes of the coordinate system coincide with the direction of principal permeability and the well produces over its entire length, ''L''<sub>''w''</sub>.
  
 
<gallery widths="300px" heights="200px">
 
<gallery widths="300px" heights="200px">
Line 32: Line 39:
 
</gallery>
 
</gallery>
  
The axes for this box are the usual ''x''-, ''y''-, and ''z''-axes. Notice that the x-axis is measured along the bottom edge of the reservoir, going from left to right in the direction perpendicular to the well. The ''y''-axis lies along the axis from front to back of the reservoir, parallel to the wellbore. The ''z''-axis is oriented in the direction of reservoir thickness.  
+
The axes for this box are the usual ''x''-, ''y''-, and ''z''-axes. Notice that the x-axis is measured along the bottom edge of the reservoir, going from left to right in the direction perpendicular to the well. The ''y''-axis lies along the axis from front to back of the reservoir, parallel to the wellbore. The ''z''-axis is oriented in the direction of reservoir thickness.
 +
 
 +
The total width of the reservoir perpendicular to the wellbore is ''a''<sub>''H''</sub>, the total length in a direction parallel to the wellbore is ''b''<sub>''H''</sub>, and the total height of the reservoir is the net pay thickness, ''h''. Notice the parameters for the distance from the well to the various borders. Along the axis of the well, the shortest distance from the end of the well to a boundary is ''d''<sub>''y''</sub>, and the longest distance from the other end of the well to the boundary is ''D''<sub>''y''</sub>. In the vertical direction, the shortest distance to a vertical boundary is ''d''<sub>''z''</sub>, and the longest distance to a vertical boundary is ''D''<sub>''z''</sub>.
  
The total width of the reservoir perpendicular to the wellbore is ''a''<sub>''H''</sub>, the total length in a direction parallel to the wellbore is ''b''<sub>''H''</sub>, and the total height of the reservoir is the net pay thickness, ''h''. Notice the parameters for the distance from the well to the various borders. Along the axis of the well, the shortest distance from the end of the well to a boundary is ''d''<sub>''y''</sub>, and the longest distance from the other end of the well to the boundary is ''D''<sub>''y''</sub>. In the vertical direction, the shortest distance to a vertical boundary is ''d''<sub>''z''</sub>, and the longest distance to a vertical boundary is ''D''<sub>''z''</sub>.
+
=== Characteristics of flow regimes ===
  
===Characteristics of flow regimes===
+
Consider a well producing at a constant rate. The early radial flow regime occurs before the area drained or the pressure transient caused by this production encounters either of the vertical boundaries of the reservoir. '''Fig. 2''' shows a radial flow pattern penetrating out into the reservoir. Actually, however, this flow pattern is likely to be elliptical, moving further into the reservoir at a given time in the higher-permeability ''x''-direction than in the lower-permeability ''z''-direction. This phenomenon causes no significant complications in our analysis.
Consider a well producing at a constant rate. The early radial flow regime occurs before the area drained or the pressure transient caused by this production encounters either of the vertical boundaries of the reservoir. '''Fig. 2''' shows a radial flow pattern penetrating out into the reservoir. Actually, however, this flow pattern is likely to be elliptical, moving further into the reservoir at a given time in the higher-permeability ''x''-direction than in the lower-permeability ''z''-direction. This phenomenon causes no significant complications in our analysis.  
 
  
 
<gallery widths="300px" heights="200px">
 
<gallery widths="300px" heights="200px">
Line 43: Line 51:
 
</gallery>
 
</gallery>
  
When the wellbore is much nearer one vertical boundary than the other, another flow regime, called half-radial or hemiradial flow ('''Fig. 3''') may exist. Hemiradial flow can occur immediately following the early radial flow regime, if the well is much nearer one of the vertical boundaries than the other. Eventually, the area affected by the production will include the entire thickness of the reservoir. When that happens, a linear flow pattern may develop, as '''Fig. 4''' shows.  
+
When the wellbore is much nearer one vertical boundary than the other, another flow regime, called half-radial or hemiradial flow ('''Fig. 3''') may exist. Hemiradial flow can occur immediately following the early radial flow regime, if the well is much nearer one of the vertical boundaries than the other. Eventually, the area affected by the production will include the entire thickness of the reservoir. When that happens, a linear flow pattern may develop, as '''Fig. 4''' shows.
  
 
<gallery widths="300px" heights="200px">
 
<gallery widths="300px" heights="200px">
Line 51: Line 59:
 
</gallery>
 
</gallery>
  
Eventually, flow will begin to come into the wellbore from beyond the ends of the well. Until these end effects become important, early linear flow continues. Once end effects become important a transition period is followed by a later pseudoradial flow regime, as '''Fig. 5''' illustrates.  
+
Eventually, flow will begin to come into the wellbore from beyond the ends of the well. Until these end effects become important, early linear flow continues. Once end effects become important a transition period is followed by a later pseudoradial flow regime, as '''Fig. 5''' illustrates.
  
 
<gallery widths="300px" heights="200px">
 
<gallery widths="300px" heights="200px">
Line 57: Line 65:
 
</gallery>
 
</gallery>
  
This flow regime continues until the area affected by the production reaches one of the sides of the reservoir. Once the area affected is the entire width of the reservoir (that is, the pressure transient has reached both sides of the reservoir), then the late-linear flow regime begins ('''Fig. 6''').  
+
This flow regime continues until the area affected by the production reaches one of the sides of the reservoir. Once the area affected is the entire width of the reservoir (that is, the pressure transient has reached both sides of the reservoir), then the late-linear flow regime begins ('''Fig. 6''').
  
 
<gallery widths="300px" heights="200px">
 
<gallery widths="300px" heights="200px">
Line 64: Line 72:
  
 
== Identifying flow regimes in horizontal wells ==
 
== Identifying flow regimes in horizontal wells ==
All of these flow regimes in a test can be identified on a diagnostic log-log plot of the pressure change, Δ''p'', and pressure derivative, ''p''′, against the logarithm of time ('''Fig. 7''').  
+
 
 +
All of these flow regimes in a test can be identified on a diagnostic log-log plot of the pressure change, Δ''p'', and pressure derivative, ''p''′, against the logarithm of time ('''Fig. 7''').
  
 
<gallery widths="300px" heights="200px">
 
<gallery widths="300px" heights="200px">
Line 70: Line 79:
 
</gallery>
 
</gallery>
  
A unit-slope line appears during wellbore storage; a horizontal derivative during early radial flow, and then, later, in pseudoradial flow; and a half-slope line in early-linear flow and then in late-linear flow. (These half-slope lines appear on the derivative but not on the pressure-change curves.) This does not imply that all these flow regimes will appear in any given test; in fact, that would be rare. But these are the shapes that identify the flow regimes that may appear in the test being analyzed.  
+
A unit-slope line appears during wellbore storage; a horizontal derivative during early radial flow, and then, later, in pseudoradial flow; and a half-slope line in early-linear flow and then in late-linear flow. (These half-slope lines appear on the derivative but not on the pressure-change curves.) This does not imply that all these flow regimes will appear in any given test; in fact, that would be rare. But these are the shapes that identify the flow regimes that may appear in the test being analyzed.
  
The shapes that may appear in a drawdown test (which is the basis of '''Fig. 7''') may not appear in a buildup test because of the complex superposition of flow regimes. For example, a test would have to be in linear flow both at time (''t''<sub>''p''</sub> + Δ''t'') and at time Δ t to ensure appearance of a derivative with half-slope; this is highly unlikely. The best way to solve the problem is to ensure that a buildup test on a horizontal well is run with a producing time, ''t''<sub>''p''</sub>, much greater than the maximum shut-in time in the test (that is, ''t''<sub>''p''</sub> > 10 Δ''t'' max ).  
+
The shapes that may appear in a drawdown test (which is the basis of '''Fig. 7''') may not appear in a buildup test because of the complex superposition of flow regimes. For example, a test would have to be in linear flow both at time (''t''<sub>''p''</sub> + Δ''t'') and at time Δ t to ensure appearance of a derivative with half-slope; this is highly unlikely. The best way to solve the problem is to ensure that a buildup test on a horizontal well is run with a producing time, ''t''<sub>''p''</sub>, much greater than the maximum shut-in time in the test (that is, ''t''<sub>''p''</sub> > 10 Δ''t'' max ).
  
'''Table 1''' summarizes the working equations for permeability, skin, and start and end of each of the recognized flow regimes. Different investigators have found different equations for start and end of various flow regimes, especially linear flow regimes. This is partly because of a difference of assumptions about flow into the wellbore. Uniform flux or infinite conductivity models are common; neither is rigorously applicable in practice. <ref name="r1" /> In this section, the equations for duration of flow regimes derived by Odeh and Babu<ref name="r2" /> are used. This model assumes uniform flux into the wellbore.  
+
'''Table 1''' summarizes the working equations for permeability, skin, and start and end of each of the recognized flow regimes. Different investigators have found different equations for start and end of various flow regimes, especially linear flow regimes. This is partly because of a difference of assumptions about flow into the wellbore. Uniform flux or infinite conductivity models are common; neither is rigorously applicable in practice. <ref name="r1">Goode, P.A. and Thambynayagam, R.K.M. 1987. Pressure Drawdown and Buildup Analysis of Horizontal Wells in Anisotropic Media. SPE Form Eval 2 (4): 683–697. SPE-14250-PA. http://dx.doi.org/10.2118/14250-PA</ref> In this section, the equations for duration of flow regimes derived by Odeh and Babu<ref name="r2">Odeh, A.S. and Babu, D.K. 1990. Transient Flow Behavior of Horizontal Wells: Pressure Drawdown and Buildup Analysis. SPE Form Eval 5 (1): 7-15. SPE-18802-PA. http://dx.doi.org/10.2118/18802-PA</ref> are used. This model assumes uniform flux into the wellbore.
  
<gallery widths=300px heights=200px>
+
<gallery widths="300px" heights="200px">
 
File:Vol5 Page 0892 Image 0001.png|'''Table 1a'''
 
File:Vol5 Page 0892 Image 0001.png|'''Table 1a'''
  
Line 84: Line 93:
 
</gallery>
 
</gallery>
  
===Early-radial flow===
+
=== Early-radial flow ===
 +
 
 
Early-radial flow is similar to the radial flow period in a vertical well ('''Fig. 2'''). The governing equation for this flow regime is
 
Early-radial flow is similar to the radial flow period in a vertical well ('''Fig. 2'''). The governing equation for this flow regime is
  
[[File:Vol5 page 0809 eq 001.png]]....................(1)
+
[[File:Vol5 page 0809 eq 001.png|RTENOTITLE]]....................(1)
  
Data for this period may be masked by wellbore storage effects, but, when present, they may be analyzed on a semilog plot.  
+
Data for this period may be masked by wellbore storage effects, but, when present, they may be analyzed on a semilog plot.
  
The early-radial flow regime may in theory start at time zero, in absence of wellbore storage effects. The end of the early-radial flow regime may occur when the transient reaches a vertical boundary or when flow comes from beyond the end of the wellbore. The end of the period is the smaller of these two values.  
+
The early-radial flow regime may in theory start at time zero, in absence of wellbore storage effects. The end of the early-radial flow regime may occur when the transient reaches a vertical boundary or when flow comes from beyond the end of the wellbore. The end of the period is the smaller of these two values.
  
'''Eq. 2'''<ref name="r2" /> says that the period must end when the transient reaches the nearest boundary, ''d''<sub>''z''</sub>, from the well. This equation includes the permeability in the vertical direction:
+
'''Eq. 2'''<ref name="r2">Odeh, A.S. and Babu, D.K. 1990. Transient Flow Behavior of Horizontal Wells: Pressure Drawdown and Buildup Analysis. SPE Form Eval 5 (1): 7-15. SPE-18802-PA. http://dx.doi.org/10.2118/18802-PA</ref> says that the period must end when the transient reaches the nearest boundary, ''d''<sub>''z''</sub>, from the well. This equation includes the permeability in the vertical direction:
  
[[File:Vol5 page 0810 eq 001.png]]....................(2)
+
[[File:Vol5 page 0810 eq 001.png|RTENOTITLE]]....................(2)
  
 
The radial flow regime may also end when flow from beyond the end of the wellbore becomes important. '''Eq. 3''' gives the time by
 
The radial flow regime may also end when flow from beyond the end of the wellbore becomes important. '''Eq. 3''' gives the time by
  
[[File:Vol5 page 0810 eq 002.png]]....................(3)
+
[[File:Vol5 page 0810 eq 002.png|RTENOTITLE]]....................(3)
  
''L''<sub>''w''</sub> is the completed length of the well, and k y is the permeability in the direction parallel to the wellbore. The actual end is the lesser of the two times calculated from '''Eqs. 2''' and '''3'''. It is helpful to check the expected duration of the early-radial flow regime after estimating the parameters necessary to make these calculations.  
+
''L''<sub>''w''</sub> is the completed length of the well, and k y is the permeability in the direction parallel to the wellbore. The actual end is the lesser of the two times calculated from '''Eqs. 2''' and '''3'''. It is helpful to check the expected duration of the early-radial flow regime after estimating the parameters necessary to make these calculations.
  
 
'''Eq. 1''' suggests that possible radial flow on the diagnostic plot be identified and then bottomhole flowing pressure be plotted against time during the appropriate time range on semilog coordinates. The slope of the straight line that results is
 
'''Eq. 1''' suggests that possible radial flow on the diagnostic plot be identified and then bottomhole flowing pressure be plotted against time during the appropriate time range on semilog coordinates. The slope of the straight line that results is
  
[[File:Vol5 page 0810 eq 003.png]]....................(4)
+
[[File:Vol5 page 0810 eq 003.png|RTENOTITLE]]....................(4)
  
The group [[File:Vol5 page 0810 inline 001.png]] can be found from the slope, ''m''<sub>erf</sub>:
+
The group [[File:Vol5 page 0810 inline 001.png|RTENOTITLE]] can be found from the slope, ''m''<sub>erf</sub>:
  
[[File:Vol5 page 0810 eq 004.png]]....................(5)
+
[[File:Vol5 page 0810 eq 004.png|RTENOTITLE]]....................(5)
  
Effective completed length of the well must be known to make this calculation. This is not necessarily the same as the perforated or completed length of the well. Some sections of the well may not produce at all.  
+
Effective completed length of the well must be known to make this calculation. This is not necessarily the same as the perforated or completed length of the well. Some sections of the well may not produce at all.
  
 
The equation for calculating the altered permeability skin, ''s''<sub>''d''</sub>, for early-radial flow is
 
The equation for calculating the altered permeability skin, ''s''<sub>''d''</sub>, for early-radial flow is
  
[[File:Vol5 page 0810 eq 005.png]]....................(6)
+
[[File:Vol5 page 0810 eq 005.png|RTENOTITLE]]....................(6)
  
When analyzing a buildup test rather than a constant-rate flow test, plot the HTR or equivalent time on the horizontal axis of the semilog plot, and then plot shut-in or equivalent time on the vertical axis. Note that this plotting is correct only if (''t''<sub>''p''</sub> + Δ''t'') and Δ''t'' appear in this time period simultaneously; that is, radial flow must exist at both time (''t''<sub>''p''</sub> + Δ''t'') and time Δ''t''. This is unlikely because radial-flow regime may exist at time Δ''t'', but a different flow regime is likely at time (''t''<sub>''p''</sub> + Δ''t'').  
+
When analyzing a buildup test rather than a constant-rate flow test, plot the HTR or equivalent time on the horizontal axis of the semilog plot, and then plot shut-in or equivalent time on the vertical axis. Note that this plotting is correct only if (''t''<sub>''p''</sub> + Δ''t'') and Δ''t'' appear in this time period simultaneously; that is, radial flow must exist at both time (''t''<sub>''p''</sub> + Δ''t'') and time Δ''t''. This is unlikely because radial-flow regime may exist at time Δ''t'', but a different flow regime is likely at time (''t''<sub>''p''</sub> + Δ''t'').
  
====Example 1: Well Erf-1====
+
==== Example 1: Well Erf-1 ====
For drawdown test data from Well Erf-1, <ref name="r2" /> the diagnostic plot indicates the data from approximately 0.24 to 24 hours may be in early-radial flow. The following information is available for this well: ''q'' = 800 STB/D, ''μ'' = 1 cp, ''B'' = 1.25 RB/STB, ''r''<sub>''w''</sub> = 0.25 ft, ''ϕ'' = 0.2, ''c''<sub>''t''</sub> = 15×10 –6 psi<sup>–1</sup>, centered in box-shaped drainage area, ''h'' = 200 ft, ''b''<sub>''H''</sub> = 4,000 ft, and ''a''<sub>''H''</sub> = 2,000 ft, ''L''<sub>''w''</sub> =1,000 ft, and, from analysis of data in early linear flow regime, ''k''<sub>''x''</sub> = 200 md. '''Table 2''' shows the pressure change data for 0.24 to 24 hours.
 
  
<gallery widths=300px heights=200px>
+
For drawdown test data from Well Erf-1, <ref name="r2">Odeh, A.S. and Babu, D.K. 1990. Transient Flow Behavior of Horizontal Wells: Pressure Drawdown and Buildup Analysis. SPE Form Eval 5 (1): 7-15. SPE-18802-PA. http://dx.doi.org/10.2118/18802-PA</ref> the diagnostic plot indicates the data from approximately 0.24 to 24 hours may be in early-radial flow. The following information is available for this well: ''q'' = 800 STB/D, ''μ'' = 1 cp, ''B'' = 1.25 RB/STB, ''r''<sub>''w''</sub> = 0.25 ft, ''ϕ'' = 0.2, ''c''<sub>''t''</sub> = 15×10 –6 psi<sup>–1</sup>, centered in box-shaped drainage area, ''h'' = 200 ft, ''b''<sub>''H''</sub> = 4,000 ft, and ''a''<sub>''H''</sub> = 2,000 ft, ''L''<sub>''w''</sub> =1,000 ft, and, from analysis of data in early linear flow regime, ''k''<sub>''x''</sub> = 200 md. '''Table 2''' shows the pressure change data for 0.24 to 24 hours.
 +
 
 +
<gallery widths="300px" heights="200px">
 
File:Vol5 Page 0812 Image 0001.png|'''Table 2'''
 
File:Vol5 Page 0812 Image 0001.png|'''Table 2'''
 
</gallery>
 
</gallery>
Line 128: Line 139:
 
Plot (''p''<sub>''i''</sub> − ''p''<sub>''wf''</sub> ) = Δ''p'' vs. ''t'' on semilog coordinates ('''Fig. 8'''). The plot results in a straight line with a slope of 8 psi/cycle. In '''Fig. 8''', at ''t'' = 2.4 hours, the points begin to deviate from the straight line, as expected from calculations for flow regime duration that follow. The pressure change at 1 hour is 39 psia. Using the slope of 8 and '''Eq. 5''',
 
Plot (''p''<sub>''i''</sub> − ''p''<sub>''wf''</sub> ) = Δ''p'' vs. ''t'' on semilog coordinates ('''Fig. 8'''). The plot results in a straight line with a slope of 8 psi/cycle. In '''Fig. 8''', at ''t'' = 2.4 hours, the points begin to deviate from the straight line, as expected from calculations for flow regime duration that follow. The pressure change at 1 hour is 39 psia. Using the slope of 8 and '''Eq. 5''',
  
[[File:Vol5 page 0811 eq 001.png]]
+
[[File:Vol5 page 0811 eq 001.png|RTENOTITLE]]
  
 
Thus, because ''k''<sub>''x''</sub> = 200 md, ''k''<sub>''z''</sub> = 2 md. Using the value of 39 for Δ''p''<sub>1''hr''</sub> from '''Fig. 8''', skin from '''Eq. 6''' is
 
Thus, because ''k''<sub>''x''</sub> = 200 md, ''k''<sub>''z''</sub> = 2 md. Using the value of 39 for Δ''p''<sub>1''hr''</sub> from '''Fig. 8''', skin from '''Eq. 6''' is
  
[[File:Vol5 page 0811 eq 002.png]]
+
[[File:Vol5 page 0811 eq 002.png|RTENOTITLE]]
  
 
The start of the early-radial flow regime is controlled by wellbore storage, which appears to have vanished at times earlier than 0.24 hours in this example. The end of the early radial flow regime is expected at the lesser of the two values derived from '''Eqs. 2''' and '''3'''. For a centered well, ''d''<sub>''z''</sub> = ''h''/2 = 100 ft, and '''Eq. 2''' gives
 
The start of the early-radial flow regime is controlled by wellbore storage, which appears to have vanished at times earlier than 0.24 hours in this example. The end of the early radial flow regime is expected at the lesser of the two values derived from '''Eqs. 2''' and '''3'''. For a centered well, ''d''<sub>''z''</sub> = ''h''/2 = 100 ft, and '''Eq. 2''' gives
  
[[File:Vol5 page 0811 eq 003.png]]
+
[[File:Vol5 page 0811 eq 003.png|RTENOTITLE]]
  
 
Assuming ''k''<sub>''y''</sub> = ''k''<sub>''x''</sub> = 200 md, '''Eq. 3''' gives
 
Assuming ''k''<sub>''y''</sub> = ''k''<sub>''x''</sub> = 200 md, '''Eq. 3''' gives
  
[[File:Vol5 page 0811 eq 004.png]]
+
[[File:Vol5 page 0811 eq 004.png|RTENOTITLE]]
  
Thus, expect the early-radial flow regime to end at approximately 1.875 hours, which is the smaller value and is consistent with observed test data.  
+
Thus, expect the early-radial flow regime to end at approximately 1.875 hours, which is the smaller value and is consistent with observed test data.
  
 
<gallery widths="300px" heights="200px">
 
<gallery widths="300px" heights="200px">
Line 149: Line 160:
  
 
=== Hemiradial flow ===
 
=== Hemiradial flow ===
The hemiradial flow period ('''Fig. 3''') will occur only when the well is close to one of the vertical boundaries (either the upper or the lower boundaries) and is analogous to a vertical well near a fault. The governing equation is<ref name="r1" />
 
  
[[File:Vol5 page 0812 eq 001.png]]....................(7)
+
The hemiradial flow period ('''Fig. 3''') will occur only when the well is close to one of the vertical boundaries (either the upper or the lower boundaries) and is analogous to a vertical well near a fault. The governing equation is<ref name="r1">Goode, P.A. and Thambynayagam, R.K.M. 1987. Pressure Drawdown and Buildup Analysis of Horizontal Wells in Anisotropic Media. SPE Form Eval 2 (4): 683–697. SPE-14250-PA. http://dx.doi.org/10.2118/14250-PA</ref>
  
A horizontal derivative on the diagnostic plot identifies hemiradial flow. If data appear to fall into this flow regime, a straight line on a semilog plot would provide more confidence that radial flow has been identified. Consistency checks in the analysis coupled with a well survey will be required to distinguish hemiradial flow from early radial flow.  
+
[[File:Vol5 page 0812 eq 001.png|RTENOTITLE]]....................(7)
 +
 
 +
A horizontal derivative on the diagnostic plot identifies hemiradial flow. If data appear to fall into this flow regime, a straight line on a semilog plot would provide more confidence that radial flow has been identified. Consistency checks in the analysis coupled with a well survey will be required to distinguish hemiradial flow from early radial flow.
  
 
The time range in which the analysis for hemiradial flow is valid begins after the closest vertical boundary, ''d''<sub>''z''</sub>, affects the data and before the farthest boundary, ''D''<sub>''z''</sub>, affects them. In the absence of wellbore storage, the start of hemiradial flow is given by
 
The time range in which the analysis for hemiradial flow is valid begins after the closest vertical boundary, ''d''<sub>''z''</sub>, affects the data and before the farthest boundary, ''D''<sub>''z''</sub>, affects them. In the absence of wellbore storage, the start of hemiradial flow is given by
  
[[File:Vol5 page 0812 eq 002.png]]....................(8)
+
[[File:Vol5 page 0812 eq 002.png|RTENOTITLE]]....................(8)
  
Note that the start of the hemiradial flow regime involves the shortest distance to a vertical boundary and the permeability in the vertical direction. However, wellbore storage will most likely determine the actual start of hemiradial flow.  
+
Note that the start of the hemiradial flow regime involves the shortest distance to a vertical boundary and the permeability in the vertical direction. However, wellbore storage will most likely determine the actual start of hemiradial flow.
  
 
The end of hemiradial flow occurs when pressure is affected by the farther vertical boundary or flow from beyond the ends of the wellbore, whichever occurs first. It is the smaller of the times calculated using '''Eqs. 9''' and '''10'''. If the hemiradial flow regime ends when pressure reaches the farthest vertical boundary, it depends on the distance, ''D''<sub>''z''</sub>, and the vertical permeability, ''k''<sub>''z''</sub>:
 
The end of hemiradial flow occurs when pressure is affected by the farther vertical boundary or flow from beyond the ends of the wellbore, whichever occurs first. It is the smaller of the times calculated using '''Eqs. 9''' and '''10'''. If the hemiradial flow regime ends when pressure reaches the farthest vertical boundary, it depends on the distance, ''D''<sub>''z''</sub>, and the vertical permeability, ''k''<sub>''z''</sub>:
  
[[File:Vol5 page 0813 eq 001.png]]....................(9)
+
[[File:Vol5 page 0813 eq 001.png|RTENOTITLE]]....................(9)
  
 
When the appearance of end effects—flow from beyond the ends of the wellbore—causes the end of the hemiradial flow regime to appear, the end of the flow regime occurs when
 
When the appearance of end effects—flow from beyond the ends of the wellbore—causes the end of the hemiradial flow regime to appear, the end of the flow regime occurs when
  
[[File:Vol5 page 0813 eq 002.png]]....................(10)
+
[[File:Vol5 page 0813 eq 002.png|RTENOTITLE]]....................(10)
  
 
The completed length of the well, ''L''<sub>''w''</sub>, and the permeability ''k''<sub>''y''</sub>, parallel to the wellbore appear in this equation. These parameters determine when enough flow has come from beyond the ends of the wellbore to distort the radial flow pattern that appeared earlier.
 
The completed length of the well, ''L''<sub>''w''</sub>, and the permeability ''k''<sub>''y''</sub>, parallel to the wellbore appear in this equation. These parameters determine when enough flow has come from beyond the ends of the wellbore to distort the radial flow pattern that appeared earlier.
  
[[File:Vol5 page 0813 eq 003.png]]....................(11)
+
[[File:Vol5 page 0813 eq 003.png|RTENOTITLE]]....................(11)
  
 
gives the slope of the semilog straight line for semiradial flow, ''m''<sub>hrf</sub>. The multiplier, 325.2, is twice the multiplier for early-radial flow. The equation to estimate the damage skin factor is also similar to that for radial flow but has a multiplier that differs by a factor of two.
 
gives the slope of the semilog straight line for semiradial flow, ''m''<sub>hrf</sub>. The multiplier, 325.2, is twice the multiplier for early-radial flow. The equation to estimate the damage skin factor is also similar to that for radial flow but has a multiplier that differs by a factor of two.
  
[[File:Vol5 page 0813 eq 004.png]]....................(12)
+
[[File:Vol5 page 0813 eq 004.png|RTENOTITLE]]....................(12)
  
The equations relating slope and permeability and the equation for skin are similar in a buildup test to those for a drawdown test. The pressure change in the equation for skin is [''p''<sub>1hr</sub> − ''p''<sub>''wf''</sub> (Δ''t'' = 0)]. Semilog plots of buildup test data from the hemiradial flow regime cannot be analyzed rigorously using data from a Horner plot unless the pressure data at (''t''<sub>''p''</sub> + Δ''t'') and at time Δ''t'' are simultaneously in this flow regime. As a practical matter, the hemiradial flow regime is likely to appear clearly in the buildup test only when the producing time is much greater that the shut-in time.  
+
The equations relating slope and permeability and the equation for skin are similar in a buildup test to those for a drawdown test. The pressure change in the equation for skin is [''p''<sub>1hr</sub> − ''p''<sub>''wf''</sub> (Δ''t'' = 0)]. Semilog plots of buildup test data from the hemiradial flow regime cannot be analyzed rigorously using data from a Horner plot unless the pressure data at (''t''<sub>''p''</sub> + Δ''t'') and at time Δ''t'' are simultaneously in this flow regime. As a practical matter, the hemiradial flow regime is likely to appear clearly in the buildup test only when the producing time is much greater that the shut-in time.
  
 
=== Early linear flow ===
 
=== Early linear flow ===
The governing equation for early-linear flow is<ref name="r1" />
 
  
[[File:Vol5 page 0813 eq 005.png]]....................(13)
+
The governing equation for early-linear flow is<ref name="r1">Goode, P.A. and Thambynayagam, R.K.M. 1987. Pressure Drawdown and Buildup Analysis of Horizontal Wells in Anisotropic Media. SPE Form Eval 2 (4): 683–697. SPE-14250-PA. http://dx.doi.org/10.2118/14250-PA</ref>
  
The "convergence skin," ''s''<sub>''c''</sub>, is discussed later in this section. The start of the early-linear flow regime ('''Fig. 4''') depends on the farthest distance to a vertical boundary, ''D''<sub>''z''</sub>, and the vertical permeability, ''k''<sub>''z''</sub>.<ref name="r2" />
+
[[File:Vol5 page 0813 eq 005.png|RTENOTITLE]]....................(13)
  
[[File:Vol5 page 0814 eq 001.png]]....................(14)
+
The "convergence skin," ''s''<sub>''c''</sub>, is discussed later in this section. The start of the early-linear flow regime ('''Fig. 4''') depends on the farthest distance to a vertical boundary, ''D''<sub>''z''</sub>, and the vertical permeability, ''k''<sub>''z''</sub>.<ref name="r2">Odeh, A.S. and Babu, D.K. 1990. Transient Flow Behavior of Horizontal Wells: Pressure Drawdown and Buildup Analysis. SPE Form Eval 5 (1): 7-15. SPE-18802-PA. http://dx.doi.org/10.2118/18802-PA</ref>
 +
 
 +
[[File:Vol5 page 0814 eq 001.png|RTENOTITLE]]....................(14)
  
 
Not until flow reaches that farthest vertical boundary can a linear flow pattern begin toward the well. This flow period ends when fluids flow from beyond the ends of the wellbore. Thus,
 
Not until flow reaches that farthest vertical boundary can a linear flow pattern begin toward the well. This flow period ends when fluids flow from beyond the ends of the wellbore. Thus,
  
[[File:Vol5 page 0814 eq 002.png]]....................(15)
+
[[File:Vol5 page 0814 eq 002.png|RTENOTITLE]]....................(15)
  
Notice that the end depends on the effective completed length of the well, ''L''<sub>''w''</sub>, and on the permeability in the direction parallel to the well. This is the time in which end effects—flow beyond the ends of the well—begin to significantly distort the linear flow pattern.  
+
Notice that the end depends on the effective completed length of the well, ''L''<sub>''w''</sub>, and on the permeability in the direction parallel to the well. This is the time in which end effects—flow beyond the ends of the well—begin to significantly distort the linear flow pattern.
  
The early-linear flow regime is identified in a drawdown test with a half-slope for the derivative. (Because of the skin effect, the pressure change curve on the diagnostic plot will only approach a half-slope asymptotically.) For data identified as being in this flow regime, plot pressure against the square root of time.  
+
The early-linear flow regime is identified in a drawdown test with a half-slope for the derivative. (Because of the skin effect, the pressure change curve on the diagnostic plot will only approach a half-slope asymptotically.) For data identified as being in this flow regime, plot pressure against the square root of time.
  
 
The slope of the straight line on such a plot, ''m''<sub>elf</sub>, can be used to estimate the square-root of ''k''<sub>''x''</sub>, the horizontal permeability perpendicular to the well:
 
The slope of the straight line on such a plot, ''m''<sub>elf</sub>, can be used to estimate the square-root of ''k''<sub>''x''</sub>, the horizontal permeability perpendicular to the well:
  
[[File:Vol5 page 0814 eq 003.png]]....................(16)
+
[[File:Vol5 page 0814 eq 003.png|RTENOTITLE]]....................(16)
  
 
To calculate the damage skin,
 
To calculate the damage skin,
  
[[File:Vol5 page 0814 eq 004.png]]....................(17)
+
[[File:Vol5 page 0814 eq 004.png|RTENOTITLE]]....................(17)
  
This equation includes a convergence skin, ''s''<sub>''c''</sub>, which is<ref name="r2" />
+
This equation includes a convergence skin, ''s''<sub>''c''</sub>, which is<ref name="r2">Odeh, A.S. and Babu, D.K. 1990. Transient Flow Behavior of Horizontal Wells: Pressure Drawdown and Buildup Analysis. SPE Form Eval 5 (1): 7-15. SPE-18802-PA. http://dx.doi.org/10.2118/18802-PA</ref>
  
[[File:Vol5 page 0814 eq 005.png]]....................(8.203)
+
[[File:Vol5 page 0814 eq 005.png|RTENOTITLE]]....................(8.203)
  
This convergence skin is an additional pressure drop that acts like a skin effect caused by flow moving from throughout the entire formation until it converges down to the small wellbore in the middle of the formation ('''Fig. 9'''). This convergence skin is defined in terms of the ratio of the permeability in the ''x''-direction, which is perpendicular to the wellbore, to the vertical permeability. It also involves the distance to the nearest vertical boundary, ''d''<sub>''z''</sub>, and the net pay thickness, ''h''.  
+
This convergence skin is an additional pressure drop that acts like a skin effect caused by flow moving from throughout the entire formation until it converges down to the small wellbore in the middle of the formation ('''Fig. 9'''). This convergence skin is defined in terms of the ratio of the permeability in the ''x''-direction, which is perpendicular to the wellbore, to the vertical permeability. It also involves the distance to the nearest vertical boundary, ''d''<sub>''z''</sub>, and the net pay thickness, ''h''.
  
 
<gallery widths="300px" heights="200px">
 
<gallery widths="300px" heights="200px">
Line 214: Line 227:
 
</gallery>
 
</gallery>
  
Kuchuk<ref name="r3" /> derived a different equation for convergence skin. As a practical matter, the Odeh-Babu and Kuchuk equations lead to the same result. When there has been a single rate preceding shut-in during early-linear flow, the buildup pressure is plotted against [[File:Vol5 page 0814 inline 001.png]] on a Cartesian plot, which is sometimes called a tandem root plot. The permeability, ''k''<sub>''x''</sub>, is calculated from the slope, ''m''<sub>elf</sub>, of the plot and '''Eq. 16'''. ''k''<sub>''x''</sub> has the same relationship to the slope that existed in a drawdown test. Skin for this flow regime is calculated with '''Eqs. 17''' and '''18'''.  
+
Kuchuk<ref name="r3">Kuchuk, F.J. 1995. Well Testing and Interpretation for Horizontal Wells. J Pet Technol 47 (1): 36–41. SPE-25232-PA. http://dx.doi.org/10.2118/25232-PA</ref> derived a different equation for convergence skin. As a practical matter, the Odeh-Babu and Kuchuk equations lead to the same result. When there has been a single rate preceding shut-in during early-linear flow, the buildup pressure is plotted against [[File:Vol5 page 0814 inline 001.png|RTENOTITLE]] on a Cartesian plot, which is sometimes called a tandem root plot. The permeability, ''k''<sub>''x''</sub>, is calculated from the slope, ''m''<sub>elf</sub>, of the plot and '''Eq. 16'''. ''k''<sub>''x''</sub> has the same relationship to the slope that existed in a drawdown test. Skin for this flow regime is calculated with '''Eqs. 17''' and '''18'''.
 +
 
 +
Plots of buildup data from the early-linear flow regime cannot be analyzed rigorously with a plot of ''p''<sub>''ws''</sub> vs. [[File:Vol5 page 0814 inline 001.png|RTENOTITLE]] (that is, the tandem-root plot) unless data at (''t''<sub>''p''</sub> + Δ''t'') and Δ''t'' are simultaneously in this flow regime—highly unlikely—or unless ''t''<sub>''p''</sub> is much greater than Δ''t'', in which case simply plot ''p''<sub>''ws''</sub> vs. [[File:Vol5 page 0815 inline 001.png|RTENOTITLE]]. Little error results from ignoring the (''t''<sub>''p''</sub> + Δ''t'') term, which is essentially constant.
  
Plots of buildup data from the early-linear flow regime cannot be analyzed rigorously with a plot of ''p''<sub>''ws''</sub> vs. [[File:Vol5 page 0814 inline 001.png]] (that is, the tandem-root plot) unless data at (''t''<sub>''p''</sub> + Δ''t'') and Δ''t'' are simultaneously in this flow regime—highly unlikely—or unless ''t''<sub>''p''</sub> is much greater than Δ''t'', in which case simply plot ''p''<sub>''ws''</sub> vs. [[File:Vol5 page 0815 inline 001.png]]. Little error results from ignoring the (''t''<sub>''p''</sub> + Δ''t'') term, which is essentially constant.
+
==== Example 2: Well Elf-2 ====
  
====Example 2: Well Elf-2====
+
The diagnostic plot for a drawdown test from Well Elf-2<ref name="r2">Odeh, A.S. and Babu, D.K. 1990. Transient Flow Behavior of Horizontal Wells: Pressure Drawdown and Buildup Analysis. SPE Form Eval 5 (1): 7-15. SPE-18802-PA. http://dx.doi.org/10.2118/18802-PA</ref> indicates data in the early-linear flow regime because the derivative has a half-slope. The following data apply to this well: ''q'' = 800 STB/D; ''μ'' = 1 cp; ''B'' = 1.25 RB/STB; ''r''<sub>''w''</sub> = 0.25 ft; ''ϕ'' = 0.2; ''c''<sub>''t''</sub> = 15×10<sup>–6</sup> psi<sup>–1</sup>; centered in box-shaped drainage area 100 ft thick, 4,000 ft long, 4,000 ft wide; ''L''<sub>''w''</sub> = 2,500 ft; and, from early radial-flow regime data, ''k''<sub>''x''</sub>''k''<sub>''z''</sub> = 8,000 md<sup>2</sup>. In addition, '''Table 3''' shows the pressure-change data for this well.
The diagnostic plot for a drawdown test from Well Elf-2<ref name="r2" /> indicates data in the early-linear flow regime because the derivative has a half-slope. The following data apply to this well: ''q'' = 800 STB/D; ''μ'' = 1 cp; ''B'' = 1.25 RB/STB; ''r''<sub>''w''</sub> = 0.25 ft; ''ϕ'' = 0.2; ''c''<sub>''t''</sub> = 15×10<sup>–6</sup> psi<sup>–1</sup>; centered in box-shaped drainage area 100 ft thick, 4,000 ft long, 4,000 ft wide; ''L''<sub>''w''</sub> = 2,500 ft; and, from early radial-flow regime data, ''k''<sub>''x''</sub>''k''<sub>''z''</sub> = 8,000 md<sup>2</sup>. In addition, '''Table 3''' shows the pressure-change data for this well.  
 
  
<gallery widths=300px heights=200px>
+
<gallery widths="300px" heights="200px">
 
File:Vol5 Page 0817 Image 0002.png|'''Table 3'''
 
File:Vol5 Page 0817 Image 0002.png|'''Table 3'''
 
</gallery>
 
</gallery>
Line 227: Line 241:
 
A plot of pressure change vs. the square root of time ('''Fig. 10''') indicates early linear flow. The final point on the straight line is at a time of approximately 24 hours, but there may have been some deviation from the straight line by this time. From the slope of the straight line, ''m''<sub>elf</sub> = 0.934 psi/hr<sup>1/2</sup> and '''Eq. 16''', calculate the permeability in the horizontal plane perpendicular to the wellbore.
 
A plot of pressure change vs. the square root of time ('''Fig. 10''') indicates early linear flow. The final point on the straight line is at a time of approximately 24 hours, but there may have been some deviation from the straight line by this time. From the slope of the straight line, ''m''<sub>elf</sub> = 0.934 psi/hr<sup>1/2</sup> and '''Eq. 16''', calculate the permeability in the horizontal plane perpendicular to the wellbore.
  
[[File:Vol5 page 0815 eq 001.png]]
+
[[File:Vol5 page 0815 eq 001.png|RTENOTITLE]]
  
or ''k''<sub>''x''</sub> = (20.1)<sup>2</sup> = 404 md.  
+
or ''k''<sub>''x''</sub> = (20.1)<sup>2</sup> = 404 md.
  
 
<gallery widths="300px" heights="200px">
 
<gallery widths="300px" heights="200px">
Line 237: Line 251:
 
Analysis of data from the early radial flow regime indicated that ''k''<sub>''x''</sub>''k''<sub>''z''</sub> is 8,000 md<sup>2</sup>; thus, ''k''<sub>''z''</sub> is approximately 19.8 md. To calculate ''s''<sub>''d''</sub>, use '''Eqs. 17''' and '''18''', noting that the value for (''p''<sub>''i''</sub> − ''p''<sub>''wf''</sub> ) = 3.1 at ''t'' = 0.
 
Analysis of data from the early radial flow regime indicated that ''k''<sub>''x''</sub>''k''<sub>''z''</sub> is 8,000 md<sup>2</sup>; thus, ''k''<sub>''z''</sub> is approximately 19.8 md. To calculate ''s''<sub>''d''</sub>, use '''Eqs. 17''' and '''18''', noting that the value for (''p''<sub>''i''</sub> − ''p''<sub>''wf''</sub> ) = 3.1 at ''t'' = 0.
  
[[File:Vol5 page 0815 eq 002.png]]
+
[[File:Vol5 page 0815 eq 002.png|RTENOTITLE]] [[File:Vol5 page 0816 eq 001.png|RTENOTITLE]]
[[File:Vol5 page 0816 eq 001.png]]
 
  
Then, ''s''<sub>''d''</sub> = 4.91 – 4.91 = 0.  
+
Then, ''s''<sub>''d''</sub> = 4.91 – 4.91 = 0.
  
Check the expected time range for the early-linear flow regime.  
+
Check the expected time range for the early-linear flow regime.
  
 
Using ''D''<sub>''z''</sub> = ''h''/2 = 50 ft and ''k''<sub>''z''</sub> = 20 md, calculate the beginning of linear flow with '''Eq. 14''':
 
Using ''D''<sub>''z''</sub> = ''h''/2 = 50 ft and ''k''<sub>''z''</sub> = 20 md, calculate the beginning of linear flow with '''Eq. 14''':
  
[[File:Vol5 page 0816 eq 002.png]]
+
[[File:Vol5 page 0816 eq 002.png|RTENOTITLE]]
  
 
Assuming ''k''<sub>''y''</sub> = ''k''<sub>''x''</sub> at 400 md, use '''Eq. 15''' to find the end of early linear flow.
 
Assuming ''k''<sub>''y''</sub> = ''k''<sub>''x''</sub> at 400 md, use '''Eq. 15''' to find the end of early linear flow.
  
[[File:Vol5 page 0816 eq 003.png]]
+
[[File:Vol5 page 0816 eq 003.png|RTENOTITLE]]
  
These limits are reasonably consistent with the time range analyzed assuming early linear flow.  
+
These limits are reasonably consistent with the time range analyzed assuming early linear flow.
  
 
=== Late pseudoradial ===
 
=== Late pseudoradial ===
  
The governing equation for late-pseudoradial flow is<ref name="r1" />
+
The governing equation for late-pseudoradial flow is<ref name="r1">Goode, P.A. and Thambynayagam, R.K.M. 1987. Pressure Drawdown and Buildup Analysis of Horizontal Wells in Anisotropic Media. SPE Form Eval 2 (4): 683–697. SPE-14250-PA. http://dx.doi.org/10.2118/14250-PA</ref>
  
[[File:Vol5 page 0816 eq 004.png]]....................(19)
+
[[File:Vol5 page 0816 eq 004.png|RTENOTITLE]]....................(19)
  
The late-pseudoradial flow period occurs only if <ref name="r2" />
+
The late-pseudoradial flow period occurs only if <ref name="r2">Odeh, A.S. and Babu, D.K. 1990. Transient Flow Behavior of Horizontal Wells: Pressure Drawdown and Buildup Analysis. SPE Form Eval 5 (1): 7-15. SPE-18802-PA. http://dx.doi.org/10.2118/18802-PA</ref>
  
[[File:Vol5 page 0816 eq 005.png]]....................(20)
+
[[File:Vol5 page 0816 eq 005.png|RTENOTITLE]]....................(20)
  
Here, ''b''<sub>''H''</sub> is the dimension of the reservoir parallel to the wellbore. As long as the completed length of the well is relatively short compared with the length of the drainage area late-pseudoradial flow can occur.  
+
Here, ''b''<sub>''H''</sub> is the dimension of the reservoir parallel to the wellbore. As long as the completed length of the well is relatively short compared with the length of the drainage area late-pseudoradial flow can occur.
  
The start of this flow period occurs when fluid flows from well beyond the ends of the wellbore ('''Fig. 5'''). It is approximated with<ref name="r2" />
+
The start of this flow period occurs when fluid flows from well beyond the ends of the wellbore ('''Fig. 5'''). It is approximated with<ref name="r2">Odeh, A.S. and Babu, D.K. 1990. Transient Flow Behavior of Horizontal Wells: Pressure Drawdown and Buildup Analysis. SPE Form Eval 5 (1): 7-15. SPE-18802-PA. http://dx.doi.org/10.2118/18802-PA</ref>
  
[[File:Vol5 page 0816 eq 006.png]]....................(21)
+
[[File:Vol5 page 0816 eq 006.png|RTENOTITLE]]....................(21)
  
 
This starting time depends on the completed length of the well, ''L''<sub>''w''</sub>, and on the permeability in the direction of the well, ''k''<sub>''y''</sub>. The end of this period, like others in this section, is approximated by the minimum of the results of two calculations. The first,
 
This starting time depends on the completed length of the well, ''L''<sub>''w''</sub>, and on the permeability in the direction of the well, ''k''<sub>''y''</sub>. The end of this period, like others in this section, is approximated by the minimum of the results of two calculations. The first,
  
[[File:Vol5 page 0817 eq 001.png]]....................(22)
+
[[File:Vol5 page 0817 eq 001.png|RTENOTITLE]]....................(22)
  
depends on ''d''<sub>''y''</sub> and the length of the wellbore along with ''k''<sub>''y''</sub>, the permeability in the direction parallel to the wellbore. This is the time at which horizontal boundary effects first appear.  
+
depends on ''d''<sub>''y''</sub> and the length of the wellbore along with ''k''<sub>''y''</sub>, the permeability in the direction parallel to the wellbore. This is the time at which horizontal boundary effects first appear.
  
 
The other equation gives a time at which the radial-flow pattern begins to be distorted depending on the shortest distance, ''d''<sub>''x''</sub>, from the well to a boundary perpendicular to the wellbore and on ''k''<sub>''x''</sub>, the permeability in that direction.
 
The other equation gives a time at which the radial-flow pattern begins to be distorted depending on the shortest distance, ''d''<sub>''x''</sub>, from the well to a boundary perpendicular to the wellbore and on ''k''<sub>''x''</sub>, the permeability in that direction.
  
[[File:Vol5 page 0818 eq 001.png]]....................(23)
+
[[File:Vol5 page 0818 eq 001.png|RTENOTITLE]]....................(23)
  
Whenever boundary effects first appear, whether in a direction that is parallel to the well or perpendicular to the axis of the well, the late-pseudoradial flow period will end.  
+
Whenever boundary effects first appear, whether in a direction that is parallel to the well or perpendicular to the axis of the well, the late-pseudoradial flow period will end.
  
 
The diagnostic plot helps identify the late-pseudoradial flow regime with the characteristic horizontal derivative. For data in the appropriate time range, prepare a semilog plot of pressure against time for a drawdown test. The slope of this plot will be ''m''<sub>prf</sub> and the relationship between that slope and the square root of ''k''<sub>''x''</sub>''k''<sub>''y''</sub>, or the permeabilities in the horizontal plane, is given by
 
The diagnostic plot helps identify the late-pseudoradial flow regime with the characteristic horizontal derivative. For data in the appropriate time range, prepare a semilog plot of pressure against time for a drawdown test. The slope of this plot will be ''m''<sub>prf</sub> and the relationship between that slope and the square root of ''k''<sub>''x''</sub>''k''<sub>''y''</sub>, or the permeabilities in the horizontal plane, is given by
  
[[File:Vol5 page 0818 eq 002.png]]....................(24)
+
[[File:Vol5 page 0818 eq 002.png|RTENOTITLE]]....................(24)
  
 
The skin equation is similar in form to those seen before:
 
The skin equation is similar in form to those seen before:
  
[[File:Vol5 page 0818 eq 003.png]]....................(25)
+
[[File:Vol5 page 0818 eq 003.png|RTENOTITLE]]....................(25)
  
Again, the total skin depends on Δ''p''<sub>1hr</sub>. The convergence skin ('''Eq. 18''') is subtracted from the "total" skin to determine the damage skin.  
+
Again, the total skin depends on Δ''p''<sub>1hr</sub>. The convergence skin ('''Eq. 18''') is subtracted from the "total" skin to determine the damage skin.
  
For a buildup test preceded by production at a single rate, plot pressure against the HTR on a semilog graph. Permeability is calculated from '''Eq. 24''', the same as for a drawdown test. The skin equation is basically the same as for a drawdown test, except that the Δ''p''<sub>1hr</sub> is now ''p''<sub>1hr</sub> − ''p''<sub>''wf''</sub>. To obtain ''p''*, the extrapolated pressure, extrapolate the semilog straight line to a HTR of unity.  
+
For a buildup test preceded by production at a single rate, plot pressure against the HTR on a semilog graph. Permeability is calculated from '''Eq. 24''', the same as for a drawdown test. The skin equation is basically the same as for a drawdown test, except that the Δ''p''<sub>1hr</sub> is now ''p''<sub>1hr</sub> − ''p''<sub>''wf''</sub>. To obtain ''p''*, the extrapolated pressure, extrapolate the semilog straight line to a HTR of unity.
  
Semilog plots of buildup test data from the late-pseudoradial flow regime cannot be analyzed rigorously using a Horner plot unless pressures at (''t''<sub>''p''</sub> + Δ''t'') and at time Δ''t'' are simultaneously in the pseudoradial flow regime, which is highly unlikely. However, little error appears if the producing time before shut-in is much greater than the maximum shut-in time achieved in the buildup test.  
+
Semilog plots of buildup test data from the late-pseudoradial flow regime cannot be analyzed rigorously using a Horner plot unless pressures at (''t''<sub>''p''</sub> + Δ''t'') and at time Δ''t'' are simultaneously in the pseudoradial flow regime, which is highly unlikely. However, little error appears if the producing time before shut-in is much greater than the maximum shut-in time achieved in the buildup test.
  
====Example 3: Well Prf-3====
+
==== Example 3: Well Prf-3 ====
The diagnostic plot suggests that a constant-rate drawdown test from Well Prf-3<ref name="r2" /> includes data in the late-pseudoradial period. The following data are available from the test: ''q'' = 800 STB/D, ''μ'' = 1 cp, ''B'' = 1.25 RB/STB, ''r''<sub>''w''</sub> = 0.25 ft, ''ϕ'' = 0.2, ''c''<sub>''t''</sub> = 15×10<sup>–6</sup> psi<sup>–1</sup>, ''h'' = 150 ft, ''L''<sub>''w''</sub> = 900 ft, ''a''<sub>''H''</sub> = 5,280 ft, ''b''<sub>''H''</sub> = 5,280 ft, well centered in drainage volume, ''k''<sub>''x''</sub> = 100 md (from analysis of early linear flow), ''k''<sub>''x''</sub>''k''<sub>''z''</sub> = 1,000 md<sup>2</sup>, and ''k''<sub>''z''</sub> = 10 md (from analysis of early radial flow). '''Table 4''' gives the pressure change, Δ''p'' = ''p''<font size= "4"><sub>''i''</sub></font> – ''p''<sub>''wf''</sub> vs. time.
 
  
<gallery widths=300px heights=200px>
+
The diagnostic plot suggests that a constant-rate drawdown test from Well Prf-3<ref name="r2">Odeh, A.S. and Babu, D.K. 1990. Transient Flow Behavior of Horizontal Wells: Pressure Drawdown and Buildup Analysis. SPE Form Eval 5 (1): 7-15. SPE-18802-PA. http://dx.doi.org/10.2118/18802-PA</ref> includes data in the late-pseudoradial period. The following data are available from the test: ''q'' = 800 STB/D, ''μ'' = 1 cp, ''B'' = 1.25 RB/STB, ''r''<sub>''w''</sub> = 0.25 ft, ''ϕ'' = 0.2, ''c''<sub>''t''</sub> = 15×10<sup>–6</sup> psi<sup>–1</sup>, ''h'' = 150 ft, ''L''<sub>''w''</sub> = 900 ft, ''a''<sub>''H''</sub> = 5,280 ft, ''b''<sub>''H''</sub> = 5,280 ft, well centered in drainage volume, ''k''<sub>''x''</sub> = 100 md (from analysis of early linear flow), ''k''<sub>''x''</sub>''k''<sub>''z''</sub> = 1,000 md<sup>2</sup>, and ''k''<sub>''z''</sub> = 10 md (from analysis of early radial flow). '''Table 4''' gives the pressure change, Δ''p'' = ''p''<font size="4"><sub>''i''</sub></font> – ''p''<sub>''wf''</sub> vs. time.
 +
 
 +
<gallery widths="300px" heights="200px">
 
File:Vol5 Page 0820 Image 0002.png|'''Table 4'''
 
File:Vol5 Page 0820 Image 0002.png|'''Table 4'''
 
</gallery>
 
</gallery>
Line 305: Line 319:
 
A plot of pressure change vs. the logarithm of time ('''Fig. 11''') confirms pseudoradial flow. A straight line fits all the data from 192 to 432 hours; the slope of the line, ''m''<sub>prf</sub>, is 15.3 psi/cycle, and Δ''p''<sub>1 hr</sub> = 18.94 psi (extrapolated). Then, from '''Eq. 24''',
 
A plot of pressure change vs. the logarithm of time ('''Fig. 11''') confirms pseudoradial flow. A straight line fits all the data from 192 to 432 hours; the slope of the line, ''m''<sub>prf</sub>, is 15.3 psi/cycle, and Δ''p''<sub>1 hr</sub> = 18.94 psi (extrapolated). Then, from '''Eq. 24''',
  
[[File:Vol5 page 0818 eq 004.png]]
+
[[File:Vol5 page 0818 eq 004.png|RTENOTITLE]]
  
Thus, [[File:Vol5 page 0818 eq 005.png]]
+
Thus, [[File:Vol5 page 0818 eq 005.png|RTENOTITLE]]
  
 
From Eqs. 25 and 18 ,
 
From Eqs. 25 and 18 ,
  
[[File:Vol5 page 0819 eq 001.png]]
+
[[File:Vol5 page 0819 eq 001.png|RTENOTITLE]]
  
 
Here,
 
Here,
  
[[File:Vol5 page 0819 eq 002.png]]
+
[[File:Vol5 page 0819 eq 002.png|RTENOTITLE]]
  
 
The pseudoradial flow regime should start at the time given by '''Eq. 21''':
 
The pseudoradial flow regime should start at the time given by '''Eq. 21''':
  
[[File:Vol5 page 0819 eq 003.png]]
+
[[File:Vol5 page 0819 eq 003.png|RTENOTITLE]]
  
 
It should end at a time given by the lesser of values from '''Eqs. 22''' and '''23'''. From '''Eq. 22''',
 
It should end at a time given by the lesser of values from '''Eqs. 22''' and '''23'''. From '''Eq. 22''',
  
[[File:Vol5 page 0819 eq 004.png]]
+
[[File:Vol5 page 0819 eq 004.png|RTENOTITLE]]
  
 
where ''d<sub>y</sub>'' = 1/2(5,200-900) = 2,190ft for this centered well. From '''Eq. 23''',
 
where ''d<sub>y</sub>'' = 1/2(5,200-900) = 2,190ft for this centered well. From '''Eq. 23''',
  
[[File:Vol5 page 0819 eq 005.png]]
+
[[File:Vol5 page 0819 eq 005.png|RTENOTITLE]]
  
The smaller of these two values is 344 hours, which is thus the expected end of pseudoradial flow. The data on the figure that lie on the straight line show the time range from 192 to 432 hours, which is generally consistent with the expected duration of the flow regime.  
+
The smaller of these two values is 344 hours, which is thus the expected end of pseudoradial flow. The data on the figure that lie on the straight line show the time range from 192 to 432 hours, which is generally consistent with the expected duration of the flow regime.
  
 
<gallery widths="300px" heights="200px">
 
<gallery widths="300px" heights="200px">
Line 337: Line 351:
 
=== Late-linear flow ===
 
=== Late-linear flow ===
  
The governing equation for late-linear flow is<ref name="r1" /><ref name="r2" />
+
The governing equation for late-linear flow is<ref name="r1">Goode, P.A. and Thambynayagam, R.K.M. 1987. Pressure Drawdown and Buildup Analysis of Horizontal Wells in Anisotropic Media. SPE Form Eval 2 (4): 683–697. SPE-14250-PA. http://dx.doi.org/10.2118/14250-PA</ref><ref name="r2">Odeh, A.S. and Babu, D.K. 1990. Transient Flow Behavior of Horizontal Wells: Pressure Drawdown and Buildup Analysis. SPE Form Eval 5 (1): 7-15. SPE-18802-PA. http://dx.doi.org/10.2118/18802-PA</ref>
  
[[File:Vol5 page 0819 eq 006.png]]....................(26)
+
[[File:Vol5 page 0819 eq 006.png|RTENOTITLE]]....................(26)
  
The late-linear flow regime starts after the pressure transient has reached the boundaries in the ''z''- and ''y''-directions, and the flow behavior with regard to these directions has become pseudosteady state, as '''Fig. 6''' shows.  
+
The late-linear flow regime starts after the pressure transient has reached the boundaries in the ''z''- and ''y''-directions, and the flow behavior with regard to these directions has become pseudosteady state, as '''Fig. 6''' shows.
  
The start of this time period is the maximum of two equations. <ref name="r2" /> The first depends on the time to reach the boundary, ''D''<sub>''y''</sub>, beyond the end of the horizontal well. It also depends on the permeability, ''k''<sub>''y''</sub>, in the direction parallel to the wellbore.
+
The start of this time period is the maximum of two equations. <ref name="r2">Odeh, A.S. and Babu, D.K. 1990. Transient Flow Behavior of Horizontal Wells: Pressure Drawdown and Buildup Analysis. SPE Form Eval 5 (1): 7-15. SPE-18802-PA. http://dx.doi.org/10.2118/18802-PA</ref> The first depends on the time to reach the boundary, ''D''<sub>''y''</sub>, beyond the end of the horizontal well. It also depends on the permeability, ''k''<sub>''y''</sub>, in the direction parallel to the wellbore.
  
[[File:Vol5 page 0820 eq 001.png]]....................(27)
+
[[File:Vol5 page 0820 eq 001.png|RTENOTITLE]]....................(27)
  
 
Another requirement for the start of the late-linear flow regime is the time to reach the maximum vertical distance, ''D''<sub>''z''</sub>, divided by the vertical permeability:
 
Another requirement for the start of the late-linear flow regime is the time to reach the maximum vertical distance, ''D''<sub>''z''</sub>, divided by the vertical permeability:
  
[[File:Vol5 page 0820 eq 002.png]]....................(28)
+
[[File:Vol5 page 0820 eq 002.png|RTENOTITLE]]....................(28)
  
 
Usually, the start of the late-linear flow regime is dictated by the time to reach the boundaries in the ''y''-direction. The end of this period is given by the equation
 
Usually, the start of the late-linear flow regime is dictated by the time to reach the boundaries in the ''y''-direction. The end of this period is given by the equation
  
[[File:Vol5 page 0821 eq 001.png]]....................(29)
+
[[File:Vol5 page 0821 eq 001.png|RTENOTITLE]]....................(29)
  
The end of the late-linear flow regime depends on reaching the nearest boundary in the direction perpendicular to the wellbore, which is the distance, ''d''<sub>''x''</sub>, away, and on the permeability ''k''<sub>''x''</sub> in that direction.  
+
The end of the late-linear flow regime depends on reaching the nearest boundary in the direction perpendicular to the wellbore, which is the distance, ''d''<sub>''x''</sub>, away, and on the permeability ''k''<sub>''x''</sub> in that direction.
  
 
Identify the late-linear flow regime by a half-slope on the derivative in the diagnostic plot of drawdown test data. (The pressure change may approach a half-slope asymptotically.) For data that appear to be in this flow regime, plot pressure against the square root of time. From the slope ''m''<sub>llf</sub> of the plot, estimate permeability in the ''x''-direction from
 
Identify the late-linear flow regime by a half-slope on the derivative in the diagnostic plot of drawdown test data. (The pressure change may approach a half-slope asymptotically.) For data that appear to be in this flow regime, plot pressure against the square root of time. From the slope ''m''<sub>llf</sub> of the plot, estimate permeability in the ''x''-direction from
  
[[File:Vol5 page 0821 eq 002.png]]....................(30)
+
[[File:Vol5 page 0821 eq 002.png|RTENOTITLE]]....................(30)
  
 
Alternatively, if ''k''<sub>''x''</sub> is known from an early-linear flow regime, estimate ''b''<sub>''H''</sub>, the length of the drainage area, from
 
Alternatively, if ''k''<sub>''x''</sub> is known from an early-linear flow regime, estimate ''b''<sub>''H''</sub>, the length of the drainage area, from
  
[[File:Vol5 page 0821 eq 003.png]]....................(31)
+
[[File:Vol5 page 0821 eq 003.png|RTENOTITLE]]....................(31)
  
 
This late-linear flow regime is the only period that provides the data to calculate the total skin, ''s'', including the partial-penetration skin, ''s''<sub>''p''</sub>, and the convergence skin, ''s''<sub>''c''</sub>. To calculate the damage skin, ''s''<sub>''d''</sub>, use
 
This late-linear flow regime is the only period that provides the data to calculate the total skin, ''s'', including the partial-penetration skin, ''s''<sub>''p''</sub>, and the convergence skin, ''s''<sub>''c''</sub>. To calculate the damage skin, ''s''<sub>''d''</sub>, use
  
[[File:Vol5 page 0821 eq 004.png]]....................(32)
+
[[File:Vol5 page 0821 eq 004.png|RTENOTITLE]]....................(32)
  
The total skin depends on Δ''p''<sub>''t''</sub> = 0. Subtracting the partial penetration skin, ''s''<sub>''p''</sub>, and the convergence skin, ''s''<sub>''c''</sub>, from the total skin yields the damage skin.  
+
The total skin depends on Δ''p''<sub>''t''</sub> = 0. Subtracting the partial penetration skin, ''s''<sub>''p''</sub>, and the convergence skin, ''s''<sub>''c''</sub>, from the total skin yields the damage skin.
  
The partial-penetration skin is a complex function that is calculated with Eqs. A-25 through A-35 in '''Table 1'''. For a buildup test, plot pressure against the HTR. From the slope, ''m''<sub>llf</sub>, calculate ''k''<sub>''x''</sub> with '''Eq. 30''', exactly the same as for drawdown tests. Or, if ''k''<sub>''x''</sub> is known, estimate the length, ''b''<sub>''H''</sub>, of the drainage area with '''Eq. 31'''. Calculate the damage skin, ''s''<sub>''d''</sub>, from a pressure buildup test from '''Eq. 32''', where Δ''p''<sub>''t''</sub> = 0 = (''p''<sub>''t''</sub> = 0)<sub>ext</sub> – ''p''<sub>''wf''(''t'' = 0)</sub>.  
+
The partial-penetration skin is a complex function that is calculated with Eqs. A-25 through A-35 in '''Table 1'''. For a buildup test, plot pressure against the HTR. From the slope, ''m''<sub>llf</sub>, calculate ''k''<sub>''x''</sub> with '''Eq. 30''', exactly the same as for drawdown tests. Or, if ''k''<sub>''x''</sub> is known, estimate the length, ''b''<sub>''H''</sub>, of the drainage area with '''Eq. 31'''. Calculate the damage skin, ''s''<sub>''d''</sub>, from a pressure buildup test from '''Eq. 32''', where Δ''p''<sub>''t''</sub> = 0 = (''p''<sub>''t''</sub> = 0)<sub>ext</sub> – ''p''<sub>''wf''(''t'' = 0)</sub>.
  
Note that the same difficulty arises in using superposition to find plotting functions plots of buildup data from the late-linear flow regime as existed with the previous flow regimes. Pressures at both time (''t''<sub>''p''</sub> + Δ''t'') and time Δ''t'' must be in the late linear flow regime for a tandem-root plot to be valid. However, if ''t''<sub>''p''</sub> >> Δ''t''<sub>max</sub>, there is little error.  
+
Note that the same difficulty arises in using superposition to find plotting functions plots of buildup data from the late-linear flow regime as existed with the previous flow regimes. Pressures at both time (''t''<sub>''p''</sub> + Δ''t'') and time Δ''t'' must be in the late linear flow regime for a tandem-root plot to be valid. However, if ''t''<sub>''p''</sub> >> Δ''t''<sub>max</sub>, there is little error.
  
====Example 4: Well Llf-4====
+
==== Example 4: Well Llf-4 ====
The diagnostic plot for a drawdown test from Well Llf-4<ref name="r2" /> appears to include data in the late-linear flow regime (derivative with half-slope). The following data applies to this well: ''q'' = 800 STB/D, ''μ'' = 1 cp, ''B'' = 1.25 RB/STB, ''r''<sub>''w''</sub> = 0.25 ft, ''ϕ'' = 0.2, ''c''<sub>''t''</sub> = 15 × 10<sup>–6</sup> psi<sup>–1</sup> , ''h'' = 150 ft, ''L''<sub>''w''</sub> = 1,000 ft, ''b''<sub>''H''</sub> = 2,000 ft (well centered), ''a''<sub>''H''</sub> = 6,968 ft (well centered), ''D''<sub>''z''</sub> = 85 ft, ''d''<sub>''z''</sub> = 65 ft, ''k''<sub>''x''</sub>''k''<sub>''z''</sub> = 1,000 md<sup>2</sup> (from analysis of early-radial flow), and ''k''<sub>''x''</sub>''k''<sub>''y''</sub> = 5,000 md<sup>2</sup> (from analysis of pseudoradial flow). '''Table 5''' gives pressure change, Δ''p'' = ''p''<sub>''i''</sub> – ''p''<sub>''wf''</sub> , data vs. time.
 
  
<gallery widths=300px heights=200px>
+
The diagnostic plot for a drawdown test from Well Llf-4<ref name="r2">Odeh, A.S. and Babu, D.K. 1990. Transient Flow Behavior of Horizontal Wells: Pressure Drawdown and Buildup Analysis. SPE Form Eval 5 (1): 7-15. SPE-18802-PA. http://dx.doi.org/10.2118/18802-PA</ref> appears to include data in the late-linear flow regime (derivative with half-slope). The following data applies to this well: ''q'' = 800 STB/D, ''μ'' = 1 cp, ''B'' = 1.25 RB/STB, ''r''<sub>''w''</sub> = 0.25 ft, ''ϕ'' = 0.2, ''c''<sub>''t''</sub> = 15 × 10<sup>–6</sup> psi<sup>–1</sup> , ''h'' = 150 ft, ''L''<sub>''w''</sub> = 1,000 ft, ''b''<sub>''H''</sub> = 2,000 ft (well centered), ''a''<sub>''H''</sub> = 6,968 ft (well centered), ''D''<sub>''z''</sub> = 85 ft, ''d''<sub>''z''</sub> = 65 ft, ''k''<sub>''x''</sub>''k''<sub>''z''</sub> = 1,000 md<sup>2</sup> (from analysis of early-radial flow), and ''k''<sub>''x''</sub>''k''<sub>''y''</sub> = 5,000 md<sup>2</sup> (from analysis of pseudoradial flow). '''Table 5''' gives pressure change, Δ''p'' = ''p''<sub>''i''</sub> – ''p''<sub>''wf''</sub> , data vs. time.
 +
 
 +
<gallery widths="300px" heights="200px">
 
File:Vol5 Page 0825 Image 0002.png|'''Table 5'''
 
File:Vol5 Page 0825 Image 0002.png|'''Table 5'''
 
</gallery>
 
</gallery>
  
'''Fig. 12''' is a plot of pressure change vs. the square root of time. The straight line on this plot for the entire time range (60 to 240 hours) confirms late-linear flow for this range. The slope of the line is 1.56 psi/hr<sup>1/2</sup>, and the intercept is Δ''p''<sub>''t'' = 0</sub> =28.4 psi.  
+
'''Fig. 12''' is a plot of pressure change vs. the square root of time. The straight line on this plot for the entire time range (60 to 240 hours) confirms late-linear flow for this range. The slope of the line is 1.56 psi/hr<sup>1/2</sup>, and the intercept is Δ''p''<sub>''t'' = 0</sub> =28.4 psi.
  
 
<gallery widths="300px" heights="200px">
 
<gallery widths="300px" heights="200px">
Line 390: Line 405:
 
From '''Eq. 30''',
 
From '''Eq. 30''',
  
[[File:Vol5 page 0822 eq 001.png]]
+
[[File:Vol5 page 0822 eq 001.png|RTENOTITLE]]
  
 
Then, ''k''<sub>''x''</sub> = 100 md. Because ''k''<sub>''k''</sub>''k''<sub>''z''</sub> = 100 md<sup>2</sup>, ''k''<sub>''z''</sub> = 10 md. Also, because ''k''<sub>''x''</sub>''k''<sub>''y''</sub> = 5,000 md<sup>2</sup>, ''k''<sub>''y''</sub> =50 md. From '''Eq. 32''',
 
Then, ''k''<sub>''x''</sub> = 100 md. Because ''k''<sub>''k''</sub>''k''<sub>''z''</sub> = 100 md<sup>2</sup>, ''k''<sub>''z''</sub> = 10 md. Also, because ''k''<sub>''x''</sub>''k''<sub>''y''</sub> = 5,000 md<sup>2</sup>, ''k''<sub>''y''</sub> =50 md. From '''Eq. 32''',
  
[[File:Vol5 page 0822 eq 002.png]]
+
[[File:Vol5 page 0822 eq 002.png|RTENOTITLE]]
  
 
From '''Eq. 18''',
 
From '''Eq. 18''',
  
[[File:Vol5 page 0822 eq 003.png]]
+
[[File:Vol5 page 0822 eq 003.png|RTENOTITLE]]
  
 
Calculate the partial penetration skin, s<sub>p</sub> , using the appropriate equation from among Eqs. A-25 through A-35 in '''Table 1'''. First, calculate
 
Calculate the partial penetration skin, s<sub>p</sub> , using the appropriate equation from among Eqs. A-25 through A-35 in '''Table 1'''. First, calculate
  
[[File:Vol5 page 0822 eq 004.png]]
+
[[File:Vol5 page 0822 eq 004.png|RTENOTITLE]]
  
 
to determine whether "Case 1" or "Case 2" applies:
 
to determine whether "Case 1" or "Case 2" applies:
  
[[File:Vol5 page 0822 eq 005.png]]
+
[[File:Vol5 page 0822 eq 005.png|RTENOTITLE]]
  
 
Because
 
Because
  
[[File:Vol5 page 0822 eq 006.png]]
+
[[File:Vol5 page 0822 eq 006.png|RTENOTITLE]]
  
 
this is Case 1 (Eq. A-26). Use Eqs. A-25 through A-31 from the table. From Eq. A-27, ''s''<sub>''p''</sub> = ''p''<sub>''xyz''</sub> + ''p''′<sub>''xy''</sub>. From Eq. A-25,
 
this is Case 1 (Eq. A-26). Use Eqs. A-25 through A-31 from the table. From Eq. A-27, ''s''<sub>''p''</sub> = ''p''<sub>''xyz''</sub> + ''p''′<sub>''xy''</sub>. From Eq. A-25,
  
[[File:Vol5 page 0823 eq 001.png]]
+
[[File:Vol5 page 0823 eq 001.png|RTENOTITLE]]
  
 
From Eq. A-28,
 
From Eq. A-28,
  
[[File:Vol5 page 0823 eq 002.png]]
+
[[File:Vol5 page 0823 eq 002.png|RTENOTITLE]]
  
 
Here, from Eq. A-29,
 
Here, from Eq. A-29,
  
[[File:Vol5 page 0823 eq 003.png]]
+
[[File:Vol5 page 0823 eq 003.png|RTENOTITLE]]
  
 
The well is centered, so ''d''<sub>''y''</sub> = ''D''<sub>''y''</sub> = 500 ft.
 
The well is centered, so ''d''<sub>''y''</sub> = ''D''<sub>''y''</sub> = 500 ft.
  
[[File:Vol5 page 0823 eq 004.png]]
+
[[File:Vol5 page 0823 eq 004.png|RTENOTITLE]]
  
 
From Eq. A-30,
 
From Eq. A-30,
  
[[File:Vol5 page 0823 eq 005.png]]
+
[[File:Vol5 page 0823 eq 005.png|RTENOTITLE]]
  
 
Also,
 
Also,
  
[[File:Vol5 page 0823 eq 006.png]]
+
[[File:Vol5 page 0823 eq 006.png|RTENOTITLE]]
  
 
From Eq. A-31,
 
From Eq. A-31,
  
[[File:Vol5 page 0823 eq 007.png]]
+
[[File:Vol5 page 0823 eq 007.png|RTENOTITLE]]
  
 
and
 
and
  
[[File:Vol5 page 0824 eq 001.png]]
+
[[File:Vol5 page 0824 eq 001.png|RTENOTITLE]]
  
 
Then,
 
Then,
  
[[File:Vol5 page 0824 eq 002.png]]
+
[[File:Vol5 page 0824 eq 002.png|RTENOTITLE]]
  
 
Then,
 
Then,
  
[[File:Vol5 page 0824 eq 003.png]]
+
[[File:Vol5 page 0824 eq 003.png|RTENOTITLE]]
  
 
Now check the expected duration of the late-linear flow regime. The start is the larger of values from '''Eqs. 27''' and '''28'''. From '''Eq. 27''',
 
Now check the expected duration of the late-linear flow regime. The start is the larger of values from '''Eqs. 27''' and '''28'''. From '''Eq. 27''',
  
[[File:Vol5 page 0824 eq 004.png]]
+
[[File:Vol5 page 0824 eq 004.png|RTENOTITLE]]
  
 
From '''Eq. 28''',
 
From '''Eq. 28''',
  
[[File:Vol5 page 0824 eq 005.png]]
+
[[File:Vol5 page 0824 eq 005.png|RTENOTITLE]]
  
 
Thus, the start is expected to be at approximately 162 hours. '''Eq. 29''' gives the end of the flow regime.
 
Thus, the start is expected to be at approximately 162 hours. '''Eq. 29''' gives the end of the flow regime.
  
[[File:Vol5 page 0824 eq 006.png]]
+
[[File:Vol5 page 0824 eq 006.png|RTENOTITLE]]
  
 
The data in this example spanned the time range from 60 to 240 hours. Some of the data that fall on the straight line are, in theory, from times before the start of the late-linear flow regime, but they appear to cause no problem in determining the slope of the line.
 
The data in this example spanned the time range from 60 to 240 hours. Some of the data that fall on the straight line are, in theory, from times before the start of the late-linear flow regime, but they appear to cause no problem in determining the slope of the line.
  
== Field examples<ref name="r4" /> ==
+
== Field examples<ref name="r4">Lichtenberger, G.J. 1994. Data Acquisition and Interpretation of Horizontal Well Pressure-Transient Tests. J Pet Technol 46 (2): 157-162. SPE-25922-PA. http://dx.doi.org/10.2118/25922-PA</ref> ==
  
The following field examples illustrate the procedures used in analyzing horizontal well-test data.  
+
The following field examples illustrate the procedures used in analyzing horizontal well-test data.
 +
 
 +
=== Field example Well A ===
  
===Field example Well A===
 
 
'''Table 6''' summarizes the reservoir and completion properties for Well A. The target for Well A, a horizontal exploration well, was vertical tectonic fracture development in a low-permeability shale. Because of the fractures, the permeability is assumed to be isotropic (''k''<sub>''h''</sub> = ''k''<sub>''z''</sub>) and a result of the fractures. '''Fig. 13''' is a diagnostic plot for Well A and includes a history match using an analytical model.
 
'''Table 6''' summarizes the reservoir and completion properties for Well A. The target for Well A, a horizontal exploration well, was vertical tectonic fracture development in a low-permeability shale. Because of the fractures, the permeability is assumed to be isotropic (''k''<sub>''h''</sub> = ''k''<sub>''z''</sub>) and a result of the fractures. '''Fig. 13''' is a diagnostic plot for Well A and includes a history match using an analytical model.
  
<gallery widths=300px heights=200px>
+
<gallery widths="300px" heights="200px">
 
File:Vol5 Page 0826 Image 0001.png|'''Table 6'''
 
File:Vol5 Page 0826 Image 0001.png|'''Table 6'''
  
Line 479: Line 495:
 
</gallery>
 
</gallery>
  
During the early part of the test there is a unit-slope line representing wellbore storage. Following that wellbore storage, there is a transition to radial flow. The final few data points may be in radial flow. On the Horner plot ('''Fig. 14'''), the last few data points fall on a straight line. From the slope of the straight line, the apparent permeability is 0.011 md and the altered zone skin is 2.9. There is no evidence of boundary effects on this Horner plot. The existence of the semilog straight line is not assured, but the data are at least on the verge of reaching it.  
+
During the early part of the test there is a unit-slope line representing wellbore storage. Following that wellbore storage, there is a transition to radial flow. The final few data points may be in radial flow. On the Horner plot ('''Fig. 14'''), the last few data points fall on a straight line. From the slope of the straight line, the apparent permeability is 0.011 md and the altered zone skin is 2.9. There is no evidence of boundary effects on this Horner plot. The existence of the semilog straight line is not assured, but the data are at least on the verge of reaching it.
  
 
<gallery widths="300px" heights="200px">
 
<gallery widths="300px" heights="200px">
Line 485: Line 501:
 
</gallery>
 
</gallery>
  
Guided by the Horner analysis results, engineers simulated these data with an analytical horizontal well model. The initial match was poor. The match was improved considerably by introducing a no-flow boundary approximately 16 ft above the wellbore, which led to a permeability estimate of 0.027 md and an apparent skin of 11.5. It was concluded that the flow regime observed in the test was hemiradial flow.  
+
Guided by the Horner analysis results, engineers simulated these data with an analytical horizontal well model. The initial match was poor. The match was improved considerably by introducing a no-flow boundary approximately 16 ft above the wellbore, which led to a permeability estimate of 0.027 md and an apparent skin of 11.5. It was concluded that the flow regime observed in the test was hemiradial flow.
  
The final match, shown on the type-curve plot in '''Fig. 13''', is still not a good match at all times, but the author stated that the poor match in the transition region could be the result of phase-redistribution effects in the wellbore. The distance to the no-flow boundary that led to the best match compared favorably with well survey data, which indicated that the well was drilled approximately 20 ft below the upper limit of the productive horizon.  
+
The final match, shown on the type-curve plot in '''Fig. 13''', is still not a good match at all times, but the author stated that the poor match in the transition region could be the result of phase-redistribution effects in the wellbore. The distance to the no-flow boundary that led to the best match compared favorably with well survey data, which indicated that the well was drilled approximately 20 ft below the upper limit of the productive horizon.
  
 
Well B is in a west Texas carbonate formation. It was expected to have isotropic permeability caused by fracturing and dissolution. '''Table 7''' gives the field data for this well.
 
Well B is in a west Texas carbonate formation. It was expected to have isotropic permeability caused by fracturing and dissolution. '''Table 7''' gives the field data for this well.
  
<gallery widths=300px heights=200px>
+
<gallery widths="300px" heights="200px">
 
File:Vol5 Page 0827 Image 0002.png|'''Table 7'''
 
File:Vol5 Page 0827 Image 0002.png|'''Table 7'''
 
</gallery>
 
</gallery>
  
'''Fig. 15''' is the diagnostic plot for this well. After wellbore storage, a short period of radial flow appears to be followed by the onset of linear flow, because p′ approaches a slope of 0.5. In the time region where the derivative is horizontal, a straight line on the Horner plot ('''Fig. 16''') yields ''k'' = 0.14 md. Using these results in the analysis shows that the end of radial flow occurs at ''t''<sub>Erf</sub> = 165 hours.  
+
'''Fig. 15''' is the diagnostic plot for this well. After wellbore storage, a short period of radial flow appears to be followed by the onset of linear flow, because p′ approaches a slope of 0.5. In the time region where the derivative is horizontal, a straight line on the Horner plot ('''Fig. 16''') yields ''k'' = 0.14 md. Using these results in the analysis shows that the end of radial flow occurs at ''t''<sub>Erf</sub> = 165 hours.
  
 
<gallery widths="300px" heights="200px">
 
<gallery widths="300px" heights="200px">
Line 503: Line 519:
 
</gallery>
 
</gallery>
  
 
+
A tandem-root plot ('''Fig. 17''') indicates linear flow and also suggests a distance to the nearest boundary of 29 ft. This is in good agreement with geological observations and helps to verify the assumption of isotropic permeability. The history match with an analytical horizontal well model, shown in '''Fig. 15''', confirms the results of the Horner and tandem-root plots.
A tandem-root plot ('''Fig. 17''') indicates linear flow and also suggests a distance to the nearest boundary of 29 ft. This is in good agreement with geological observations and helps to verify the assumption of isotropic permeability. The history match with an analytical horizontal well model, shown in '''Fig. 15''', confirms the results of the Horner and tandem-root plots.  
 
  
 
<gallery widths="300px" heights="200px">
 
<gallery widths="300px" heights="200px">
Line 510: Line 525:
 
</gallery>
 
</gallery>
  
===Field example Well C===
+
=== Field example Well C ===
 +
 
 
Well C data are from a buildup test of a horizontal well in a high-permeability sandstone where a 54-ft oil column overlies an extensive aquifer estimated to be approximately 180 ft thick. '''Table 8''' shows the available data.
 
Well C data are from a buildup test of a horizontal well in a high-permeability sandstone where a 54-ft oil column overlies an extensive aquifer estimated to be approximately 180 ft thick. '''Table 8''' shows the available data.
  
<gallery widths=300px heights=200px>
+
<gallery widths="300px" heights="200px">
 
File:Vol5 Page 0829 Image 0002.png|'''Table 8'''
 
File:Vol5 Page 0829 Image 0002.png|'''Table 8'''
 
</gallery>
 
</gallery>
  
The diagnostic plot ('''Fig. 18''') shows essentially no wellbore storage and a constant derivative, indicating radial, hemiradial, or elliptical flow at early times. The rapid decline of the derivative at the end of the test is caused by the aquifer underlying the oil column, which is acting like a constant-pressure boundary. A history match with an analytical horizontal well model with one no-flow and one constant-pressure boundary (the lower boundary), yielded ''k''<sub>''h''</sub> = 313 md, ''k''<sub>''z''</sub> = 7.5 md, ''s''<sub>''a''</sub> = 1.5, and ''L''<sub>''w''</sub> = 356 ft. The no-flow boundary was estimated to be approximately 112 ft below the wellbore.  
+
The diagnostic plot ('''Fig. 18''') shows essentially no wellbore storage and a constant derivative, indicating radial, hemiradial, or elliptical flow at early times. The rapid decline of the derivative at the end of the test is caused by the aquifer underlying the oil column, which is acting like a constant-pressure boundary. A history match with an analytical horizontal well model with one no-flow and one constant-pressure boundary (the lower boundary), yielded ''k''<sub>''h''</sub> = 313 md, ''k''<sub>''z''</sub> = 7.5 md, ''s''<sub>''a''</sub> = 1.5, and ''L''<sub>''w''</sub> = 356 ft. The no-flow boundary was estimated to be approximately 112 ft below the wellbore.
  
 
<gallery widths="300px" heights="200px">
 
<gallery widths="300px" heights="200px">
Line 523: Line 539:
 
</gallery>
 
</gallery>
  
The time at which the early radial flow regime ends—the time where the derivative ceases to be flat on the diagnostic plot—is approximately 1.5 hours. For a wellbore with a volume of 130 bbl filled with a fluid of compressibility of 3.5 × 10<sup>–6</sup> psi, the duration of wellbore storage (the unit-slope line) is estimated to be 0.0005 hours. With the gauge sampling rate set at 0.017 hours, the wellbore-storage unit slope simply could not be detected and does not appear at all on this plot. '''Fig. 19''' is the Horner plot for this test. A straight line appears in the same range as the flat derivative on a diagnostic plot. From the slope of the line, the permeability is estimated to be 53 md, close to the regression analysis match value of 48 md.  
+
The time at which the early radial flow regime ends—the time where the derivative ceases to be flat on the diagnostic plot—is approximately 1.5 hours. For a wellbore with a volume of 130 bbl filled with a fluid of compressibility of 3.5 × 10<sup>–6</sup> psi, the duration of wellbore storage (the unit-slope line) is estimated to be 0.0005 hours. With the gauge sampling rate set at 0.017 hours, the wellbore-storage unit slope simply could not be detected and does not appear at all on this plot. '''Fig. 19''' is the Horner plot for this test. A straight line appears in the same range as the flat derivative on a diagnostic plot. From the slope of the line, the permeability is estimated to be 53 md, close to the regression analysis match value of 48 md.
  
 
<gallery widths="300px" heights="200px">
 
<gallery widths="300px" heights="200px">
Line 530: Line 546:
  
 
== Running horizontal well tests ==
 
== Running horizontal well tests ==
The measurements in horizontal wells are usually made above the wellbore with the pressure gauge still in the vertical section.
 
  
The test string may often be too rigid to pass through the wellbore. However, in most cases, conventional hardware can be used for horizontal well tests. With longer horizontal wellbores, wellbore storage is an inherent problem for testing, even for buildup tests with downhole shut-in. As mentioned previously, problems arise in conducting buildup tests with short-duration production periods because superposition is inappropriate; therefore, Horner plots and tandem-root plots, which depend on superposition being applicable, are often inappropriate.  
+
The measurements in horizontal wells are usually made above the wellbore with the pressure gauge still in the vertical section.
  
Another problem in conducting buildup test following a short production period is that significant pressure gradients along the length of the wellbore may cause crossflow within the wellbore during shut-in, so fluid may flow from one region to another in the wellbore. These gradients can be removed and this crossflow eliminated by a longer-duration flow period preceding a buildup test.  
+
The test string may often be too rigid to pass through the wellbore. However, in most cases, conventional hardware can be used for horizontal well tests. With longer horizontal wellbores, wellbore storage is an inherent problem for testing, even for buildup tests with downhole shut-in. As mentioned previously, problems arise in conducting buildup tests with short-duration production periods because superposition is inappropriate; therefore, Horner plots and tandem-root plots, which depend on superposition being applicable, are often inappropriate.
  
===Factors that affect transient responses===
+
Another problem in conducting buildup test following a short production period is that significant pressure gradients along the length of the wellbore may cause crossflow within the wellbore during shut-in, so fluid may flow from one region to another in the wellbore. These gradients can be removed and this crossflow eliminated by a longer-duration flow period preceding a buildup test.
A number of factors may affect the transient response of a horizontal well test: horizontal permeability (normal and parallel to well trajectory), vertical permeability, drilling damage, completion damage, producing interval that may be effectively much less than drilled length, and variations in standoff along length of well.  
 
  
In summary, seven or more factors may affect interpretation for horizontal wells in homogeneous reservoirs before the effects from boundaries. The problem is complex, so test results are frequently inconclusive. Furthermore, wellbore storage inhibits determination of properties associated with early-time transient behavior such as vertical permeability and damage from drilling and completion. Middle- and late-time behavior may require several hours, days, or months to appear in transient data.
+
=== Factors that affect transient responses ===
  
Some practical steps will help ensure interpretable test data. First, it is helpful to run tests in the pilot hole before kicking off to drill the horizontal borehole section. From a test in the vertical section, it is possible to get usable estimates of horizontal and vertical permeabilities using modern wireline test tools. Second, a good directional drilling survey can frequently provide an adequate estimate of standoff. A production-log flow survey conducted with coiled tubing can determine what part of the wellbore is actually producing and, therefore, help provide an estimate of effective productive length.  
+
A number of factors may affect the transient response of a horizontal well test: horizontal permeability (normal and parallel to well trajectory), vertical permeability, drilling damage, completion damage, producing interval that may be effectively much less than drilled length, and variations in standoff along length of well.
 +
 
 +
In summary, seven or more factors may affect interpretation for horizontal wells in homogeneous reservoirs before the effects from boundaries. The problem is complex, so test results are frequently inconclusive. Furthermore, wellbore storage inhibits determination of properties associated with early-time transient behavior such as vertical permeability and damage from drilling and completion. Middle- and late-time behavior may require several hours, days, or months to appear in transient data.
 +
 
 +
Some practical steps will help ensure interpretable test data. First, it is helpful to run tests in the pilot hole before kicking off to drill the horizontal borehole section. From a test in the vertical section, it is possible to get usable estimates of horizontal and vertical permeabilities using modern wireline test tools. Second, a good directional drilling survey can frequently provide an adequate estimate of standoff. A production-log flow survey conducted with coiled tubing can determine what part of the wellbore is actually producing and, therefore, help provide an estimate of effective productive length.
  
 
Wells in developed reservoirs should be flowed long enough to bring pressures along the wellbore to equilibrium and thus minimize crossflow. For high-rate wells, continuous borehole pressure and flow-rate measurements acquired during production can be used to interpret the pressure-drawdown transient response. If the downhole rates are not measured, the buildup test should be conducted with downhole shut-in to minimize wellbore storage distortion of test data.
 
Wells in developed reservoirs should be flowed long enough to bring pressures along the wellbore to equilibrium and thus minimize crossflow. For high-rate wells, continuous borehole pressure and flow-rate measurements acquired during production can be used to interpret the pressure-drawdown transient response. If the downhole rates are not measured, the buildup test should be conducted with downhole shut-in to minimize wellbore storage distortion of test data.
  
 
== Nomenclature ==
 
== Nomenclature ==
 +
 
{|
 
{|
|''a''
 
|=
 
|[[File:Vol5 page 0879 inline 001.png]], stabilized deliverability coefficient, psia<sup>2</sup>-cp/MMscf-D
 
 
|-
 
|-
|''a''  
+
| ''a''
|=  
+
| =
|total length of reservoir perpendicular to wellbore, ft  
+
| [[File:Vol5 page 0879 inline 001.png|RTENOTITLE]], stabilized deliverability coefficient, psia<sup>2</sup>-cp/MMscf-D
 +
|-
 +
| ''a''
 +
| =
 +
| total length of reservoir perpendicular to wellbore, ft
 
|-
 
|-
|''a''<sub>''h''</sub>  
+
| ''a''<sub>''h''</sub>
|=  
+
| =
|length of reservoir perpendicular to horizontal well, ft  
+
| length of reservoir perpendicular to horizontal well, ft
 
|-
 
|-
|''a''<sub>''f''</sub>  
+
| ''a''<sub>''f''</sub>
|=  
+
| =
|[[File:Vol5 page 0879 inline 002.png]], depth of investigation along major axis in fractured well, ft  
+
| [[File:Vol5 page 0879 inline 002.png|RTENOTITLE]], depth of investigation along major axis in fractured well, ft
 
|-
 
|-
|''a''<sub>''t''</sub>  
+
| ''a''<sub>''t''</sub>
|=  
+
| =
|[[File:Vol5 page 0879 inline 002.png]], transient deliverability coefficient, psia<sup>2</sup>-cp/MMscf-D  
+
| [[File:Vol5 page 0879 inline 002.png|RTENOTITLE]], transient deliverability coefficient, psia<sup>2</sup>-cp/MMscf-D
 
|-
 
|-
|''a''<sub>''H''</sub>  
+
| ''a''<sub>''H''</sub>
|=  
+
| =
|total width of reservoir perpendicular to the wellbore, ft  
+
| total width of reservoir perpendicular to the wellbore, ft
 
|-
 
|-
|''a''<sub>''H''</sub>′  
+
| ''a''<sub>''H''</sub>′
|=  
+
| =
|modified total width of reservoir perpendicular to the wellbore, ft  
+
| modified total width of reservoir perpendicular to the wellbore, ft
 
|-
 
|-
|''A''  
+
| ''A''
|=  
+
| =
|drainage area, sq ft  
+
| drainage area, sq ft
 
|-
 
|-
|''A''  
+
| ''A''
|=  
+
| =
|''πa''<sub>''f''</sub>''b''<sub>''f''</sub> , area of investigation in fractured well, ft<sup>2</sup>  
+
| ''πa''<sub>''f''</sub>''b''<sub>''f''</sub> , area of investigation in fractured well, ft<sup>2</sup>
 
|-
 
|-
|''A''<sub>''f''</sub>  
+
| ''A''<sub>''f''</sub>
|=  
+
| =
|cross-sectional area perpendicular to flow, sq ft  
+
| cross-sectional area perpendicular to flow, sq ft
 
|-
 
|-
|''A''<sub>''wb''</sub>  
+
| ''A''<sub>''wb''</sub>
|=  
+
| =
|wellbore area, sq ft  
+
| wellbore area, sq ft
 
|-
 
|-
|''b''  
+
| ''b''
|=  
+
| =
|[[File:Vol5 page 0880 inline 001.png]] (gas flow equation)  
+
| [[File:Vol5 page 0880 inline 001.png|RTENOTITLE]] (gas flow equation)
 
|-
 
|-
|''b''<sub>''f''</sub>  
+
| ''b''<sub>''f''</sub>
|=  
+
| =
|[[File:Vol5 page 0880 inline 002.png]], depth of investigation of along minor axis in fractured well, ft  
+
| [[File:Vol5 page 0880 inline 002.png|RTENOTITLE]], depth of investigation of along minor axis in fractured well, ft
 
|-
 
|-
|''b''<sub>''B''</sub>  
+
| ''b''<sub>''B''</sub>
|=  
+
| =
|intercept of Cartesian plot of bilinear flow data, psi  
+
| intercept of Cartesian plot of bilinear flow data, psi
 
|-
 
|-
|''b''<sub>''H''</sub>  
+
| ''b''<sub>''H''</sub>
|=  
+
| =
|length in direction parallel to wellbore, ft  
+
| length in direction parallel to wellbore, ft
 
|-
 
|-
|''b''<sub>''H''</sub>′  
+
| ''b''<sub>''H''</sub>′
|=  
+
| =
|modified length in direction parallel to wellbore, ft  
+
| modified length in direction parallel to wellbore, ft
 
|-
 
|-
|''b''<sub>''L''</sub>  
+
| ''b''<sub>''L''</sub>
|=  
+
| =
|intercept of Cartesian plot of linear flow data, psi  
+
| intercept of Cartesian plot of linear flow data, psi
 
|-
 
|-
|''b''<sub>''V''</sub>  
+
| ''b''<sub>''V''</sub>
|=  
+
| =
|intercept of Cartesian plot of data during volumetric behavior, psi  
+
| intercept of Cartesian plot of data during volumetric behavior, psi
 
|-
 
|-
|''B''  
+
| ''B''
|=  
+
| =
|formation volume factor, res vol/surface vol  
+
| formation volume factor, res vol/surface vol
 
|-
 
|-
|''B''<sub>''g''</sub>  
+
| ''B''<sub>''g''</sub>
|=  
+
| =
|gas formation volume factor, RB/STB  
+
| gas formation volume factor, RB/STB
 
|-
 
|-
|''B''<sub>''gi''</sub>  
+
| ''B''<sub>''gi''</sub>
|=  
+
| =
|gas formation volume factor evaluated at ''p''<sub>''i''</sub>, RB/Mscf  
+
| gas formation volume factor evaluated at ''p''<sub>''i''</sub>, RB/Mscf
 
|-
 
|-
|''B''<sub>''o''</sub>  
+
| ''B''<sub>''o''</sub>
|=  
+
| =
|oil formation volume factor, RB/Mscf  
+
| oil formation volume factor, RB/Mscf
 
|-
 
|-
|''B''<sub>''w''</sub>  
+
| ''B''<sub>''w''</sub>
|=  
+
| =
|water formation volume factor, RB/STB  
+
| water formation volume factor, RB/STB
 
|-
 
|-
|[[File:Vol5 page 0880 inline 003.png]]  
+
| [[File:Vol5 page 0880 inline 003.png|RTENOTITLE]]
|=  
+
| =
|gas formation volume factor evaluated at average drainage area pressure, RB/Mscf  
+
| gas formation volume factor evaluated at average drainage area pressure, RB/Mscf
 
|-
 
|-
|''B''<sub>ND</sub>  
+
| ''B''<sub>ND</sub>
|=  
+
| =
|1,422 ''μ'' ''z'' TD/''kh'', non-Darcy flow coefficient  
+
| 1,422 ''μ'' ''z'' TD/''kh'', non-Darcy flow coefficient
 
|-
 
|-
|''c''  
+
| ''c''
|=  
+
| =
|compressibility, psi<sup>–1</sup>  
+
| compressibility, psi<sup>–1</sup>
 
|-
 
|-
|''c''<sub>''f''</sub>  
+
| ''c''<sub>''f''</sub>
|=  
+
| =
|formation compressibility, psi<sup>–1</sup>  
+
| formation compressibility, psi<sup>–1</sup>
 
|-
 
|-
|''c''<sub>''g''</sub>  
+
| ''c''<sub>''g''</sub>
|=  
+
| =
|gas compressibility, psi<sup>–1</sup>  
+
| gas compressibility, psi<sup>–1</sup>
 
|-
 
|-
|''c''<sub>''o''</sub>  
+
| ''c''<sub>''o''</sub>
|=  
+
| =
|oil compressibility, psi<sup>–1</sup>  
+
| oil compressibility, psi<sup>–1</sup>
 
|-
 
|-
|''c''<sub>''t''</sub>  
+
| ''c''<sub>''t''</sub>
|=  
+
| =
|''S''<sub>''o''</sub>''c''<sub>''o''</sub> + ''S''<sub>''w''</sub>''c''<sub>''w''</sub> + ''S''<sub>''g''</sub>''c''<sub>''g''</sub> + ''c''<sub>''f''</sub> = total compressibility, psi<sup>–1</sup>  
+
| ''S''<sub>''o''</sub>''c''<sub>''o''</sub> + ''S''<sub>''w''</sub>''c''<sub>''w''</sub> + ''S''<sub>''g''</sub>''c''<sub>''g''</sub> + ''c''<sub>''f''</sub> = total compressibility, psi<sup>–1</sup>
 
|-
 
|-
|''c''<sub>''w''</sub>  
+
| ''c''<sub>''w''</sub>
|=  
+
| =
|water compressibility, psi<sup>–1</sup>  
+
| water compressibility, psi<sup>–1</sup>
 
|-
 
|-
|[[File:Vol5 page 0789 inline 002.png]]  
+
| [[File:Vol5 page 0789 inline 002.png|RTENOTITLE]]
|=  
+
| =
|total compressibility evaluated at average drainage area pressure, psi<sup>–1</sup>  
+
| total compressibility evaluated at average drainage area pressure, psi<sup>–1</sup>
 
|-
 
|-
|''c''<sub>''tf''</sub>  
+
| ''c''<sub>''tf''</sub>
|=  
+
| =
|total compressibility of pore space and fluids in fracture porosity, psi<sup>–1</sup>  
+
| total compressibility of pore space and fluids in fracture porosity, psi<sup>–1</sup>
 
|-
 
|-
|''c''<sub>''tm''</sub>  
+
| ''c''<sub>''tm''</sub>
|=  
+
| =
|total compressibility of pore space and fluids in matrix porosity, psi<sup>–1</sup>  
+
| total compressibility of pore space and fluids in matrix porosity, psi<sup>–1</sup>
 
|-
 
|-
|''c''<sub>''wb''</sub>  
+
| ''c''<sub>''wb''</sub>
|=  
+
| =
|compressibility of fluid in wellbore, psi<sup>–1</sup>  
+
| compressibility of fluid in wellbore, psi<sup>–1</sup>
 
|-
 
|-
|''C''  
+
| ''C''
|=  
+
| =
|performance coefficient in gas-well deliverability equation, or wellbore storage coefficient, bbl/psi  
+
| performance coefficient in gas-well deliverability equation, or wellbore storage coefficient, bbl/psi
 
|-
 
|-
|''C''<sub>''A''</sub>  
+
| ''C''<sub>''A''</sub>
|=  
+
| =
|shape factor or constant  
+
| shape factor or constant
 
|-
 
|-
|''C''<sub>''D''</sub>  
+
| ''C''<sub>''D''</sub>
|=  
+
| =
|0.8936 ''C''/''ϕc''<sub>''t''</sub>''hr''<sub>''w''</sub><sup>2</sup> , dimensionless wellbore storage coefficient  
+
| 0.8936 ''C''/''ϕc''<sub>''t''</sub>''hr''<sub>''w''</sub><sup>2</sup> , dimensionless wellbore storage coefficient
 
|-
 
|-
|(''C''<sub>''D''</sub>''e''<sup>2s</sup>)<sub>''f''</sub>  
+
| (''C''<sub>''D''</sub>''e''<sup>2s</sup>)<sub>''f''</sub>
|=  
+
| =
|type-curve parameter value for the formation  
+
| type-curve parameter value for the formation
 
|-
 
|-
|(''C''<sub>''D''</sub>''e''<sup>2s</sup>)<sub>''f''</sub> <sub>+''m''</sub>  
+
| (''C''<sub>''D''</sub>''e''<sup>2s</sup>)<sub>''f''</sub> <sub>+''m''</sub>
|=  
+
| =
|type-curve parameter value for the formation plus the matrix  
+
| type-curve parameter value for the formation plus the matrix
 
|-
 
|-
|''C''<sub>''LfD''</sub>  
+
| ''C''<sub>''LfD''</sub>
|=  
+
| =
|0.8936 [[File:Vol5 page 0880 inline 004.png]], dimensionless wellbore storage coefficient in fractured well  
+
| 0.8936 [[File:Vol5 page 0880 inline 004.png|RTENOTITLE]], dimensionless wellbore storage coefficient in fractured well
 
|-
 
|-
|''C''<sub>''r''</sub>  
+
| ''C''<sub>''r''</sub>
|=  
+
| =
|''w''<sub>''f''</sub>''k''<sub>''f''</sub>/''πkL''<sub>''f''</sub>, fracture conductivity, dimensionless  
+
| ''w''<sub>''f''</sub>''k''<sub>''f''</sub>/''πkL''<sub>''f''</sub>, fracture conductivity, dimensionless
 
|-
 
|-
|''d''<sub>''x''</sub>  
+
| ''d''<sub>''x''</sub>
|=  
+
| =
|shortest distance between horizontal well and ''x'' boundary, ft  
+
| shortest distance between horizontal well and ''x'' boundary, ft
 
|-
 
|-
|''d''<sub>''y''</sub>  
+
| ''d''<sub>''y''</sub>
|=  
+
| =
|shortest distance between tip of horizontal well and ''y'' boundary, ft  
+
| shortest distance between tip of horizontal well and ''y'' boundary, ft
 
|-
 
|-
|''d''<sub>''z''</sub>  
+
| ''d''<sub>''z''</sub>
|=  
+
| =
|shortest distance between horizontal well and ''z'' boundary, ft  
+
| shortest distance between horizontal well and ''z'' boundary, ft
 
|-
 
|-
|''D''<sub>''x''</sub>  
+
| ''D''<sub>''x''</sub>
|=  
+
| =
|longest distance between horizontal well and x boundary, ft  
+
| longest distance between horizontal well and x boundary, ft
 
|-
 
|-
|''D''<sub>''y''</sub>  
+
| ''D''<sub>''y''</sub>
|=  
+
| =
|longest distance between tip of horizontal well and y boundary, ft  
+
| longest distance between tip of horizontal well and y boundary, ft
 
|-
 
|-
|''D''<sub>''z''</sub>  
+
| ''D''<sub>''z''</sub>
|=  
+
| =
|longest distance between horizontal well and z boundary, ft  
+
| longest distance between horizontal well and z boundary, ft
 
|-
 
|-
|''D''  
+
| ''D''
|=  
+
| =
|non-Darcy flow constant, D/Mscf  
+
| non-Darcy flow constant, D/Mscf
 
|-
 
|-
|''e''<sup>–''bt''</sup>  
+
| ''e''<sup>–''bt''</sup>
|=  
+
| =
|exponential decline with a constant ''b'' and elapsed time, ''t''  
+
| exponential decline with a constant ''b'' and elapsed time, ''t''
 
|-
 
|-
|''E''<sub>''f''</sub>  
+
| ''E''<sub>''f''</sub>
|=  
+
| =
|flow efficiency, dimensionless  
+
| flow efficiency, dimensionless
 
|-
 
|-
|''Ei''(–''x'')  
+
| ''Ei''(–''x'')
|=  
+
| =
|[[File:Vol5 page 0881 inline 001.png]], the exponential integral  
+
| [[File:Vol5 page 0881 inline 001.png|RTENOTITLE]], the exponential integral
 
|-
 
|-
|''F''(''u'')  
+
| ''F''(''u'')
|=  
+
| =
|function used in horizontal well analysis  
+
| function used in horizontal well analysis
 
|-
 
|-
|''F''<sub>''CD''</sub>  
+
| ''F''<sub>''CD''</sub>
|=  
+
| =
|''w''<sub>''f''</sub>''k''<sub>''f''</sub>/''kL''<sub>''f''</sub>, fracture conductivity, dimensionless  
+
| ''w''<sub>''f''</sub>''k''<sub>''f''</sub>/''kL''<sub>''f''</sub>, fracture conductivity, dimensionless
 
|-
 
|-
|''g''  
+
| ''g''
|=  
+
| =
|acceleration due to gravity, ft/sec<sup>2</sup>  
+
| acceleration due to gravity, ft/sec<sup>2</sup>
 
|-
 
|-
|''g''<sub>''c''</sub>  
+
| ''g''<sub>''c''</sub>
|=  
+
| =
|gravitational units conversion factor, 32.17 (lbm/ft)(lbf-s<sup>2</sup>)  
+
| gravitational units conversion factor, 32.17 (lbm/ft)(lbf-s<sup>2</sup>)
 
|-
 
|-
|''h''  
+
| ''h''
|=  
+
| =
|net formation thickness, ft  
+
| net formation thickness, ft
 
|-
 
|-
|''h''<sub>''D''</sub>  
+
| ''h''<sub>''D''</sub>
|=  
+
| =
|(''h''/''r''<sub>''w''</sub>)(''k''<sub>''h''</sub>/''k''<sub>''v''</sub>)<sup>1/2</sup>, dimensionless  
+
| (''h''/''r''<sub>''w''</sub>)(''k''<sub>''h''</sub>/''k''<sub>''v''</sub>)<sup>1/2</sup>, dimensionless
 
|-
 
|-
|''h''<sub>''f''</sub>  
+
| ''h''<sub>''f''</sub>
|=  
+
| =
|fracture height, ft  
+
| fracture height, ft
 
|-
 
|-
|''h''<sub>''m''</sub>  
+
| ''h''<sub>''m''</sub>
|=  
+
| =
|thickness of matrix, ft  
+
| thickness of matrix, ft
 
|-
 
|-
|''h''<sub>''p''</sub>  
+
| ''h''<sub>''p''</sub>
|=  
+
| =
|perforated interval thickness, ft  
+
| perforated interval thickness, ft
 
|-
 
|-
|''h''<sub>''pD''</sub>  
+
| ''h''<sub>''pD''</sub>
|=  
+
| =
|''h''<sub>''p''</sub>/''h''<sub>''t''</sub>  
+
| ''h''<sub>''p''</sub>/''h''<sub>''t''</sub>
 
|-
 
|-
|''h''<sub>''t''</sub>  
+
| ''h''<sub>''t''</sub>
|=  
+
| =
|total formation thickness, ft  
+
| total formation thickness, ft
 
|-
 
|-
|''h''<sub>''1''</sub>  
+
| ''h''<sub>''1''</sub>
|=  
+
| =
|distance from top of formation to top of perforations, ft  
+
| distance from top of formation to top of perforations, ft
 
|-
 
|-
|''h''<sub>''1D''</sub>  
+
| ''h''<sub>''1D''</sub>
|=  
+
| =
|''h''<sub>''1''</sub>/''h''<sub>''t''</sub>  
+
| ''h''<sub>''1''</sub>/''h''<sub>''t''</sub>
 
|-
 
|-
|HTR<sub>avg</sub>  
+
| HTR<sub>avg</sub>
|=  
+
| =
|HTR at average drainage area pressure  
+
| HTR at average drainage area pressure
 
|-
 
|-
|''J''  
+
| ''J''
|=  
+
| =
|productivity index, STB/D, psi  
+
| productivity index, STB/D, psi
 
|-
 
|-
|''J''<sub>actual</sub>  
+
| ''J''<sub>actual</sub>
|=  
+
| =
|actual well productivity index, STB/D-psi  
+
| actual well productivity index, STB/D-psi
 
|-
 
|-
|''J''<sub>ideal</sub>  
+
| ''J''<sub>ideal</sub>
|=  
+
| =
|ideal productivity index (''s'' = 0), STB/D-psi  
+
| ideal productivity index (''s'' = 0), STB/D-psi
 
|-
 
|-
|''k''  
+
| ''k''
|=  
+
| =
|matrix permeability, md  
+
| matrix permeability, md
 
|-
 
|-
|[[File:Vol5 page 0800 inline 001.png]]  
+
| [[File:Vol5 page 0800 inline 001.png|RTENOTITLE]]
|=  
+
| =
|average permeability, md  
+
| average permeability, md
 
|-
 
|-
|''k''<sub>''f''</sub>  
+
| ''k''<sub>''f''</sub>
|=  
+
| =
|permeability of the proppant in the fracture, md  
+
| permeability of the proppant in the fracture, md
 
|-
 
|-
|''k''<sub>''fs''</sub>  
+
| ''k''<sub>''fs''</sub>
|=  
+
| =
|permeability near the wellbore, md  
+
| permeability near the wellbore, md
 
|-
 
|-
|''k''<sub>''g''</sub>  
+
| ''k''<sub>''g''</sub>
|=  
+
| =
|permeability to gas, md  
+
| permeability to gas, md
 
|-
 
|-
|''k''<sub>''gp''</sub>  
+
| ''k''<sub>''gp''</sub>
|=  
+
| =
|permeability of the gravel in the gravel pack, md  
+
| permeability of the gravel in the gravel pack, md
 
|-
 
|-
|''k''<sub>''h''</sub>  
+
| ''k''<sub>''h''</sub>
|=  
+
| =
|horizontal permeability, md  
+
| horizontal permeability, md
 
|-
 
|-
|''k''<sub>''m''</sub>  
+
| ''k''<sub>''m''</sub>
|=  
+
| =
|matrix permeability, md  
+
| matrix permeability, md
 
|-
 
|-
|''k''<sub>''o''</sub>  
+
| ''k''<sub>''o''</sub>
|=  
+
| =
|permeability to oil, md  
+
| permeability to oil, md
 
|-
 
|-
|''k''<sub>''r''</sub>  
+
| ''k''<sub>''r''</sub>
|=  
+
| =
|permeability in horizontal radial direction, md  
+
| permeability in horizontal radial direction, md
 
|-
 
|-
|''k''<sub>''s''</sub>  
+
| ''k''<sub>''s''</sub>
|=  
+
| =
|permeability of altered zone, md  
+
| permeability of altered zone, md
 
|-
 
|-
|''k''<sub>''w''</sub>  
+
| ''k''<sub>''w''</sub>
|=  
+
| =
|permeability to water, md  
+
| permeability to water, md
 
|-
 
|-
|''k''<sub>''x''</sub>  
+
| ''k''<sub>''x''</sub>
|=  
+
| =
|permeability in ''x''-direction, md  
+
| permeability in ''x''-direction, md
 
|-
 
|-
|''k''<sub>''y''</sub>  
+
| ''k''<sub>''y''</sub>
|=  
+
| =
|permeability in ''y''-direction, md  
+
| permeability in ''y''-direction, md
 
|-
 
|-
|''k''<sub>''z''</sub>  
+
| ''k''<sub>''z''</sub>
|=  
+
| =
|permeability in ''z''-direction, md  
+
| permeability in ''z''-direction, md
 
|-
 
|-
|''L''  
+
| ''L''
|=  
+
| =
|distance from well to no-flow boundary, ft  
+
| distance from well to no-flow boundary, ft
 
|-
 
|-
|''L''<sub>''d''</sub>  
+
| ''L''<sub>''d''</sub>
|=  
+
| =
|drilled length of horizontal well, ft  
+
| drilled length of horizontal well, ft
 
|-
 
|-
|''L''<sub>''f''</sub>  
+
| ''L''<sub>''f''</sub>
|=  
+
| =
|fracture half length, ft  
+
| fracture half length, ft
 
|-
 
|-
|''L''<sub>''g''</sub>  
+
| ''L''<sub>''g''</sub>
|=  
+
| =
|length of flow path through gravel pack, ft  
+
| length of flow path through gravel pack, ft
 
|-
 
|-
|''L''<sub>''m''</sub>  
+
| ''L''<sub>''m''</sub>
|=  
+
| =
|length of matrix, ft  
+
| length of matrix, ft
 
|-
 
|-
|''L''<sub>''p''</sub>  
+
| ''L''<sub>''p''</sub>
|=  
+
| =
|length of perforation tunnel, ft  
+
| length of perforation tunnel, ft
 
|-
 
|-
|''L''<sub>''s''</sub>  
+
| ''L''<sub>''s''</sub>
|=  
+
| =
|length of damaged zone in fracture, ft  
+
| length of damaged zone in fracture, ft
 
|-
 
|-
|''L''<sub>''w''</sub>  
+
| ''L''<sub>''w''</sub>
|=  
+
| =
|completed length of horizontal well, ft  
+
| completed length of horizontal well, ft
 
|-
 
|-
|''L''<sub>''x''</sub>  
+
| ''L''<sub>''x''</sub>
|=  
+
| =
|distance from boundary, ft  
+
| distance from boundary, ft
 
|-
 
|-
|''m''  
+
| ''m''
|=  
+
| =
|162.2 ''qBμ''/''kh'' = slope of middle-time line, psi/cycle  
+
| 162.2 ''qBμ''/''kh'' = slope of middle-time line, psi/cycle
 
|-
 
|-
|''m''<sub>''B''</sub>  
+
| ''m''<sub>''B''</sub>
|=  
+
| =
|slope of bilinear flow graph, psi/hr<sup>1/4</sup>  
+
| slope of bilinear flow graph, psi/hr<sup>1/4</sup>
 
|-
 
|-
|''m''<sub>''L''</sub>  
+
| ''m''<sub>''L''</sub>
|=  
+
| =
|slope of linear flow graph, psi/hr<sup>1/2</sup>  
+
| slope of linear flow graph, psi/hr<sup>1/2</sup>
 
|-
 
|-
|''m''<sub>''s''</sub>  
+
| ''m''<sub>''s''</sub>
|=  
+
| =
|[[File:Vol5 page 0882 inline 001.png]], slope of spherical flow plot, psi-hr<sup>1/2</sup>  
+
| [[File:Vol5 page 0882 inline 001.png|RTENOTITLE]], slope of spherical flow plot, psi-hr<sup>1/2</sup>
 
|-
 
|-
|''m''<sub>''V''</sub>  
+
| ''m''<sub>''V''</sub>
|=  
+
| =
|slope of volumetric flow graph, psi/hr  
+
| slope of volumetric flow graph, psi/hr
 
|-
 
|-
|''m''<sub>hrf</sub>  
+
| ''m''<sub>hrf</sub>
|=  
+
| =
|slope of semilog plot for hemiradial flow, psi/log cycle  
+
| slope of semilog plot for hemiradial flow, psi/log cycle
 
|-
 
|-
|''m''<sub>elf</sub>  
+
| ''m''<sub>elf</sub>
|=  
+
| =
|slope of square-root-of-time plot for early linear flow, psi/[[File:Vol5 page 0882 inline 002.png]]  
+
| slope of square-root-of-time plot for early linear flow, psi/[[File:Vol5 page 0882 inline 002.png|RTENOTITLE]]
 
|-
 
|-
|''m''<sub>erf</sub>  
+
| ''m''<sub>erf</sub>
|=  
+
| =
|slope of semilog plot of early radial flow, psi/log cycle  
+
| slope of semilog plot of early radial flow, psi/log cycle
 
|-
 
|-
|''m''<sub>llf</sub>  
+
| ''m''<sub>llf</sub>
|=  
+
| =
|slope of square-root-of-time plot for late linear flow, psi/[[File:Vol5 page 0882 inline 002.png]]  
+
| slope of square-root-of-time plot for late linear flow, psi/[[File:Vol5 page 0882 inline 002.png|RTENOTITLE]]
 
|-
 
|-
|''m''<sub>prf</sub>  
+
| ''m''<sub>prf</sub>
|=  
+
| =
|slope of semilog plot for pseudoradial flow, psi/log cycle  
+
| slope of semilog plot for pseudoradial flow, psi/log cycle
 
|-
 
|-
|M  
+
| M
|=  
+
| =
|Molecular weight of gas  
+
| Molecular weight of gas
 
|-
 
|-
|MTR  
+
| MTR
|=  
+
| =
|middle-time region  
+
| middle-time region
 
|-
 
|-
|''n''  
+
| ''n''
|=  
+
| =
|inverse slope of the line on a log-log plot of the change in pressure squared or pseudopressure vs. gas flow rate  
+
| inverse slope of the line on a log-log plot of the change in pressure squared or pseudopressure vs. gas flow rate
 
|-
 
|-
|''p''  
+
| ''p''
|=  
+
| =
|pressure, psi  
+
| pressure, psi
 
|-
 
|-
|''p''<sub>avg</sub>  
+
| ''p''<sub>avg</sub>
|=  
+
| =
|average pressure, psi  
+
| average pressure, psi
 
|-
 
|-
|''p''<sub>''b''</sub>  
+
| ''p''<sub>''b''</sub>
|=  
+
| =
|base (atmospheric) pressure, psia  
+
| base (atmospheric) pressure, psia
 
|-
 
|-
|''p''<sub>0</sub>  
+
| ''p''<sub>0</sub>
|=  
+
| =
|arbitrary reference or base pressure, psi  
+
| arbitrary reference or base pressure, psi
 
|-
 
|-
|[[File:Vol5 page 0781 inline 001.png]]  
+
| [[File:Vol5 page 0781 inline 001.png|RTENOTITLE]]
|=  
+
| =
|volumetric average or static drainage-area pressure, psi  
+
| volumetric average or static drainage-area pressure, psi
 
|-
 
|-
|''p''<sub>''a''</sub>  
+
| ''p''<sub>''a''</sub>
|=  
+
| =
|adjusted or normalized pseudopressure, (''μz''/''p'')''p''<sub>''p''</sub>, psia  
+
| adjusted or normalized pseudopressure, (''μz''/''p'')''p''<sub>''p''</sub>, psia
 
|-
 
|-
|''p''<sub>''awf''</sub>  
+
| ''p''<sub>''awf''</sub>
|=  
+
| =
|adjusted flowing bottomhole pressure, psia  
+
| adjusted flowing bottomhole pressure, psia
 
|-
 
|-
|''p''<sub>''aws''</sub>  
+
| ''p''<sub>''aws''</sub>
|=  
+
| =
|adjusted shut-in bottomhole pressure, psia  
+
| adjusted shut-in bottomhole pressure, psia
 
|-
 
|-
|''p''<sub>''f''</sub>  
+
| ''p''<sub>''f''</sub>
|=  
+
| =
|formation pressure, psi  
+
| formation pressure, psi
 
|-
 
|-
|''p''<sub>''i''</sub>  
+
| ''p''<sub>''i''</sub>
|=  
+
| =
|original reservoir pressure, psi  
+
| original reservoir pressure, psi
 
|-
 
|-
|''p''<sub>''m''</sub>  
+
| ''p''<sub>''m''</sub>
|=  
+
| =
|matrix pressure, psi  
+
| matrix pressure, psi
 
|-
 
|-
|''p''<sub>''p''</sub>  
+
| ''p''<sub>''p''</sub>
|=  
+
| =
|pseudopressure, psia<sup>2</sup>/cp  
+
| pseudopressure, psia<sup>2</sup>/cp
 
|-
 
|-
|''p''<sub>''s''</sub>  
+
| ''p''<sub>''s''</sub>
|=  
+
| =
|stabilized shut-in BHP measured just before start of a deliverability test, psia  
+
| stabilized shut-in BHP measured just before start of a deliverability test, psia
 
|-
 
|-
|''p''<sub>''sc''</sub>  
+
| ''p''<sub>''sc''</sub>
|=  
+
| =
|standard-condition pressure, psia  
+
| standard-condition pressure, psia
 
|-
 
|-
|''p''<sub>''t''</sub>  
+
| ''p''<sub>''t''</sub>
|=  
+
| =
|surface pressure in tubing, psi  
+
| surface pressure in tubing, psi
 
|-
 
|-
|''p''<sub>''w''</sub>  
+
| ''p''<sub>''w''</sub>
|=  
+
| =
|BHP in wellbore, psi  
+
| BHP in wellbore, psi
 
|-
 
|-
|''p''<sub>''wf''</sub>  
+
| ''p''<sub>''wf''</sub>
|=  
+
| =
|flowing BHP, psi  
+
| flowing BHP, psi
 
|-
 
|-
|''p''<sub>''ws''</sub>  
+
| ''p''<sub>''ws''</sub>
|=  
+
| =
|shut-in BHP, psi  
+
| shut-in BHP, psi
 
|-
 
|-
|''p''<sub>''xy''</sub>  
+
| ''p''<sub>''xy''</sub>
|=  
+
| =
|parameter in horizontal well analysis equations  
+
| parameter in horizontal well analysis equations
 
|-
 
|-
|''p''<sub>''xyz''</sub>  
+
| ''p''<sub>''xyz''</sub>
|=  
+
| =
|parameter in horizontal well analysis equations  
+
| parameter in horizontal well analysis equations
 
|-
 
|-
|''p''<sub>''y''</sub>  
+
| ''p''<sub>''y''</sub>
|=  
+
| =
|parameter in horizontal well analysis equations  
+
| parameter in horizontal well analysis equations
 
|-
 
|-
|''p''<sub>1hr</sub>  
+
| ''p''<sub>1hr</sub>
|=  
+
| =
|pressure at 1-hour shut-in (flow) time on MTR line or its extrapolation, psi  
+
| pressure at 1-hour shut-in (flow) time on MTR line or its extrapolation, psi
 
|-
 
|-
|''p''′  
+
| ''p''′
|=  
+
| =
|pressure derivative
+
| pressure derivative
 
|-
 
|-
|''p''*  
+
| ''p''*
|=  
+
| =
|MTR pressure trend extrapolated to infinite shut-in time, psi  
+
| MTR pressure trend extrapolated to infinite shut-in time, psi
 
|-
 
|-
|''p''<sub>''D''</sub>  
+
| ''p''<sub>''D''</sub>
|=  
+
| =
|0.00708 ''kh''(''p''<sub>''i''</sub> – ''p'')/''qBμ'', dimensionless pressure as defined for constant-rate production  
+
| 0.00708 ''kh''(''p''<sub>''i''</sub> – ''p'')/''qBμ'', dimensionless pressure as defined for constant-rate production
 
|-
 
|-
|''p''<sub>MBHD</sub>  
+
| ''p''<sub>MBHD</sub>
|=  
+
| =
|Matthews-Brons-Hazebroek pressure, dimensionless  
+
| Matthews-Brons-Hazebroek pressure, dimensionless
 
|-
 
|-
|(''p''<sub>''D''</sub>)<sub>MP</sub>  
+
| (''p''<sub>''D''</sub>)<sub>MP</sub>
|=  
+
| =
|dimensionless pressure at match point  
+
| dimensionless pressure at match point
 
|-
 
|-
|''q''  
+
| ''q''
|=  
+
| =
|flow rate at surface, STB/D  
+
| flow rate at surface, STB/D
 
|-
 
|-
|''q''<sub>AOF</sub>  
+
| ''q''<sub>AOF</sub>
|=  
+
| =
|absolute-open-flow potential, MMscf/D  
+
| absolute-open-flow potential, MMscf/D
 
|-
 
|-
|''q''<sub>''g''</sub>  
+
| ''q''<sub>''g''</sub>
|=  
+
| =
|gas flow rate, Mscf/D  
+
| gas flow rate, Mscf/D
 
|-
 
|-
|''q''<sub>''o''</sub>  
+
| ''q''<sub>''o''</sub>
|=  
+
| =
|water flow rate, STB/D  
+
| water flow rate, STB/D
 
|-
 
|-
|''q''<sub>''Rt''</sub>  
+
| ''q''<sub>''Rt''</sub>
|=  
+
| =
|total flow rate at reservoir conditions, RB/D  
+
| total flow rate at reservoir conditions, RB/D
 
|-
 
|-
|''q''<sub>''sf''</sub>  
+
| ''q''<sub>''sf''</sub>
|=  
+
| =
|flow rate at formation (sand) face, STB/D  
+
| flow rate at formation (sand) face, STB/D
 
|-
 
|-
|''q''<sub>''w''</sub>  
+
| ''q''<sub>''w''</sub>
|=  
+
| =
|water flow rate, STB/D  
+
| water flow rate, STB/D
 
|-
 
|-
|''r''  
+
| ''r''
|=  
+
| =
|distance from the center of wellbore, ft  
+
| distance from the center of wellbore, ft
 
|-
 
|-
|''r''<sub>''a''</sub>  
+
| ''r''<sub>''a''</sub>
|=  
+
| =
|radius of altered zone (skin effect), ft  
+
| radius of altered zone (skin effect), ft
 
|-
 
|-
|''r''<sub>''d''</sub>  
+
| ''r''<sub>''d''</sub>
|=  
+
| =
|effective drainage radius, ft  
+
| effective drainage radius, ft
 
|-
 
|-
|''r''<sub>''dp''</sub>  
+
| ''r''<sub>''dp''</sub>
|=  
+
| =
|radius of damage zone around perforation tunnel, ft  
+
| radius of damage zone around perforation tunnel, ft
 
|-
 
|-
|''r''<sub>''e''</sub>  
+
| ''r''<sub>''e''</sub>
|=  
+
| =
|external drainage radius, ft  
+
| external drainage radius, ft
 
|-
 
|-
|''r''<sub>''i''</sub>  
+
| ''r''<sub>''i''</sub>
|=  
+
| =
|radius of investigation, ft  
+
| radius of investigation, ft
 
|-
 
|-
|''r''<sub>''p''</sub>  
+
| ''r''<sub>''p''</sub>
|=  
+
| =
|radius of perforation tunnel, ft  
+
| radius of perforation tunnel, ft
 
|-
 
|-
|''r''<sub>''s''</sub>  
+
| ''r''<sub>''s''</sub>
|=  
+
| =
|outer radius of the altered zone, ft  
+
| outer radius of the altered zone, ft
 
|-
 
|-
|''r''<sub>''sp''</sub>  
+
| ''r''<sub>''sp''</sub>
|=  
+
| =
|radius of source or inner boundary of spherical flow pattern, ft  
+
| radius of source or inner boundary of spherical flow pattern, ft
 
|-
 
|-
|''r''<sub>''w''</sub>  
+
| ''r''<sub>''w''</sub>
|=  
+
| =
|wellbore radius, ft  
+
| wellbore radius, ft
 
|-
 
|-
|''r''<sub>''wa''</sub>  
+
| ''r''<sub>''wa''</sub>
|=  
+
| =
|apparent or effective wellbore radius, ft  
+
| apparent or effective wellbore radius, ft
 
|-
 
|-
|''r''<sub>''D''</sub>  
+
| ''r''<sub>''D''</sub>
|=  
+
| =
|''r''/''r''<sub>''w''</sub>, dimensionless radius  
+
| ''r''/''r''<sub>''w''</sub>, dimensionless radius
 
|-
 
|-
|''R''<sub>''s''</sub>  
+
| ''R''<sub>''s''</sub>
|=  
+
| =
|dissolved GOR, scf/STB  
+
| dissolved GOR, scf/STB
 
|-
 
|-
|''s''  
+
| ''s''
|=  
+
| =
|skin factor, dimensionless  
+
| skin factor, dimensionless
 
|-
 
|-
|''s''<sub>''a''</sub>  
+
| ''s''<sub>''a''</sub>
|=  
+
| =
|skin caused by alteration of permeability around wellbore, dimensionless  
+
| skin caused by alteration of permeability around wellbore, dimensionless
 
|-
 
|-
|''s''<sub>''c''</sub>  
+
| ''s''<sub>''c''</sub>
|=  
+
| =
|convergence skin, dimensionless  
+
| convergence skin, dimensionless
 
|-
 
|-
|''s''<sub>''d''</sub>  
+
| ''s''<sub>''d''</sub>
|=  
+
| =
|skin caused by formation damage, dimensionless  
+
| skin caused by formation damage, dimensionless
 
|-
 
|-
|''s''<sub>''e''</sub>  
+
| ''s''<sub>''e''</sub>
|=  
+
| =
|skin caused by eccentric effects, dimensionless  
+
| skin caused by eccentric effects, dimensionless
 
|-
 
|-
|''s''<sub>''dp''</sub>  
+
| ''s''<sub>''dp''</sub>
|=  
+
| =
|perforation damage skin, dimensionless  
+
| perforation damage skin, dimensionless
 
|-
 
|-
|''s''<sub>''f''</sub>  
+
| ''s''<sub>''f''</sub>
|=  
+
| =
|skin of hydraulically fractured well, dimensionless  
+
| skin of hydraulically fractured well, dimensionless
 
|-
 
|-
|''s''<sub>''gp''</sub>  
+
| ''s''<sub>''gp''</sub>
|=  
+
| =
|skin factor from to Darcy flow through gravel pack, dimensionless  
+
| skin factor from to Darcy flow through gravel pack, dimensionless
 
|-
 
|-
|''s''<sub>min</sub>  
+
| ''s''<sub>min</sub>
|=  
+
| =
|minimum skin factor, dimensionless  
+
| minimum skin factor, dimensionless
 
|-
 
|-
|''s''<sub>''p''</sub>  
+
| ''s''<sub>''p''</sub>
|=  
+
| =
|skin resulting from an incompletely perforated interval, dimensionless  
+
| skin resulting from an incompletely perforated interval, dimensionless
 
|-
 
|-
|''s''<sub>''t''</sub>  
+
| ''s''<sub>''t''</sub>
|=  
+
| =
|total skin, dimensionless  
+
| total skin, dimensionless
 
|-
 
|-
|''s''<sub>''θ''</sub>  
+
| ''s''<sub>''θ''</sub>
|=  
+
| =
|skin factor resulting from well inclination, dimensionless  
+
| skin factor resulting from well inclination, dimensionless
 
|-
 
|-
|''s''′  
+
| ''s''′
|=  
+
| =
|''s'' + ''Dq'' = apparent skin factor, dimensionless  
+
| ''s'' + ''Dq'' = apparent skin factor, dimensionless
 
|-
 
|-
|''S''<sub>''g''</sub>  
+
| ''S''<sub>''g''</sub>
|=  
+
| =
|gas saturation, fraction of pore volume  
+
| gas saturation, fraction of pore volume
 
|-
 
|-
|''S''<sub>''o''</sub>  
+
| ''S''<sub>''o''</sub>
|=  
+
| =
|oil saturation, fraction of pore volume  
+
| oil saturation, fraction of pore volume
 
|-
 
|-
|''S''<sub>''w''</sub>  
+
| ''S''<sub>''w''</sub>
|=  
+
| =
|water saturation, fraction of pore volume  
+
| water saturation, fraction of pore volume
 
|-
 
|-
|''t''  
+
| ''t''
|=  
+
| =
|elapsed time, hours  
+
| elapsed time, hours
 
|-
 
|-
|''t''<sub>''a''</sub>  
+
| ''t''<sub>''a''</sub>
|=  
+
| =
|''μc''<sub>''t''</sub>''t''<sub>''ap''</sub>, adjusted or normalized pseudotime, hours  
+
| ''μc''<sub>''t''</sub>''t''<sub>''ap''</sub>, adjusted or normalized pseudotime, hours
 
|-
 
|-
|''t''<sub>''ap''</sub>  
+
| ''t''<sub>''ap''</sub>
|=  
+
| =
|pseudotime, hours  
+
| pseudotime, hours
 
|-
 
|-
|''t''<sub>''bD''</sub>  
+
| ''t''<sub>''bD''</sub>
|=  
+
| =
|dimensionless time in linear flow, hours  
+
| dimensionless time in linear flow, hours
 
|-
 
|-
|''t''<sub>''D''</sub>  
+
| ''t''<sub>''D''</sub>
|=  
+
| =
|0.0002637''kt''/''ϕμc''<sub>''t''</sub>''r''<sub>''w''</sub><sup>2</sup>, dimensionless time  
+
| 0.0002637''kt''/''ϕμc''<sub>''t''</sub>''r''<sub>''w''</sub><sup>2</sup>, dimensionless time
 
|-
 
|-
|''t''<sub>''DA''</sub>  
+
| ''t''<sub>''DA''</sub>
|=  
+
| =
|0.0002637 ''kt''/''ϕμc''<sub>''t''</sub>''A'' = dimensionless time based on drainage area, A  
+
| 0.0002637 ''kt''/''ϕμc''<sub>''t''</sub>''A'' = dimensionless time based on drainage area, A
 
|-
 
|-
|''t''<sub>''eqB''</sub>  
+
| ''t''<sub>''eqB''</sub>
|=  
+
| =
|equivalent time for bilinear flow, hours  
+
| equivalent time for bilinear flow, hours
 
|-
 
|-
|''t''<sub>''e''</sub>  
+
| ''t''<sub>''e''</sub>
|=  
+
| =
|equivalent time, hours  
+
| equivalent time, hours
 
|-
 
|-
|''t''<sub>''LfD''</sub>  
+
| ''t''<sub>''LfD''</sub>
|=  
+
| =
|0.0002637 ''kt''/''ϕμe''<sub>''t''</sub>''L''<sub>''f''</sub><sup>2</sup>, dimensionless time for fractured wells  
+
| 0.0002637 ''kt''/''ϕμe''<sub>''t''</sub>''L''<sub>''f''</sub><sup>2</sup>, dimensionless time for fractured wells
 
|-
 
|-
|''t''<sub>''p''</sub>  
+
| ''t''<sub>''p''</sub>
|=  
+
| =
|pseudoproducing time, hours  
+
| pseudoproducing time, hours
 
|-
 
|-
|''t''<sub>''pD''</sub>  
+
| ''t''<sub>''pD''</sub>
|=  
+
| =
|pseudoproducing time, dimensionless  
+
| pseudoproducing time, dimensionless
 
|-
 
|-
|''t''<sub>prf</sub>  
+
| ''t''<sub>prf</sub>
|=  
+
| =
|time required to reach the pseudoradial flow regime, hours  
+
| time required to reach the pseudoradial flow regime, hours
 
|-
 
|-
|''t''<sub>Eelf</sub>  
+
| ''t''<sub>Eelf</sub>
|=  
+
| =
|end of early linear flow, t, hours  
+
| end of early linear flow, t, hours
 
|-
 
|-
|''t''<sub>Eerf</sub>  
+
| ''t''<sub>Eerf</sub>
|=  
+
| =
|end of early radial flow, t, hours  
+
| end of early radial flow, t, hours
 
|-
 
|-
|''t''<sub>Elf</sub>  
+
| ''t''<sub>Elf</sub>
|=  
+
| =
|end of linear flow, hours  
+
| end of linear flow, hours
 
|-
 
|-
|''t''<sub>Ellf</sub>  
+
| ''t''<sub>Ellf</sub>
|=  
+
| =
|time to end of late linear flow regime, hours  
+
| time to end of late linear flow regime, hours
 
|-
 
|-
|''t''<sub>Ehrf</sub>  
+
| ''t''<sub>Ehrf</sub>
|=  
+
| =
|end of hemiradial flow, hours  
+
| end of hemiradial flow, hours
 
|-
 
|-
|''t''<sub>Erf</sub>  
+
| ''t''<sub>Erf</sub>
|=  
+
| =
|end of early radial flow, hours  
+
| end of early radial flow, hours
 
|-
 
|-
|''t''<sub>Eprf</sub>  
+
| ''t''<sub>Eprf</sub>
|=  
+
| =
|end of pseudoradial flow, hours  
+
| end of pseudoradial flow, hours
 
|-
 
|-
|''t''<sub>''p''</sub>  
+
| ''t''<sub>''p''</sub>
|=  
+
| =
|constant-rate production period, t, hours  
+
| constant-rate production period, t, hours
 
|-
 
|-
|''t''<sub>''pAD''</sub>  
+
| ''t''<sub>''pAD''</sub>
|=  
+
| =
|dimensionless producing time, hours  
+
| dimensionless producing time, hours
 
|-
 
|-
|''t''<sub>''pss''</sub>  
+
| ''t''<sub>''pss''</sub>
|=  
+
| =
|time required to reach pseudosteady state, hours  
+
| time required to reach pseudosteady state, hours
 
|-
 
|-
|''t''<sub>Self</sub>  
+
| ''t''<sub>Self</sub>
|=  
+
| =
|start of early linear flow, hours  
+
| start of early linear flow, hours
 
|-
 
|-
|''t''<sub>Sllf</sub>  
+
| ''t''<sub>Sllf</sub>
|=  
+
| =
|start of late linear flow, hours  
+
| start of late linear flow, hours
 
|-
 
|-
|''t''<sub>Shrf</sub>  
+
| ''t''<sub>Shrf</sub>
|=  
+
| =
|start of hemiradial flow, t, hours  
+
| start of hemiradial flow, t, hours
 
|-
 
|-
|''t''<sub>Sprf</sub>  
+
| ''t''<sub>Sprf</sub>
|=  
+
| =
|start of pseudoradial flow, t, hours  
+
| start of pseudoradial flow, t, hours
 
|-
 
|-
|''t''<sub>''s''</sub>  
+
| ''t''<sub>''s''</sub>
|=  
+
| =
|time required for stabilization, hours  
+
| time required for stabilization, hours
 
|-
 
|-
|''T''  
+
| ''T''
|=  
+
| =
|reservoir temperature, °R  
+
| reservoir temperature, °R
 
|-
 
|-
|''T''<sub>''sc''</sub>  
+
| ''T''<sub>''sc''</sub>
|=  
+
| =
|standard condition temperature, °R  
+
| standard condition temperature, °R
 
|-
 
|-
|''u''  
+
| ''u''
|=  
+
| =
|dummy variable  
+
| dummy variable
 
|-
 
|-
|''V''  
+
| ''V''
|=  
+
| =
|volume, bbl  
+
| volume, bbl
 
|-
 
|-
|''V''<sub>''f''</sub>  
+
| ''V''<sub>''f''</sub>
|=  
+
| =
|fraction of bulk volume occupied by fractures  
+
| fraction of bulk volume occupied by fractures
 
|-
 
|-
|''V''<sub>''m''</sub>  
+
| ''V''<sub>''m''</sub>
|=  
+
| =
|fraction of bulk volume occupied by matrix  
+
| fraction of bulk volume occupied by matrix
 
|-
 
|-
|''V''<sub>''w''</sub>  
+
| ''V''<sub>''w''</sub>
|=  
+
| =
|''V''<sub>''wb''</sub> = wellbore volume, bbl  
+
| ''V''<sub>''wb''</sub> = wellbore volume, bbl
 
|-
 
|-
|''w''  
+
| ''w''
|=  
+
| =
|width of channel reservoir, ft  
+
| width of channel reservoir, ft
 
|-
 
|-
|''w''<sub>''f''</sub>  
+
| ''w''<sub>''f''</sub>
|=  
+
| =
|fracture width, ft  
+
| fracture width, ft
 
|-
 
|-
|''wk''<sub>''f''</sub>  
+
| ''wk''<sub>''f''</sub>
|=  
+
| =
|fracture conductivity, md-ft  
+
| fracture conductivity, md-ft
 
|-
 
|-
|''w''<sub>''s''</sub>  
+
| ''w''<sub>''s''</sub>
|=  
+
| =
|width of damaged zone around fracture face, ft  
+
| width of damaged zone around fracture face, ft
 
|-
 
|-
|WBS  
+
| WBS
|=  
+
| =
|wellbore storage  
+
| wellbore storage
 
|-
 
|-
|''z''  
+
| ''z''
|=  
+
| =
|gas-law deviation factor, dimensionless  
+
| gas-law deviation factor, dimensionless
 
|-
 
|-
|[[File:Vol5 page 0862 inline 002.png]]  
+
| [[File:Vol5 page 0862 inline 002.png|RTENOTITLE]]
|=  
+
| =
|gas-law deviation factor at average reservoir pressure, dimensionless  
+
| gas-law deviation factor at average reservoir pressure, dimensionless
 
|-
 
|-
|Δ''p''  
+
| Δ''p''
|=  
+
| =
|pressure change since start of transient test, psi  
+
| pressure change since start of transient test, psi
 
|-
 
|-
|(Δ''p'')<sub>MP</sub>  
+
| (Δ''p'')<sub>MP</sub>
|=  
+
| =
|pressure change at match point  
+
| pressure change at match point
 
|-
 
|-
|Δ''p''<sub>''D''</sub>  
+
| Δ''p''<sub>''D''</sub>
|=  
+
| =
|dimensionless pressure change  
+
| dimensionless pressure change
 
|-
 
|-
|Δ''p''<sub>''p''</sub>  
+
| Δ''p''<sub>''p''</sub>
|=  
+
| =
|pseudopressure change since start of test, psia<sup>2</sup>/cp  
+
| pseudopressure change since start of test, psia<sup>2</sup>/cp
 
|-
 
|-
|Δ''p''<sub>''s''</sub>  
+
| Δ''p''<sub>''s''</sub>
|=  
+
| =
|additional pressure drop due to skin, psi  
+
| additional pressure drop due to skin, psi
 
|-
 
|-
|Δ''p''<sub>''t''=0</sub>  
+
| Δ''p''<sub>''t''=0</sub>
|=  
+
| =
|pressure drop at time zero, psi  
+
| pressure drop at time zero, psi
 
|-
 
|-
|Δ''p''<sub>1hr</sub>  
+
| Δ''p''<sub>1hr</sub>
|=  
+
| =
|pressure change from start of test to one hour elapsed time, psi  
+
| pressure change from start of test to one hour elapsed time, psi
 
|-
 
|-
|Δ''t''  
+
| Δ''t''
|=  
+
| =
|time elapsed since start of test, hours  
+
| time elapsed since start of test, hours
 
|-
 
|-
|Δ''t''<sub>''a''</sub>  
+
| Δ''t''<sub>''a''</sub>
|=  
+
| =
|[[File:Vol5 page 0884 inline 001.png]], normalized or adjusted pseudotime, hours  
+
| [[File:Vol5 page 0884 inline 001.png|RTENOTITLE]], normalized or adjusted pseudotime, hours
 
|-
 
|-
|Δ''t''<sub>''ap''</sub>  
+
| Δ''t''<sub>''ap''</sub>
|=  
+
| =
|[[File:Vol5 page 0884 inline 002.png]], pseudotime, hr-psia/cp  
+
| [[File:Vol5 page 0884 inline 002.png|RTENOTITLE]], pseudotime, hr-psia/cp
 
|-
 
|-
|Δ''t''<sub>''Be''</sub>  
+
| Δ''t''<sub>''Be''</sub>
|=  
+
| =
|bilinear equivalent time, hours  
+
| bilinear equivalent time, hours
 
|-
 
|-
|Δ''t''<sub>''e''</sub>  
+
| Δ''t''<sub>''e''</sub>
|=  
+
| =
|radial equivalent time, hours  
+
| radial equivalent time, hours
 
|-
 
|-
|Δ''t''<sub>''Le''</sub>  
+
| Δ''t''<sub>''Le''</sub>
|=  
+
| =
|linear equivalent time, hours  
+
| linear equivalent time, hours
 
|-
 
|-
|Δ''t''<sub>max</sub>  
+
| Δ''t''<sub>max</sub>
|=  
+
| =
|maximum shut-in time in pressure buildup test, hours  
+
| maximum shut-in time in pressure buildup test, hours
 
|-
 
|-
|Δ''V''  
+
| Δ''V''
|=  
+
| =
|change in volume, bbl  
+
| change in volume, bbl
 
|-
 
|-
|''η''  
+
| ''η''
|=  
+
| =
|0.0002637 ''k''/''ϕμc''<sub>t</sub>, hydraulic diffusivity, ft<sup>2</sup>/hr  
+
| 0.0002637 ''k''/''ϕμc''<sub>t</sub>, hydraulic diffusivity, ft<sup>2</sup>/hr
 
|-
 
|-
|''η''<sub>''fD''</sub>  
+
| ''η''<sub>''fD''</sub>
|=  
+
| =
|hydraulic diffusivity, dimensionless  
+
| hydraulic diffusivity, dimensionless
 
|-
 
|-
|''λ''  
+
| ''λ''
|=  
+
| =
|interporosity flow coefficient  
+
| interporosity flow coefficient
 
|-
 
|-
|''λ''<sub>''t''</sub>  
+
| ''λ''<sub>''t''</sub>
|=  
+
| =
|[[File:Vol5 page 0885 inline 001.png]], total mobility, md/cp  
+
| [[File:Vol5 page 0885 inline 001.png|RTENOTITLE]], total mobility, md/cp
 
|-
 
|-
|''α''  
+
| ''α''
|=  
+
| =
|exponent in deliverability equation  
+
| exponent in deliverability equation
 
|-
 
|-
|''α''  
+
| ''α''
|=  
+
| =
|parameter characteristic of system geometry in dual-porosity system  
+
| parameter characteristic of system geometry in dual-porosity system
 
|-
 
|-
|''β''  
+
| ''β''
|=  
+
| =
|turbulence factor  
+
| turbulence factor
 
|-
 
|-
|''β''′  
+
| ''β''′
|=  
+
| =
|transition parameter  
+
| transition parameter
 
|-
 
|-
|''γ''  
+
| ''γ''
|=  
+
| =
|Euler’s constant, = 1.781, dimensionless  
+
| Euler’s constant, = 1.781, dimensionless
 
|-
 
|-
|''γ''<sub>''g''</sub>  
+
| ''γ''<sub>''g''</sub>
|=  
+
| =
|gas gravity (air = 1.0)  
+
| gas gravity (air = 1.0)
 
|-
 
|-
|''γ''<sub>''m''</sub>  
+
| ''γ''<sub>''m''</sub>
|=  
+
| =
|matrix density  
+
| matrix density
 
|-
 
|-
|''ω''  
+
| ''ω''
|=  
+
| =
|storativity ratio in dual porosity reservoir  
+
| storativity ratio in dual porosity reservoir
 
|-
 
|-
|''μ''  
+
| ''μ''
|=  
+
| =
|viscosity, cp  
+
| viscosity, cp
 
|-
 
|-
|''μ''<sub>''i''</sub>  
+
| ''μ''<sub>''i''</sub>
|=  
+
| =
|viscosity evaluated at ''p''<sub>''i''</sub>, cp  
+
| viscosity evaluated at ''p''<sub>''i''</sub>, cp
 
|-
 
|-
|''μ''<sub>''g''</sub>  
+
| ''μ''<sub>''g''</sub>
|=  
+
| =
|gas viscosity, cp  
+
| gas viscosity, cp
 
|-
 
|-
|''μ''<sub>''o''</sub>  
+
| ''μ''<sub>''o''</sub>
|=  
+
| =
|oil viscosity, cp  
+
| oil viscosity, cp
 
|-
 
|-
|''μ''<sub>''w''</sub>  
+
| ''μ''<sub>''w''</sub>
|=  
+
| =
|water viscosity, cp  
+
| water viscosity, cp
 
|-
 
|-
|[[File:Vol5 page 0844 inline 002.png]]  
+
| [[File:Vol5 page 0844 inline 002.png|RTENOTITLE]]
|=  
+
| =
|gas viscosity evaluated at average pressure, cp  
+
| gas viscosity evaluated at average pressure, cp
 
|-
 
|-
|''μ''<sub>''gwf''</sub>  
+
| ''μ''<sub>''gwf''</sub>
|=  
+
| =
|gas viscosity evaluated at ''p''<sub>''wf''</sub> , cp  
+
| gas viscosity evaluated at ''p''<sub>''wf''</sub> , cp
 
|-
 
|-
|[[File:Vol5 page 0789 inline 001.png]]  
+
| [[File:Vol5 page 0789 inline 001.png|RTENOTITLE]]
|=  
+
| =
|viscosity evaluated at [[File:Vol5 page 0781 inline 001.png]], cp  
+
| viscosity evaluated at [[File:Vol5 page 0781 inline 001.png|RTENOTITLE]], cp
 
|-
 
|-
|''ρ''  
+
| ''ρ''
|=  
+
| =
|density, lbm/ft<sup>3</sup> or g/cm<sup>3</sup>  
+
| density, lbm/ft<sup>3</sup> or g/cm<sup>3</sup>
 
|-
 
|-
|''ρ''<sub>''wb''</sub>  
+
| ''ρ''<sub>''wb''</sub>
|=  
+
| =
|density of liquid in wellbore, lbm/ft<sup>3</sup>  
+
| density of liquid in wellbore, lbm/ft<sup>3</sup>
 
|-
 
|-
|''ϕ''<sub>''f''</sub>  
+
| ''ϕ''<sub>''f''</sub>
|=  
+
| =
|fraction of fracture volume occupied by pore space, ≅ 1  
+
| fraction of fracture volume occupied by pore space, ≅ 1
 
|-
 
|-
|''ϕ''<sub>''m''</sub>  
+
| ''ϕ''<sub>''m''</sub>
|=  
+
| =
|fraction of matrix volume occupied by pore space  
+
| fraction of matrix volume occupied by pore space
 
|-
 
|-
|(''ϕV'')<sub>''f''</sub>  
+
| (''ϕV'')<sub>''f''</sub>
|=  
+
| =
|fraction of bulk volume occupied by pore space in fractures  
+
| fraction of bulk volume occupied by pore space in fractures
 
|-
 
|-
|(''ϕVc''<sub>''t''</sub>)<sub>''f''</sub>  
+
| (''ϕVc''<sub>''t''</sub>)<sub>''f''</sub>
|=  
+
| =
|fracture "storativity" for dual porosity reservoir  
+
| fracture "storativity" for dual porosity reservoir
 
|-
 
|-
|(''ϕVc''<sub>''t''</sub>)<sub>''f+m''</sub>  
+
| (''ϕVc''<sub>''t''</sub>)<sub>''f+m''</sub>
|=  
+
| =
|total "storativity" for dual porosity reservoir  
+
| total "storativity" for dual porosity reservoir
 
|-
 
|-
|(''ϕV'')<sub>''m''</sub>  
+
| (''ϕV'')<sub>''m''</sub>
|=  
+
| =
|fraction of bulk volume occupied by pore space in matrix  
+
| fraction of bulk volume occupied by pore space in matrix
 
|-
 
|-
|''ϕ''  
+
| ''ϕ''
|=  
+
| =
|porosity, dimensionless  
+
| porosity, dimensionless
 
|-
 
|-
|Σ''s''  
+
| Σ''s''
|=  
+
| =
|sum of damage skin, turbulence, and other pseudoskin factors  
+
| sum of damage skin, turbulence, and other pseudoskin factors
 
|}
 
|}
  
==References==
+
== References ==
<references>
 
<ref name="r1">Goode, P.A. and Thambynayagam, R.K.M. 1987. Pressure Drawdown and Buildup Analysis of Horizontal Wells in Anisotropic Media. ''SPE Form Eval'' '''2''' (4): 683–697. SPE-14250-PA. http://dx.doi.org/10.2118/14250-PA </ref>
 
  
<ref name="r2">Odeh, A.S. and Babu, D.K. 1990. Transient Flow Behavior of Horizontal Wells: Pressure Drawdown and Buildup Analysis. SPE Form Eval 5 (1): 7-15. SPE-18802-PA. http://dx.doi.org/10.2118/18802-PA </ref>
+
<references />
  
<ref name="r3">Kuchuk, F.J. 1995. Well Testing and Interpretation for Horizontal Wells. ''J Pet Technol'' '''47''' (1): 36–41. SPE-25232-PA. http://dx.doi.org/10.2118/25232-PA </ref>
+
== Noteworthy papers in OnePetro ==
  
<ref name="r4">Lichtenberger, G.J. 1994. Data Acquisition and Interpretation of Horizontal Well Pressure-Transient Tests. ''J Pet Technol'' '''46''' (2): 157-162. SPE-25922-PA. http://dx.doi.org/10.2118/25922-PA </ref>
+
Use this section to list papers in OnePetro that a reader who wants to learn more should definitely read
</references>
 
  
==Noteworthy papers in OnePetro==
+
== External links ==
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
 
Use this section to provide links to relevant material on websites other than PetroWiki and OnePetro
  
==See also==
+
== See also ==
[[Type curves]]
+
 
 +
[[Type_curves|Type curves]]
 +
 
 +
[[Estimating_horizontal_well_productivity|Estimating horizontal well productivity]]
  
[[Estimating horizontal well productivity]]
+
[[Coning|Coning]]
  
[[Coning]]
+
[[Pressure_transient_testing|Pressure transient testing]]
  
[[Pressure transient testing]]
+
[[Horizontal_wells|Horizontal wells]]
  
[[Horizontal wells]]
+
[[Fluid_flow_through_permeable_media|Fluid flow through permeable media]]
  
[[Fluid flow through permeable media]]
+
[[PEH:Fluid_Flow_Through_Permeable_Media]]
  
[[PEH:Fluid Flow Through Permeable Media]]
+
[[Category:5.3 Reservoir fluid dynamics]]

Latest revision as of 19:01, 4 June 2015

Productivity estimates in horizontal wells are subject to more uncertainty than comparable estimates in vertical wells. Further, it is much more difficult to interpret well test data because of 3D flow geometry. The radial symmetry usually present in a vertical well does not exist. Several flow regimes can potentially occur and need to be considered in analyzing test data from horizontal wells. Wellbore storage effects can be much more significant and partial penetration and end effects commonly complicate interpretation.

In vertical wells, variables such as average permeability, net vertical thickness, and skin are used. Horizontal wells need more detail. Not only is vertical thickness important, but the horizontal dimensions of the reservoir, relative to the horizontal wellbore, need to be known.

Steps in evaluating horizontal well-test data

There are three basic steps in evaluating pressure-transient data from a horizontal well:

  1. Identify the specific flow regimes in the test data
  2. Apply the proper analytical and graphical procedures to the data
  3. Evaluate the uniqueness and sensitivity of the results to properties derived from analysis or simply assumed

Identify flow regimes

Evaluation of data from a vertical wellbore will generally center on a single flow regime, such as infinite-acting radial flow, known as the MTR. However, a pressure-transient test in a horizontal well can involve as many as five major and distinct regimes that need to be identified. These regimes may or may not occur in a given test and may or may not be obscured by wellbore storage effects.

Apply the proper procedures

Each flow regime can be modeled by an equation that can be used to estimate important reservoir properties. At best, only groups of analytical parameters can be determined directly from equations. It is imperative that the proper analytical and graphical procedures be applied to the data. In many cases, when solving for specific parameters, the application of these analytical expressions may involve a complex iterative procedure.

Evaluate uniqueness and sensitivity

Experience indicates that results of horizontal well test analysis are seldom unique, so it is important that the uniqueness and sensitivity of the results to assumed properties be evaluated. Simulation of the test using properties that have been determined from the test can confirm that at least the analysis is consistent with the test data. A simulator can also determine whether other sets of formation properties will also lead to a fit of the data.

Horizontal well flow regimes

Different formation properties can be calculated from the data in each of the five different flow regimes. Any flow regime may be absent from a plot of test data because of geometry, wellbore storage, or other factors. Nor does the fact that they can appear mean that they do appear. The five different flow regimes that can occur are:

  • Early radial
  • Hemi-radial
  • Early linear
  • Late pseudoradial
  • Late linear

Fig. 1 shows a horizontal well with length, Lw, within a reservoir that is assumed to be a rectangular parallelepiped or a "box reservoir" drainage area. In this discussion, it is assumed that the axes of the coordinate system coincide with the direction of principal permeability and the well produces over its entire length, Lw.

The axes for this box are the usual x-, y-, and z-axes. Notice that the x-axis is measured along the bottom edge of the reservoir, going from left to right in the direction perpendicular to the well. The y-axis lies along the axis from front to back of the reservoir, parallel to the wellbore. The z-axis is oriented in the direction of reservoir thickness.

The total width of the reservoir perpendicular to the wellbore is aH, the total length in a direction parallel to the wellbore is bH, and the total height of the reservoir is the net pay thickness, h. Notice the parameters for the distance from the well to the various borders. Along the axis of the well, the shortest distance from the end of the well to a boundary is dy, and the longest distance from the other end of the well to the boundary is Dy. In the vertical direction, the shortest distance to a vertical boundary is dz, and the longest distance to a vertical boundary is Dz.

Characteristics of flow regimes

Consider a well producing at a constant rate. The early radial flow regime occurs before the area drained or the pressure transient caused by this production encounters either of the vertical boundaries of the reservoir. Fig. 2 shows a radial flow pattern penetrating out into the reservoir. Actually, however, this flow pattern is likely to be elliptical, moving further into the reservoir at a given time in the higher-permeability x-direction than in the lower-permeability z-direction. This phenomenon causes no significant complications in our analysis.

When the wellbore is much nearer one vertical boundary than the other, another flow regime, called half-radial or hemiradial flow (Fig. 3) may exist. Hemiradial flow can occur immediately following the early radial flow regime, if the well is much nearer one of the vertical boundaries than the other. Eventually, the area affected by the production will include the entire thickness of the reservoir. When that happens, a linear flow pattern may develop, as Fig. 4 shows.

Eventually, flow will begin to come into the wellbore from beyond the ends of the well. Until these end effects become important, early linear flow continues. Once end effects become important a transition period is followed by a later pseudoradial flow regime, as Fig. 5 illustrates.

This flow regime continues until the area affected by the production reaches one of the sides of the reservoir. Once the area affected is the entire width of the reservoir (that is, the pressure transient has reached both sides of the reservoir), then the late-linear flow regime begins (Fig. 6).

Identifying flow regimes in horizontal wells

All of these flow regimes in a test can be identified on a diagnostic log-log plot of the pressure change, Δp, and pressure derivative, p′, against the logarithm of time (Fig. 7).

A unit-slope line appears during wellbore storage; a horizontal derivative during early radial flow, and then, later, in pseudoradial flow; and a half-slope line in early-linear flow and then in late-linear flow. (These half-slope lines appear on the derivative but not on the pressure-change curves.) This does not imply that all these flow regimes will appear in any given test; in fact, that would be rare. But these are the shapes that identify the flow regimes that may appear in the test being analyzed.

The shapes that may appear in a drawdown test (which is the basis of Fig. 7) may not appear in a buildup test because of the complex superposition of flow regimes. For example, a test would have to be in linear flow both at time (tp + Δt) and at time Δ t to ensure appearance of a derivative with half-slope; this is highly unlikely. The best way to solve the problem is to ensure that a buildup test on a horizontal well is run with a producing time, tp, much greater than the maximum shut-in time in the test (that is, tp > 10 Δt max ).

Table 1 summarizes the working equations for permeability, skin, and start and end of each of the recognized flow regimes. Different investigators have found different equations for start and end of various flow regimes, especially linear flow regimes. This is partly because of a difference of assumptions about flow into the wellbore. Uniform flux or infinite conductivity models are common; neither is rigorously applicable in practice. [1] In this section, the equations for duration of flow regimes derived by Odeh and Babu[2] are used. This model assumes uniform flux into the wellbore.

Early-radial flow

Early-radial flow is similar to the radial flow period in a vertical well (Fig. 2). The governing equation for this flow regime is

RTENOTITLE....................(1)

Data for this period may be masked by wellbore storage effects, but, when present, they may be analyzed on a semilog plot.

The early-radial flow regime may in theory start at time zero, in absence of wellbore storage effects. The end of the early-radial flow regime may occur when the transient reaches a vertical boundary or when flow comes from beyond the end of the wellbore. The end of the period is the smaller of these two values.

Eq. 2[2] says that the period must end when the transient reaches the nearest boundary, dz, from the well. This equation includes the permeability in the vertical direction:

RTENOTITLE....................(2)

The radial flow regime may also end when flow from beyond the end of the wellbore becomes important. Eq. 3 gives the time by

RTENOTITLE....................(3)

Lw is the completed length of the well, and k y is the permeability in the direction parallel to the wellbore. The actual end is the lesser of the two times calculated from Eqs. 2 and 3. It is helpful to check the expected duration of the early-radial flow regime after estimating the parameters necessary to make these calculations.

Eq. 1 suggests that possible radial flow on the diagnostic plot be identified and then bottomhole flowing pressure be plotted against time during the appropriate time range on semilog coordinates. The slope of the straight line that results is

RTENOTITLE....................(4)

The group RTENOTITLE can be found from the slope, merf:

RTENOTITLE....................(5)

Effective completed length of the well must be known to make this calculation. This is not necessarily the same as the perforated or completed length of the well. Some sections of the well may not produce at all.

The equation for calculating the altered permeability skin, sd, for early-radial flow is

RTENOTITLE....................(6)

When analyzing a buildup test rather than a constant-rate flow test, plot the HTR or equivalent time on the horizontal axis of the semilog plot, and then plot shut-in or equivalent time on the vertical axis. Note that this plotting is correct only if (tp + Δt) and Δt appear in this time period simultaneously; that is, radial flow must exist at both time (tp + Δt) and time Δt. This is unlikely because radial-flow regime may exist at time Δt, but a different flow regime is likely at time (tp + Δt).

Example 1: Well Erf-1

For drawdown test data from Well Erf-1, [2] the diagnostic plot indicates the data from approximately 0.24 to 24 hours may be in early-radial flow. The following information is available for this well: q = 800 STB/D, μ = 1 cp, B = 1.25 RB/STB, rw = 0.25 ft, ϕ = 0.2, ct = 15×10 –6 psi–1, centered in box-shaped drainage area, h = 200 ft, bH = 4,000 ft, and aH = 2,000 ft, Lw =1,000 ft, and, from analysis of data in early linear flow regime, kx = 200 md. Table 2 shows the pressure change data for 0.24 to 24 hours.

Plot (pipwf ) = Δp vs. t on semilog coordinates (Fig. 8). The plot results in a straight line with a slope of 8 psi/cycle. In Fig. 8, at t = 2.4 hours, the points begin to deviate from the straight line, as expected from calculations for flow regime duration that follow. The pressure change at 1 hour is 39 psia. Using the slope of 8 and Eq. 5,

RTENOTITLE

Thus, because kx = 200 md, kz = 2 md. Using the value of 39 for Δp1hr from Fig. 8, skin from Eq. 6 is

RTENOTITLE

The start of the early-radial flow regime is controlled by wellbore storage, which appears to have vanished at times earlier than 0.24 hours in this example. The end of the early radial flow regime is expected at the lesser of the two values derived from Eqs. 2 and 3. For a centered well, dz = h/2 = 100 ft, and Eq. 2 gives

RTENOTITLE

Assuming ky = kx = 200 md, Eq. 3 gives

RTENOTITLE

Thus, expect the early-radial flow regime to end at approximately 1.875 hours, which is the smaller value and is consistent with observed test data.

Hemiradial flow

The hemiradial flow period (Fig. 3) will occur only when the well is close to one of the vertical boundaries (either the upper or the lower boundaries) and is analogous to a vertical well near a fault. The governing equation is[1]

RTENOTITLE....................(7)

A horizontal derivative on the diagnostic plot identifies hemiradial flow. If data appear to fall into this flow regime, a straight line on a semilog plot would provide more confidence that radial flow has been identified. Consistency checks in the analysis coupled with a well survey will be required to distinguish hemiradial flow from early radial flow.

The time range in which the analysis for hemiradial flow is valid begins after the closest vertical boundary, dz, affects the data and before the farthest boundary, Dz, affects them. In the absence of wellbore storage, the start of hemiradial flow is given by

RTENOTITLE....................(8)

Note that the start of the hemiradial flow regime involves the shortest distance to a vertical boundary and the permeability in the vertical direction. However, wellbore storage will most likely determine the actual start of hemiradial flow.

The end of hemiradial flow occurs when pressure is affected by the farther vertical boundary or flow from beyond the ends of the wellbore, whichever occurs first. It is the smaller of the times calculated using Eqs. 9 and 10. If the hemiradial flow regime ends when pressure reaches the farthest vertical boundary, it depends on the distance, Dz, and the vertical permeability, kz:

RTENOTITLE....................(9)

When the appearance of end effects—flow from beyond the ends of the wellbore—causes the end of the hemiradial flow regime to appear, the end of the flow regime occurs when

RTENOTITLE....................(10)

The completed length of the well, Lw, and the permeability ky, parallel to the wellbore appear in this equation. These parameters determine when enough flow has come from beyond the ends of the wellbore to distort the radial flow pattern that appeared earlier.

RTENOTITLE....................(11)

gives the slope of the semilog straight line for semiradial flow, mhrf. The multiplier, 325.2, is twice the multiplier for early-radial flow. The equation to estimate the damage skin factor is also similar to that for radial flow but has a multiplier that differs by a factor of two.

RTENOTITLE....................(12)

The equations relating slope and permeability and the equation for skin are similar in a buildup test to those for a drawdown test. The pressure change in the equation for skin is [p1hrpwft = 0)]. Semilog plots of buildup test data from the hemiradial flow regime cannot be analyzed rigorously using data from a Horner plot unless the pressure data at (tp + Δt) and at time Δt are simultaneously in this flow regime. As a practical matter, the hemiradial flow regime is likely to appear clearly in the buildup test only when the producing time is much greater that the shut-in time.

Early linear flow

The governing equation for early-linear flow is[1]

RTENOTITLE....................(13)

The "convergence skin," sc, is discussed later in this section. The start of the early-linear flow regime (Fig. 4) depends on the farthest distance to a vertical boundary, Dz, and the vertical permeability, kz.[2]

RTENOTITLE....................(14)

Not until flow reaches that farthest vertical boundary can a linear flow pattern begin toward the well. This flow period ends when fluids flow from beyond the ends of the wellbore. Thus,

RTENOTITLE....................(15)

Notice that the end depends on the effective completed length of the well, Lw, and on the permeability in the direction parallel to the well. This is the time in which end effects—flow beyond the ends of the well—begin to significantly distort the linear flow pattern.

The early-linear flow regime is identified in a drawdown test with a half-slope for the derivative. (Because of the skin effect, the pressure change curve on the diagnostic plot will only approach a half-slope asymptotically.) For data identified as being in this flow regime, plot pressure against the square root of time.

The slope of the straight line on such a plot, melf, can be used to estimate the square-root of kx, the horizontal permeability perpendicular to the well:

RTENOTITLE....................(16)

To calculate the damage skin,

RTENOTITLE....................(17)

This equation includes a convergence skin, sc, which is[2]

RTENOTITLE....................(8.203)

This convergence skin is an additional pressure drop that acts like a skin effect caused by flow moving from throughout the entire formation until it converges down to the small wellbore in the middle of the formation (Fig. 9). This convergence skin is defined in terms of the ratio of the permeability in the x-direction, which is perpendicular to the wellbore, to the vertical permeability. It also involves the distance to the nearest vertical boundary, dz, and the net pay thickness, h.

Kuchuk[3] derived a different equation for convergence skin. As a practical matter, the Odeh-Babu and Kuchuk equations lead to the same result. When there has been a single rate preceding shut-in during early-linear flow, the buildup pressure is plotted against RTENOTITLE on a Cartesian plot, which is sometimes called a tandem root plot. The permeability, kx, is calculated from the slope, melf, of the plot and Eq. 16. kx has the same relationship to the slope that existed in a drawdown test. Skin for this flow regime is calculated with Eqs. 17 and 18.

Plots of buildup data from the early-linear flow regime cannot be analyzed rigorously with a plot of pws vs. RTENOTITLE (that is, the tandem-root plot) unless data at (tp + Δt) and Δt are simultaneously in this flow regime—highly unlikely—or unless tp is much greater than Δt, in which case simply plot pws vs. RTENOTITLE. Little error results from ignoring the (tp + Δt) term, which is essentially constant.

Example 2: Well Elf-2

The diagnostic plot for a drawdown test from Well Elf-2[2] indicates data in the early-linear flow regime because the derivative has a half-slope. The following data apply to this well: q = 800 STB/D; μ = 1 cp; B = 1.25 RB/STB; rw = 0.25 ft; ϕ = 0.2; ct = 15×10–6 psi–1; centered in box-shaped drainage area 100 ft thick, 4,000 ft long, 4,000 ft wide; Lw = 2,500 ft; and, from early radial-flow regime data, kxkz = 8,000 md2. In addition, Table 3 shows the pressure-change data for this well.

A plot of pressure change vs. the square root of time (Fig. 10) indicates early linear flow. The final point on the straight line is at a time of approximately 24 hours, but there may have been some deviation from the straight line by this time. From the slope of the straight line, melf = 0.934 psi/hr1/2 and Eq. 16, calculate the permeability in the horizontal plane perpendicular to the wellbore.

RTENOTITLE

or kx = (20.1)2 = 404 md.

Analysis of data from the early radial flow regime indicated that kxkz is 8,000 md2; thus, kz is approximately 19.8 md. To calculate sd, use Eqs. 17 and 18, noting that the value for (pipwf ) = 3.1 at t = 0.

RTENOTITLE RTENOTITLE

Then, sd = 4.91 – 4.91 = 0.

Check the expected time range for the early-linear flow regime.

Using Dz = h/2 = 50 ft and kz = 20 md, calculate the beginning of linear flow with Eq. 14:

RTENOTITLE

Assuming ky = kx at 400 md, use Eq. 15 to find the end of early linear flow.

RTENOTITLE

These limits are reasonably consistent with the time range analyzed assuming early linear flow.

Late pseudoradial

The governing equation for late-pseudoradial flow is[1]

RTENOTITLE....................(19)

The late-pseudoradial flow period occurs only if [2]

RTENOTITLE....................(20)

Here, bH is the dimension of the reservoir parallel to the wellbore. As long as the completed length of the well is relatively short compared with the length of the drainage area late-pseudoradial flow can occur.

The start of this flow period occurs when fluid flows from well beyond the ends of the wellbore (Fig. 5). It is approximated with[2]

RTENOTITLE....................(21)

This starting time depends on the completed length of the well, Lw, and on the permeability in the direction of the well, ky. The end of this period, like others in this section, is approximated by the minimum of the results of two calculations. The first,

RTENOTITLE....................(22)

depends on dy and the length of the wellbore along with ky, the permeability in the direction parallel to the wellbore. This is the time at which horizontal boundary effects first appear.

The other equation gives a time at which the radial-flow pattern begins to be distorted depending on the shortest distance, dx, from the well to a boundary perpendicular to the wellbore and on kx, the permeability in that direction.

RTENOTITLE....................(23)

Whenever boundary effects first appear, whether in a direction that is parallel to the well or perpendicular to the axis of the well, the late-pseudoradial flow period will end.

The diagnostic plot helps identify the late-pseudoradial flow regime with the characteristic horizontal derivative. For data in the appropriate time range, prepare a semilog plot of pressure against time for a drawdown test. The slope of this plot will be mprf and the relationship between that slope and the square root of kxky, or the permeabilities in the horizontal plane, is given by

RTENOTITLE....................(24)

The skin equation is similar in form to those seen before:

RTENOTITLE....................(25)

Again, the total skin depends on Δp1hr. The convergence skin (Eq. 18) is subtracted from the "total" skin to determine the damage skin.

For a buildup test preceded by production at a single rate, plot pressure against the HTR on a semilog graph. Permeability is calculated from Eq. 24, the same as for a drawdown test. The skin equation is basically the same as for a drawdown test, except that the Δp1hr is now p1hrpwf. To obtain p*, the extrapolated pressure, extrapolate the semilog straight line to a HTR of unity.

Semilog plots of buildup test data from the late-pseudoradial flow regime cannot be analyzed rigorously using a Horner plot unless pressures at (tp + Δt) and at time Δt are simultaneously in the pseudoradial flow regime, which is highly unlikely. However, little error appears if the producing time before shut-in is much greater than the maximum shut-in time achieved in the buildup test.

Example 3: Well Prf-3

The diagnostic plot suggests that a constant-rate drawdown test from Well Prf-3[2] includes data in the late-pseudoradial period. The following data are available from the test: q = 800 STB/D, μ = 1 cp, B = 1.25 RB/STB, rw = 0.25 ft, ϕ = 0.2, ct = 15×10–6 psi–1, h = 150 ft, Lw = 900 ft, aH = 5,280 ft, bH = 5,280 ft, well centered in drainage volume, kx = 100 md (from analysis of early linear flow), kxkz = 1,000 md2, and kz = 10 md (from analysis of early radial flow). Table 4 gives the pressure change, Δp = pipwf vs. time.

A plot of pressure change vs. the logarithm of time (Fig. 11) confirms pseudoradial flow. A straight line fits all the data from 192 to 432 hours; the slope of the line, mprf, is 15.3 psi/cycle, and Δp1 hr = 18.94 psi (extrapolated). Then, from Eq. 24,

RTENOTITLE

Thus, RTENOTITLE

From Eqs. 25 and 18 ,

RTENOTITLE

Here,

RTENOTITLE

The pseudoradial flow regime should start at the time given by Eq. 21:

RTENOTITLE

It should end at a time given by the lesser of values from Eqs. 22 and 23. From Eq. 22,

RTENOTITLE

where dy = 1/2(5,200-900) = 2,190ft for this centered well. From Eq. 23,

RTENOTITLE

The smaller of these two values is 344 hours, which is thus the expected end of pseudoradial flow. The data on the figure that lie on the straight line show the time range from 192 to 432 hours, which is generally consistent with the expected duration of the flow regime.

Late-linear flow

The governing equation for late-linear flow is[1][2]

RTENOTITLE....................(26)

The late-linear flow regime starts after the pressure transient has reached the boundaries in the z- and y-directions, and the flow behavior with regard to these directions has become pseudosteady state, as Fig. 6 shows.

The start of this time period is the maximum of two equations. [2] The first depends on the time to reach the boundary, Dy, beyond the end of the horizontal well. It also depends on the permeability, ky, in the direction parallel to the wellbore.

RTENOTITLE....................(27)

Another requirement for the start of the late-linear flow regime is the time to reach the maximum vertical distance, Dz, divided by the vertical permeability:

RTENOTITLE....................(28)

Usually, the start of the late-linear flow regime is dictated by the time to reach the boundaries in the y-direction. The end of this period is given by the equation

RTENOTITLE....................(29)

The end of the late-linear flow regime depends on reaching the nearest boundary in the direction perpendicular to the wellbore, which is the distance, dx, away, and on the permeability kx in that direction.

Identify the late-linear flow regime by a half-slope on the derivative in the diagnostic plot of drawdown test data. (The pressure change may approach a half-slope asymptotically.) For data that appear to be in this flow regime, plot pressure against the square root of time. From the slope mllf of the plot, estimate permeability in the x-direction from

RTENOTITLE....................(30)

Alternatively, if kx is known from an early-linear flow regime, estimate bH, the length of the drainage area, from

RTENOTITLE....................(31)

This late-linear flow regime is the only period that provides the data to calculate the total skin, s, including the partial-penetration skin, sp, and the convergence skin, sc. To calculate the damage skin, sd, use

RTENOTITLE....................(32)

The total skin depends on Δpt = 0. Subtracting the partial penetration skin, sp, and the convergence skin, sc, from the total skin yields the damage skin.

The partial-penetration skin is a complex function that is calculated with Eqs. A-25 through A-35 in Table 1. For a buildup test, plot pressure against the HTR. From the slope, mllf, calculate kx with Eq. 30, exactly the same as for drawdown tests. Or, if kx is known, estimate the length, bH, of the drainage area with Eq. 31. Calculate the damage skin, sd, from a pressure buildup test from Eq. 32, where Δpt = 0 = (pt = 0)extpwf(t = 0).

Note that the same difficulty arises in using superposition to find plotting functions plots of buildup data from the late-linear flow regime as existed with the previous flow regimes. Pressures at both time (tp + Δt) and time Δt must be in the late linear flow regime for a tandem-root plot to be valid. However, if tp >> Δtmax, there is little error.

Example 4: Well Llf-4

The diagnostic plot for a drawdown test from Well Llf-4[2] appears to include data in the late-linear flow regime (derivative with half-slope). The following data applies to this well: q = 800 STB/D, μ = 1 cp, B = 1.25 RB/STB, rw = 0.25 ft, ϕ = 0.2, ct = 15 × 10–6 psi–1 , h = 150 ft, Lw = 1,000 ft, bH = 2,000 ft (well centered), aH = 6,968 ft (well centered), Dz = 85 ft, dz = 65 ft, kxkz = 1,000 md2 (from analysis of early-radial flow), and kxky = 5,000 md2 (from analysis of pseudoradial flow). Table 5 gives pressure change, Δp = pipwf , data vs. time.

Fig. 12 is a plot of pressure change vs. the square root of time. The straight line on this plot for the entire time range (60 to 240 hours) confirms late-linear flow for this range. The slope of the line is 1.56 psi/hr1/2, and the intercept is Δpt = 0 =28.4 psi.

From Eq. 30,

RTENOTITLE

Then, kx = 100 md. Because kkkz = 100 md2, kz = 10 md. Also, because kxky = 5,000 md2, ky =50 md. From Eq. 32,

RTENOTITLE

From Eq. 18,

RTENOTITLE

Calculate the partial penetration skin, sp , using the appropriate equation from among Eqs. A-25 through A-35 in Table 1. First, calculate

RTENOTITLE

to determine whether "Case 1" or "Case 2" applies:

RTENOTITLE

Because

RTENOTITLE

this is Case 1 (Eq. A-26). Use Eqs. A-25 through A-31 from the table. From Eq. A-27, sp = pxyz + pxy. From Eq. A-25,

RTENOTITLE

From Eq. A-28,

RTENOTITLE

Here, from Eq. A-29,

RTENOTITLE

The well is centered, so dy = Dy = 500 ft.

RTENOTITLE

From Eq. A-30,

RTENOTITLE

Also,

RTENOTITLE

From Eq. A-31,

RTENOTITLE

and

RTENOTITLE

Then,

RTENOTITLE

Then,

RTENOTITLE

Now check the expected duration of the late-linear flow regime. The start is the larger of values from Eqs. 27 and 28. From Eq. 27,

RTENOTITLE

From Eq. 28,

RTENOTITLE

Thus, the start is expected to be at approximately 162 hours. Eq. 29 gives the end of the flow regime.

RTENOTITLE

The data in this example spanned the time range from 60 to 240 hours. Some of the data that fall on the straight line are, in theory, from times before the start of the late-linear flow regime, but they appear to cause no problem in determining the slope of the line.

Field examples[4]

The following field examples illustrate the procedures used in analyzing horizontal well-test data.

Field example Well A

Table 6 summarizes the reservoir and completion properties for Well A. The target for Well A, a horizontal exploration well, was vertical tectonic fracture development in a low-permeability shale. Because of the fractures, the permeability is assumed to be isotropic (kh = kz) and a result of the fractures. Fig. 13 is a diagnostic plot for Well A and includes a history match using an analytical model.

During the early part of the test there is a unit-slope line representing wellbore storage. Following that wellbore storage, there is a transition to radial flow. The final few data points may be in radial flow. On the Horner plot (Fig. 14), the last few data points fall on a straight line. From the slope of the straight line, the apparent permeability is 0.011 md and the altered zone skin is 2.9. There is no evidence of boundary effects on this Horner plot. The existence of the semilog straight line is not assured, but the data are at least on the verge of reaching it.

Guided by the Horner analysis results, engineers simulated these data with an analytical horizontal well model. The initial match was poor. The match was improved considerably by introducing a no-flow boundary approximately 16 ft above the wellbore, which led to a permeability estimate of 0.027 md and an apparent skin of 11.5. It was concluded that the flow regime observed in the test was hemiradial flow.

The final match, shown on the type-curve plot in Fig. 13, is still not a good match at all times, but the author stated that the poor match in the transition region could be the result of phase-redistribution effects in the wellbore. The distance to the no-flow boundary that led to the best match compared favorably with well survey data, which indicated that the well was drilled approximately 20 ft below the upper limit of the productive horizon.

Well B is in a west Texas carbonate formation. It was expected to have isotropic permeability caused by fracturing and dissolution. Table 7 gives the field data for this well.

Fig. 15 is the diagnostic plot for this well. After wellbore storage, a short period of radial flow appears to be followed by the onset of linear flow, because p′ approaches a slope of 0.5. In the time region where the derivative is horizontal, a straight line on the Horner plot (Fig. 16) yields k = 0.14 md. Using these results in the analysis shows that the end of radial flow occurs at tErf = 165 hours.

A tandem-root plot (Fig. 17) indicates linear flow and also suggests a distance to the nearest boundary of 29 ft. This is in good agreement with geological observations and helps to verify the assumption of isotropic permeability. The history match with an analytical horizontal well model, shown in Fig. 15, confirms the results of the Horner and tandem-root plots.

Field example Well C

Well C data are from a buildup test of a horizontal well in a high-permeability sandstone where a 54-ft oil column overlies an extensive aquifer estimated to be approximately 180 ft thick. Table 8 shows the available data.

The diagnostic plot (Fig. 18) shows essentially no wellbore storage and a constant derivative, indicating radial, hemiradial, or elliptical flow at early times. The rapid decline of the derivative at the end of the test is caused by the aquifer underlying the oil column, which is acting like a constant-pressure boundary. A history match with an analytical horizontal well model with one no-flow and one constant-pressure boundary (the lower boundary), yielded kh = 313 md, kz = 7.5 md, sa = 1.5, and Lw = 356 ft. The no-flow boundary was estimated to be approximately 112 ft below the wellbore.

The time at which the early radial flow regime ends—the time where the derivative ceases to be flat on the diagnostic plot—is approximately 1.5 hours. For a wellbore with a volume of 130 bbl filled with a fluid of compressibility of 3.5 × 10–6 psi, the duration of wellbore storage (the unit-slope line) is estimated to be 0.0005 hours. With the gauge sampling rate set at 0.017 hours, the wellbore-storage unit slope simply could not be detected and does not appear at all on this plot. Fig. 19 is the Horner plot for this test. A straight line appears in the same range as the flat derivative on a diagnostic plot. From the slope of the line, the permeability is estimated to be 53 md, close to the regression analysis match value of 48 md.

Running horizontal well tests

The measurements in horizontal wells are usually made above the wellbore with the pressure gauge still in the vertical section.

The test string may often be too rigid to pass through the wellbore. However, in most cases, conventional hardware can be used for horizontal well tests. With longer horizontal wellbores, wellbore storage is an inherent problem for testing, even for buildup tests with downhole shut-in. As mentioned previously, problems arise in conducting buildup tests with short-duration production periods because superposition is inappropriate; therefore, Horner plots and tandem-root plots, which depend on superposition being applicable, are often inappropriate.

Another problem in conducting buildup test following a short production period is that significant pressure gradients along the length of the wellbore may cause crossflow within the wellbore during shut-in, so fluid may flow from one region to another in the wellbore. These gradients can be removed and this crossflow eliminated by a longer-duration flow period preceding a buildup test.

Factors that affect transient responses

A number of factors may affect the transient response of a horizontal well test: horizontal permeability (normal and parallel to well trajectory), vertical permeability, drilling damage, completion damage, producing interval that may be effectively much less than drilled length, and variations in standoff along length of well.

In summary, seven or more factors may affect interpretation for horizontal wells in homogeneous reservoirs before the effects from boundaries. The problem is complex, so test results are frequently inconclusive. Furthermore, wellbore storage inhibits determination of properties associated with early-time transient behavior such as vertical permeability and damage from drilling and completion. Middle- and late-time behavior may require several hours, days, or months to appear in transient data.

Some practical steps will help ensure interpretable test data. First, it is helpful to run tests in the pilot hole before kicking off to drill the horizontal borehole section. From a test in the vertical section, it is possible to get usable estimates of horizontal and vertical permeabilities using modern wireline test tools. Second, a good directional drilling survey can frequently provide an adequate estimate of standoff. A production-log flow survey conducted with coiled tubing can determine what part of the wellbore is actually producing and, therefore, help provide an estimate of effective productive length.

Wells in developed reservoirs should be flowed long enough to bring pressures along the wellbore to equilibrium and thus minimize crossflow. For high-rate wells, continuous borehole pressure and flow-rate measurements acquired during production can be used to interpret the pressure-drawdown transient response. If the downhole rates are not measured, the buildup test should be conducted with downhole shut-in to minimize wellbore storage distortion of test data.

Nomenclature

a = RTENOTITLE, stabilized deliverability coefficient, psia2-cp/MMscf-D
a = total length of reservoir perpendicular to wellbore, ft
ah = length of reservoir perpendicular to horizontal well, ft
af = RTENOTITLE, depth of investigation along major axis in fractured well, ft
at = RTENOTITLE, transient deliverability coefficient, psia2-cp/MMscf-D
aH = total width of reservoir perpendicular to the wellbore, ft
aH = modified total width of reservoir perpendicular to the wellbore, ft
A = drainage area, sq ft
A = πafbf , area of investigation in fractured well, ft2
Af = cross-sectional area perpendicular to flow, sq ft
Awb = wellbore area, sq ft
b = RTENOTITLE (gas flow equation)
bf = RTENOTITLE, depth of investigation of along minor axis in fractured well, ft
bB = intercept of Cartesian plot of bilinear flow data, psi
bH = length in direction parallel to wellbore, ft
bH = modified length in direction parallel to wellbore, ft
bL = intercept of Cartesian plot of linear flow data, psi
bV = intercept of Cartesian plot of data during volumetric behavior, psi
B = formation volume factor, res vol/surface vol
Bg = gas formation volume factor, RB/STB
Bgi = gas formation volume factor evaluated at pi, RB/Mscf
Bo = oil formation volume factor, RB/Mscf
Bw = water formation volume factor, RB/STB
RTENOTITLE = gas formation volume factor evaluated at average drainage area pressure, RB/Mscf
BND = 1,422 μ z TD/kh, non-Darcy flow coefficient
c = compressibility, psi–1
cf = formation compressibility, psi–1
cg = gas compressibility, psi–1
co = oil compressibility, psi–1
ct = Soco + Swcw + Sgcg + cf = total compressibility, psi–1
cw = water compressibility, psi–1
RTENOTITLE = total compressibility evaluated at average drainage area pressure, psi–1
ctf = total compressibility of pore space and fluids in fracture porosity, psi–1
ctm = total compressibility of pore space and fluids in matrix porosity, psi–1
cwb = compressibility of fluid in wellbore, psi–1
C = performance coefficient in gas-well deliverability equation, or wellbore storage coefficient, bbl/psi
CA = shape factor or constant
CD = 0.8936 C/ϕcthrw2 , dimensionless wellbore storage coefficient
(CDe2s)f = type-curve parameter value for the formation
(CDe2s)f +m = type-curve parameter value for the formation plus the matrix
CLfD = 0.8936 RTENOTITLE, dimensionless wellbore storage coefficient in fractured well
Cr = wfkf/πkLf, fracture conductivity, dimensionless
dx = shortest distance between horizontal well and x boundary, ft
dy = shortest distance between tip of horizontal well and y boundary, ft
dz = shortest distance between horizontal well and z boundary, ft
Dx = longest distance between horizontal well and x boundary, ft
Dy = longest distance between tip of horizontal well and y boundary, ft
Dz = longest distance between horizontal well and z boundary, ft
D = non-Darcy flow constant, D/Mscf
ebt = exponential decline with a constant b and elapsed time, t
Ef = flow efficiency, dimensionless
Ei(–x) = RTENOTITLE, the exponential integral
F(u) = function used in horizontal well analysis
FCD = wfkf/kLf, fracture conductivity, dimensionless
g = acceleration due to gravity, ft/sec2
gc = gravitational units conversion factor, 32.17 (lbm/ft)(lbf-s2)
h = net formation thickness, ft
hD = (h/rw)(kh/kv)1/2, dimensionless
hf = fracture height, ft
hm = thickness of matrix, ft
hp = perforated interval thickness, ft
hpD = hp/ht
ht = total formation thickness, ft
h1 = distance from top of formation to top of perforations, ft
h1D = h1/ht
HTRavg = HTR at average drainage area pressure
J = productivity index, STB/D, psi
Jactual = actual well productivity index, STB/D-psi
Jideal = ideal productivity index (s = 0), STB/D-psi
k = matrix permeability, md
RTENOTITLE = average permeability, md
kf = permeability of the proppant in the fracture, md
kfs = permeability near the wellbore, md
kg = permeability to gas, md
kgp = permeability of the gravel in the gravel pack, md
kh = horizontal permeability, md
km = matrix permeability, md
ko = permeability to oil, md
kr = permeability in horizontal radial direction, md
ks = permeability of altered zone, md
kw = permeability to water, md
kx = permeability in x-direction, md
ky = permeability in y-direction, md
kz = permeability in z-direction, md
L = distance from well to no-flow boundary, ft
Ld = drilled length of horizontal well, ft
Lf = fracture half length, ft
Lg = length of flow path through gravel pack, ft
Lm = length of matrix, ft
Lp = length of perforation tunnel, ft
Ls = length of damaged zone in fracture, ft
Lw = completed length of horizontal well, ft
Lx = distance from boundary, ft
m = 162.2 qBμ/kh = slope of middle-time line, psi/cycle
mB = slope of bilinear flow graph, psi/hr1/4
mL = slope of linear flow graph, psi/hr1/2
ms = RTENOTITLE, slope of spherical flow plot, psi-hr1/2
mV = slope of volumetric flow graph, psi/hr
mhrf = slope of semilog plot for hemiradial flow, psi/log cycle
melf = slope of square-root-of-time plot for early linear flow, psi/RTENOTITLE
merf = slope of semilog plot of early radial flow, psi/log cycle
mllf = slope of square-root-of-time plot for late linear flow, psi/RTENOTITLE
mprf = slope of semilog plot for pseudoradial flow, psi/log cycle
M = Molecular weight of gas
MTR = middle-time region
n = inverse slope of the line on a log-log plot of the change in pressure squared or pseudopressure vs. gas flow rate
p = pressure, psi
pavg = average pressure, psi
pb = base (atmospheric) pressure, psia
p0 = arbitrary reference or base pressure, psi
RTENOTITLE = volumetric average or static drainage-area pressure, psi
pa = adjusted or normalized pseudopressure, (μz/p)pp, psia
pawf = adjusted flowing bottomhole pressure, psia
paws = adjusted shut-in bottomhole pressure, psia
pf = formation pressure, psi
pi = original reservoir pressure, psi
pm = matrix pressure, psi
pp = pseudopressure, psia2/cp
ps = stabilized shut-in BHP measured just before start of a deliverability test, psia
psc = standard-condition pressure, psia
pt = surface pressure in tubing, psi
pw = BHP in wellbore, psi
pwf = flowing BHP, psi
pws = shut-in BHP, psi
pxy = parameter in horizontal well analysis equations
pxyz = parameter in horizontal well analysis equations
py = parameter in horizontal well analysis equations
p1hr = pressure at 1-hour shut-in (flow) time on MTR line or its extrapolation, psi
p = pressure derivative
p* = MTR pressure trend extrapolated to infinite shut-in time, psi
pD = 0.00708 kh(pip)/qBμ, dimensionless pressure as defined for constant-rate production
pMBHD = Matthews-Brons-Hazebroek pressure, dimensionless
(pD)MP = dimensionless pressure at match point
q = flow rate at surface, STB/D
qAOF = absolute-open-flow potential, MMscf/D
qg = gas flow rate, Mscf/D
qo = water flow rate, STB/D
qRt = total flow rate at reservoir conditions, RB/D
qsf = flow rate at formation (sand) face, STB/D
qw = water flow rate, STB/D
r = distance from the center of wellbore, ft
ra = radius of altered zone (skin effect), ft
rd = effective drainage radius, ft
rdp = radius of damage zone around perforation tunnel, ft
re = external drainage radius, ft
ri = radius of investigation, ft
rp = radius of perforation tunnel, ft
rs = outer radius of the altered zone, ft
rsp = radius of source or inner boundary of spherical flow pattern, ft
rw = wellbore radius, ft
rwa = apparent or effective wellbore radius, ft
rD = r/rw, dimensionless radius
Rs = dissolved GOR, scf/STB
s = skin factor, dimensionless
sa = skin caused by alteration of permeability around wellbore, dimensionless
sc = convergence skin, dimensionless
sd = skin caused by formation damage, dimensionless
se = skin caused by eccentric effects, dimensionless
sdp = perforation damage skin, dimensionless
sf = skin of hydraulically fractured well, dimensionless
sgp = skin factor from to Darcy flow through gravel pack, dimensionless
smin = minimum skin factor, dimensionless
sp = skin resulting from an incompletely perforated interval, dimensionless
st = total skin, dimensionless
sθ = skin factor resulting from well inclination, dimensionless
s = s + Dq = apparent skin factor, dimensionless
Sg = gas saturation, fraction of pore volume
So = oil saturation, fraction of pore volume
Sw = water saturation, fraction of pore volume
t = elapsed time, hours
ta = μcttap, adjusted or normalized pseudotime, hours
tap = pseudotime, hours
tbD = dimensionless time in linear flow, hours
tD = 0.0002637kt/ϕμctrw2, dimensionless time
tDA = 0.0002637 kt/ϕμctA = dimensionless time based on drainage area, A
teqB = equivalent time for bilinear flow, hours
te = equivalent time, hours
tLfD = 0.0002637 kt/ϕμetLf2, dimensionless time for fractured wells
tp = pseudoproducing time, hours
tpD = pseudoproducing time, dimensionless
tprf = time required to reach the pseudoradial flow regime, hours
tEelf = end of early linear flow, t, hours
tEerf = end of early radial flow, t, hours
tElf = end of linear flow, hours
tEllf = time to end of late linear flow regime, hours
tEhrf = end of hemiradial flow, hours
tErf = end of early radial flow, hours
tEprf = end of pseudoradial flow, hours
tp = constant-rate production period, t, hours
tpAD = dimensionless producing time, hours
tpss = time required to reach pseudosteady state, hours
tSelf = start of early linear flow, hours
tSllf = start of late linear flow, hours
tShrf = start of hemiradial flow, t, hours
tSprf = start of pseudoradial flow, t, hours
ts = time required for stabilization, hours
T = reservoir temperature, °R
Tsc = standard condition temperature, °R
u = dummy variable
V = volume, bbl
Vf = fraction of bulk volume occupied by fractures
Vm = fraction of bulk volume occupied by matrix
Vw = Vwb = wellbore volume, bbl
w = width of channel reservoir, ft
wf = fracture width, ft
wkf = fracture conductivity, md-ft
ws = width of damaged zone around fracture face, ft
WBS = wellbore storage
z = gas-law deviation factor, dimensionless
RTENOTITLE = gas-law deviation factor at average reservoir pressure, dimensionless
Δp = pressure change since start of transient test, psi
p)MP = pressure change at match point
ΔpD = dimensionless pressure change
Δpp = pseudopressure change since start of test, psia2/cp
Δps = additional pressure drop due to skin, psi
Δpt=0 = pressure drop at time zero, psi
Δp1hr = pressure change from start of test to one hour elapsed time, psi
Δt = time elapsed since start of test, hours
Δta = RTENOTITLE, normalized or adjusted pseudotime, hours
Δtap = RTENOTITLE, pseudotime, hr-psia/cp
ΔtBe = bilinear equivalent time, hours
Δte = radial equivalent time, hours
ΔtLe = linear equivalent time, hours
Δtmax = maximum shut-in time in pressure buildup test, hours
ΔV = change in volume, bbl
η = 0.0002637 k/ϕμct, hydraulic diffusivity, ft2/hr
ηfD = hydraulic diffusivity, dimensionless
λ = interporosity flow coefficient
λt = RTENOTITLE, total mobility, md/cp
α = exponent in deliverability equation
α = parameter characteristic of system geometry in dual-porosity system
β = turbulence factor
β = transition parameter
γ = Euler’s constant, = 1.781, dimensionless
γg = gas gravity (air = 1.0)
γm = matrix density
ω = storativity ratio in dual porosity reservoir
μ = viscosity, cp
μi = viscosity evaluated at pi, cp
μg = gas viscosity, cp
μo = oil viscosity, cp
μw = water viscosity, cp
RTENOTITLE = gas viscosity evaluated at average pressure, cp
μgwf = gas viscosity evaluated at pwf , cp
RTENOTITLE = viscosity evaluated at RTENOTITLE, cp
ρ = density, lbm/ft3 or g/cm3
ρwb = density of liquid in wellbore, lbm/ft3
ϕf = fraction of fracture volume occupied by pore space, ≅ 1
ϕm = fraction of matrix volume occupied by pore space
(ϕV)f = fraction of bulk volume occupied by pore space in fractures
(ϕVct)f = fracture "storativity" for dual porosity reservoir
(ϕVct)f+m = total "storativity" for dual porosity reservoir
(ϕV)m = fraction of bulk volume occupied by pore space in matrix
ϕ = porosity, dimensionless
Σs = sum of damage skin, turbulence, and other pseudoskin factors

References

  1. 1.0 1.1 1.2 1.3 1.4 Goode, P.A. and Thambynayagam, R.K.M. 1987. Pressure Drawdown and Buildup Analysis of Horizontal Wells in Anisotropic Media. SPE Form Eval 2 (4): 683–697. SPE-14250-PA. http://dx.doi.org/10.2118/14250-PA
  2. 2.00 2.01 2.02 2.03 2.04 2.05 2.06 2.07 2.08 2.09 2.10 2.11 Odeh, A.S. and Babu, D.K. 1990. Transient Flow Behavior of Horizontal Wells: Pressure Drawdown and Buildup Analysis. SPE Form Eval 5 (1): 7-15. SPE-18802-PA. http://dx.doi.org/10.2118/18802-PA
  3. Kuchuk, F.J. 1995. Well Testing and Interpretation for Horizontal Wells. J Pet Technol 47 (1): 36–41. SPE-25232-PA. http://dx.doi.org/10.2118/25232-PA
  4. Lichtenberger, G.J. 1994. Data Acquisition and Interpretation of Horizontal Well Pressure-Transient Tests. J Pet Technol 46 (2): 157-162. SPE-25922-PA. http://dx.doi.org/10.2118/25922-PA

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

Type curves

Estimating horizontal well productivity

Coning

Pressure transient testing

Horizontal wells

Fluid flow through permeable media

PEH:Fluid_Flow_Through_Permeable_Media