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 "Deliverability testing of gas wells"

Line 1: | Line 1: | ||

− | This article summarizes the fundamental gas-flow equations, both theoretical and empirical, used to analyze deliverability tests in terms of pseudopressure. The four most common types of gas-well deliverability tests are discussed in separate articles: [[Flow-after- | + | This article summarizes the fundamental gas-flow equations, both theoretical and empirical, used to analyze deliverability tests in terms of pseudopressure. The four most common types of gas-well deliverability tests are discussed in separate articles: [[Flow-after-flow_tests_for_gas_wells|flow-after-flow]], [[Single-point_tests_for_gas_wells|single-point]], [[Isochronal_tests_for_gas_wells|isochronal]], and [[Modified_isochronal_tests_for_gas_wells|modified isochronal]] tests. |

== Types and purposes of deliverability tests == | == Types and purposes of deliverability tests == | ||

− | |||

− | + | Deliverability testing refers to the testing of a gas well to measure its production capabilities under specific conditions of reservoir and bottomhole flowing pressures (BHFPs). A common productivity indicator obtained from these tests is the absolute open flow (AOF) potential. The AOF is the maximum rate at which a well could flow against a theoretical atmospheric backpressure at the sandface. Although in practice the well cannot produce at this rate, regulatory agencies sometimes use the AOF to allocate allowable production among wells or to set maximum production rates for individual wells. | |

− | + | Another application of deliverability testing is to generate a [[Reservoir_inflow_performance|reservoir inflow performance]] relationship (IPR) or gas backpressure curve. The IPR curve describes the relationship between surface production rate and BHFP for a specific value of reservoir pressure (that is, either the original pressure or the current average value). The IPR curve can be used to evaluate gas-well current deliverability potential under a variety of surface conditions, such as production against a fixed backpressure. In addition, the IPR can be used to forecast future production at any stage in the reservoir’s life. | |

− | [[ | + | Several deliverability testing methods have been developed for gas wells. [[Flow-after-flow_tests_for_gas_wells|Flow-after-flow tests]] are conducted by producing the well at a series of different stabilized flow rates and measuring the stabilized BHP. Each flow rate is established in succession without an intermediate shut-in period. A [[Single-point_tests_for_gas_wells|single-point test]] is conducted by flowing the well at a single rate until the BHFP is stabilized. This type of test was developed to overcome the limitation of long testing times required to reach stabilization at each rate in the flow-after-flow test. |

+ | |||

+ | [[Isochronal_tests_for_gas_wells|Isochronal]] and [[Modified_isochronal_tests_for_gas_wells|modified isochronal]] tests were developed to shorten tests times for wells that need long times to stabilize. An isochronal test consists of a series of single-point tests usually conducted by alternately producing at a slowly declining sandface rate without pressure stabilization and then shutting in and allowing the well to build to the average reservoir pressure before the next flow period. The modified isochronal test is conducted similarly, except the flow periods are of equal duration and the shut-in periods are of equal duration (but not necessarily the same as the flow periods). | ||

== Theory of deliverability test analysis == | == Theory of deliverability test analysis == | ||

− | |||

− | ===Theoretical deliverability equations=== | + | The theoretical equations developed by Houpeurt<ref name="r1">Houpeurt, A. 1959. On the Flow of Gases in Porous Media. Revue de L’lnstitut Francais du Petrole 14 (11): 1468.</ref> are exact solutions to the generalized radial-flow diffusivity equation, while the Rawlins and Schellhardt<ref name="r2">Rawlins, E.L. and Schellhardt, M.A. 1935. Backpressure Data on Natural Gas Wells and Their Application to Production Practices, Vol. 7. Monograph Series, USBM.</ref> equation was developed empirically. All basic equations presented here assume radial flow in a homogeneous, isotropic reservoir and therefore may not be applicable to the analysis of deliverability tests from reservoirs with heterogeneities, such as natural fractures or layered pay zones. These equations should not be used to analyze tests from hydraulically fractured wells during the fracture-dominated linear or bilinear flow periods. Finally, these equations assume that wellbore-storage effects have ceased. Unfortunately, wellbore-storage distortion may affect the entire test period in short tests, especially those conducted in low-permeability reservoirs. |

− | The early-time transient solution to the diffusivity equation for gases for constant-rate production from a well in a reservoir with closed outer boundaries, written in terms of pseudopressure, ''p''<sub>''p''</sub>,<ref name="r3" /> is | + | |

+ | === Theoretical deliverability equations === | ||

+ | |||

+ | The early-time transient solution to the diffusivity equation for gases for constant-rate production from a well in a reservoir with closed outer boundaries, written in terms of pseudopressure, ''p''<sub>''p''</sub>,<ref name="r3">Al-Hussainy, R., Jr., H.J.R., and Crawford, P.B. 1966. The Flow of Real Gases Through Porous Media. J Pet Technol 18 (5): 624-636. http://dx.doi.org/10.2118/1243-A-PA.</ref> is | ||

− | [[File:Vol5 page 0841 eq 001.png]]....................(1) | + | [[File:Vol5 page 0841 eq 001.png|RTENOTITLE]]....................(1) |

− | where ''p''<sub>''s''</sub> is the stabilized shut-in BHP measured before the deliverability test. In new reservoirs with little or no pressure depletion, this shut-in pressure equals the initial reservoir pressure, ''p''<sub>''s''</sub> = ''p''<sub>''i''</sub>, while in developed reservoirs, ''p''<sub>''s''</sub> < ''p''<sub>''i''</sub>. | + | where ''p''<sub>''s''</sub> is the stabilized shut-in BHP measured before the deliverability test. In new reservoirs with little or no pressure depletion, this shut-in pressure equals the initial reservoir pressure, ''p''<sub>''s''</sub> = ''p''<sub>''i''</sub>, while in developed reservoirs, ''p''<sub>''s''</sub> < ''p''<sub>''i''</sub>. |

The late-time or pseudosteady-state solution is | The late-time or pseudosteady-state solution is | ||

− | [[File:Vol5 page 0841 eq 002.png]]....................(2) | + | [[File:Vol5 page 0841 eq 002.png|RTENOTITLE]]....................(2) |

− | where [[File:Vol5 page 0781 inline 001.png]] is current drainage-area pressure. Gas wells cannot reach true pseudosteady state because ''μ''<sub>''g''</sub>(''p'')''c''<sub>''t''</sub>(''p'') changes as [[File:Vol5 page 0781 inline 001.png]] decreases. Note that, unlike [[File:Vol5 page 0781 inline 001.png]], which decreases during pseudosteady-state flow, ''p''<sub>''s''</sub> is a constant. | + | where [[File:Vol5 page 0781 inline 001.png|RTENOTITLE]] is current drainage-area pressure. Gas wells cannot reach true pseudosteady state because ''μ''<sub>''g''</sub>(''p'')''c''<sub>''t''</sub>(''p'') changes as [[File:Vol5 page 0781 inline 001.png|RTENOTITLE]] decreases. Note that, unlike [[File:Vol5 page 0781 inline 001.png|RTENOTITLE]], which decreases during pseudosteady-state flow, ''p''<sub>''s''</sub> is a constant. |

− | '''Eqs. 1''' and '''2''' are quadratic in terms of the gas flow rate, ''q''. For convenience, Houpeurt<ref name="r1" /> wrote the transient flow equation as | + | '''Eqs. 1''' and '''2''' are quadratic in terms of the gas flow rate, ''q''. For convenience, Houpeurt<ref name="r1">Houpeurt, A. 1959. On the Flow of Gases in Porous Media. Revue de L’lnstitut Francais du Petrole 14 (11): 1468.</ref> wrote the transient flow equation as |

− | [[File:Vol5 page 0841 eq 003.png]]....................(3) | + | [[File:Vol5 page 0841 eq 003.png|RTENOTITLE]]....................(3) |

and the pseudosteady-state flow equation as | and the pseudosteady-state flow equation as | ||

− | [[File:Vol5 page 0841 eq 004.png]]....................(4) | + | [[File:Vol5 page 0841 eq 004.png|RTENOTITLE]]....................(4) |

where | where | ||

− | [[File:Vol5 page 0841 eq 005.png]]....................(5) | + | [[File:Vol5 page 0841 eq 005.png|RTENOTITLE]]....................(5) |

− | [[File:Vol5 page 0841 eq 006.png]]....................(6) | + | [[File:Vol5 page 0841 eq 006.png|RTENOTITLE]]....................(6) |

− | [[File:Vol5 page 0841 eq 007.png]]....................(7) | + | [[File:Vol5 page 0841 eq 007.png|RTENOTITLE]]....................(7) |

− | The coefficients of ''q'' (''a''<sub>''t''</sub> for transient flow and a for pseudosteady-state flow) include the Darcy flow and skin effects and are measured in (psia<sup>2</sup>/cp)/(MMscf/D) when ''q'' is in MMscf/D. The coefficient of ''q''<sup>2</sup> represents the inertial and turbulent flow effects and is measured in (psia<sup>2</sup>/cp)/(MMscf/D)<sup>2</sup> when ''q'' is in MMscf/D. | + | The coefficients of ''q'' (''a''<sub>''t''</sub> for transient flow and a for pseudosteady-state flow) include the Darcy flow and skin effects and are measured in (psia<sup>2</sup>/cp)/(MMscf/D) when ''q'' is in MMscf/D. The coefficient of ''q''<sup>2</sup> represents the inertial and turbulent flow effects and is measured in (psia<sup>2</sup>/cp)/(MMscf/D)<sup>2</sup> when ''q'' is in MMscf/D. |

The Houpeurt equations also can be written in terms of pressure squared and are derived directly from the solutions to the gas-diffusivity equation, assuming that ''μ''<sub>''g''</sub>''z'' is constant over the pressure range considered. For transient flow, | The Houpeurt equations also can be written in terms of pressure squared and are derived directly from the solutions to the gas-diffusivity equation, assuming that ''μ''<sub>''g''</sub>''z'' is constant over the pressure range considered. For transient flow, | ||

− | [[File:Vol5 page 0841 eq 008.png]]....................(8) | + | [[File:Vol5 page 0841 eq 008.png|RTENOTITLE]]....................(8) |

and for pseudosteady-state flow, | and for pseudosteady-state flow, | ||

− | [[File:Vol5 page 0842 eq 001.png]]....................(9) | + | [[File:Vol5 page 0842 eq 001.png|RTENOTITLE]]....................(9) |

The flow coefficients are | The flow coefficients are | ||

− | [[File:Vol5 page 0842 eq 002.png]]....................(10) | + | [[File:Vol5 page 0842 eq 002.png|RTENOTITLE]]....................(10) |

− | [[File:Vol5 page 0842 eq 003.png]]....................(11) | + | [[File:Vol5 page 0842 eq 003.png|RTENOTITLE]]....................(11) |

− | and [[File:Vol5 page 0842 eq 004.png]]....................(12) | + | and [[File:Vol5 page 0842 eq 004.png|RTENOTITLE]]....................(12) |

− | When the Houpeurt equation is presented in terms of pressure squared, the coefficients of ''q'' are measured in psia<sup>2</sup>/(MMscf/D) when ''q'' is in MMscf/D, while the coefficient of ''q''<sup>2</sup> is measured in units of psia<sup>2</sup>/(MMscf/D)<sup>2</sup> when ''q'' is in MMscf/D. For convenience, all equations and examples in this section are presented with ''q'' measured in MMScf/D. | + | When the Houpeurt equation is presented in terms of pressure squared, the coefficients of ''q'' are measured in psia<sup>2</sup>/(MMscf/D) when ''q'' is in MMscf/D, while the coefficient of ''q''<sup>2</sup> is measured in units of psia<sup>2</sup>/(MMscf/D)<sup>2</sup> when ''q'' is in MMscf/D. For convenience, all equations and examples in this section are presented with ''q'' measured in MMScf/D. |

− | The pressure-squared form of the equation should be used only for gas reservoirs at low pressures (less than 2,000 psia) and high temperatures. To eliminate doubt about which equations to choose, use of the pseudopressure equations, which are applicable at all pressures and temperatures, is recommended. Consequently, all the analysis procedures in this section are presented in terms of pseudopressure. | + | The pressure-squared form of the equation should be used only for gas reservoirs at low pressures (less than 2,000 psia) and high temperatures. To eliminate doubt about which equations to choose, use of the pseudopressure equations, which are applicable at all pressures and temperatures, is recommended. Consequently, all the analysis procedures in this section are presented in terms of pseudopressure. |

− | An advantage of the pseudopressure form of the theoretical deliverability equation is that the flow coefficients are independent of the average reservoir pressure and, therefore, do not change as [[File:Vol5 page 0781 inline 001.png]] decreases during a flow test conducted under pseudosteady-state flow unless ''s'', ''k'', or ''A'' changes. Because the non-Darcy flow coefficient is a function of ''μ''<sub>''g''</sub>(''p''<sub>''wf''</sub> ), the coefficient b will change slightly if the BHFP is changed. In contrast, because of the pressure dependency of the gas properties on average reservoir pressure, the flow coefficients for the pressure-squared form of the deliverability equation must be recalculated for every new [[File:Vol5 page 0781 inline 001.png]] value. When ''s'', ''k'', or ''A'' changes with time, the only way to update the deliverability curve is to retest the well. | + | An advantage of the pseudopressure form of the theoretical deliverability equation is that the flow coefficients are independent of the average reservoir pressure and, therefore, do not change as [[File:Vol5 page 0781 inline 001.png|RTENOTITLE]] decreases during a flow test conducted under pseudosteady-state flow unless ''s'', ''k'', or ''A'' changes. Because the non-Darcy flow coefficient is a function of ''μ''<sub>''g''</sub>(''p''<sub>''wf''</sub> ), the coefficient b will change slightly if the BHFP is changed. In contrast, because of the pressure dependency of the gas properties on average reservoir pressure, the flow coefficients for the pressure-squared form of the deliverability equation must be recalculated for every new [[File:Vol5 page 0781 inline 001.png|RTENOTITLE]] value. When ''s'', ''k'', or ''A'' changes with time, the only way to update the deliverability curve is to retest the well. |

− | ===Empirical deliverability equations=== | + | === Empirical deliverability equations === |

− | |||

− | [[File:Vol5 page 0842 eq 005.png]]....................(13) | + | In 1935, Rawlins and Schellhardt<ref name="r2">Rawlins, E.L. and Schellhardt, M.A. 1935. Backpressure Data on Natural Gas Wells and Their Application to Production Practices, Vol. 7. Monograph Series, USBM.</ref> presented an empirical relationship that is used frequently in deliverability test analysis. The original form of their relation, given by '''Eq. 13''' in terms of pressure squared, is applicable only at low pressures: |

+ | |||

+ | [[File:Vol5 page 0842 eq 005.png|RTENOTITLE]]....................(13) | ||

In terms of pseudopressure, '''Eq. 13''' becomes | In terms of pseudopressure, '''Eq. 13''' becomes | ||

− | [[File:Vol5 page 0842 eq 006.png]]....................(14) | + | [[File:Vol5 page 0842 eq 006.png|RTENOTITLE]]....................(14) |

− | which is applicable over all pressure ranges. In '''Eqs. 13''' and '''14''', ''C'' is the stabilized performance coefficient and ''n'' is the inverse slope of the line on a log-log plot of the change in pressure squared or pseudopressure vs. gas flow rate. Depending on the flowing conditions, the theoretical value of ''n'' ranges from 0.5, indicating turbulent flow throughout a well’s drainage area, to 1.0, indicating laminar flow behavior modeled by Darcy’s law. The value of ''C'' changes depending on the units of flow rate and whether '''Eq. 13''' or '''14''' is used. All equations and examples in this section are presented with q measured in MMscf/D. | + | which is applicable over all pressure ranges. In '''Eqs. 13''' and '''14''', ''C'' is the stabilized performance coefficient and ''n'' is the inverse slope of the line on a log-log plot of the change in pressure squared or pseudopressure vs. gas flow rate. Depending on the flowing conditions, the theoretical value of ''n'' ranges from 0.5, indicating turbulent flow throughout a well’s drainage area, to 1.0, indicating laminar flow behavior modeled by Darcy’s law. The value of ''C'' changes depending on the units of flow rate and whether '''Eq. 13''' or '''14''' is used. All equations and examples in this section are presented with q measured in MMscf/D. |

− | Houpeurt proved that neither '''Eq. 13''' nor '''Eq. 14''' can be derived from the generalized diffusivity equation for radial flow of real gas through porous media. Although the Rawlins and Shellhardt equation is not theoretically rigorous, it is still widely used in deliverability analysis and has worked well over the years, especially when the test rates approach the AOF potential of the well and the extrapolation from test rates to AOF is minimal. | + | Houpeurt proved that neither '''Eq. 13''' nor '''Eq. 14''' can be derived from the generalized diffusivity equation for radial flow of real gas through porous media. Although the Rawlins and Shellhardt equation is not theoretically rigorous, it is still widely used in deliverability analysis and has worked well over the years, especially when the test rates approach the AOF potential of the well and the extrapolation from test rates to AOF is minimal. |

== Stabilization time == | == Stabilization time == | ||

− | |||

− | Stabilization time is defined as the time when the flowing pressure is no longer changing or is no longer changing significantly. Physically, stabilized flow can be interpreted as the time when the pressure transient is affected by the no-flow boundaries, either natural reservoir boundaries or an artificial boundaries created by active wells surrounding the tested well. Consider a graph of pressure as a function of radius for constant-rate flow at various times since the beginning of flow. As '''Fig. 1''' shows, the pressure in the wellbore continues to decrease as flow time increases. Simultaneously, the area from which fluid is drained increases, and the pressure transient moves farther out into the reservoir. | + | Unlike [[Pressure_transient_testing|pressure-transient tests]], the analysis techniques for conventional flow-after-flow and single-point tests require data obtained under stabilized flowing conditions. Although isochronal and modified isochronal tests were developed to circumvent the requirement of stabilized flow, they may still require a single, stabilized flow period at the end of the test. Consequently, there is a need to understand the meaning of stabilization time and have a method to estimate its value. |

+ | |||

+ | Stabilization time is defined as the time when the flowing pressure is no longer changing or is no longer changing significantly. Physically, stabilized flow can be interpreted as the time when the pressure transient is affected by the no-flow boundaries, either natural reservoir boundaries or an artificial boundaries created by active wells surrounding the tested well. Consider a graph of pressure as a function of radius for constant-rate flow at various times since the beginning of flow. As '''Fig. 1''' shows, the pressure in the wellbore continues to decrease as flow time increases. Simultaneously, the area from which fluid is drained increases, and the pressure transient moves farther out into the reservoir. | ||

<gallery widths="300px" heights="200px"> | <gallery widths="300px" heights="200px"> | ||

Line 88: | Line 93: | ||

</gallery> | </gallery> | ||

− | The [[ | + | The [[Fluid_flow_through_permeable_media#Radius_of_investigation_and_stabilization_time|radius of investigation]], the point in the formation beyond which the pressure drawdown is negligible, is a measure of how far a transient has moved into a formation following any rate change in a well. The approximate position of the radius of investigation at any time for a gas well is estimated by '''Eq. 15'''<ref name="r4">Lee, W.J. 1977. Well Testing, Vol. 1. Richardson, Texas: Textbook Series, SPE.</ref>: |

− | [[File:Vol5 page 0843 eq 001.png]]....................(15) | + | [[File:Vol5 page 0843 eq 001.png|RTENOTITLE]]....................(15) |

Stabilized flowing conditions occur when the calculated radius of investigation equals or exceeds the distance to the drainage boundaries of the well (i.e., ''r''<sub>''i''</sub> ≥ ''r''<sub>''e''</sub>). Substituting r e and rearranging '''Eq. 15''', yields an equation for estimating the stabilization time, ''t''<sub>''s''</sub>, for a gas well centered in a circular drainage area: | Stabilized flowing conditions occur when the calculated radius of investigation equals or exceeds the distance to the drainage boundaries of the well (i.e., ''r''<sub>''i''</sub> ≥ ''r''<sub>''e''</sub>). Substituting r e and rearranging '''Eq. 15''', yields an equation for estimating the stabilization time, ''t''<sub>''s''</sub>, for a gas well centered in a circular drainage area: | ||

− | [[File:Vol5 page 0843 eq 002.png]]....................(16) | + | [[File:Vol5 page 0843 eq 002.png|RTENOTITLE]]....................(16) |

− | As long as the radius of investigation is less than the distance to the no-flow boundary, stabilization has not been attained and the pressure behavior is transient. To illustrate the importance of stabilization times in deliverability testing, stabilization times were calculated as a function of permeability and drainage area for a well producing a gas with a specific gravity of 0.6 from a formation at 210°F and an average pressure of 3,500 psia [[File:Vol5 page 0843 inline 001.png]], with a porosity of 10%. Table 1 shows that, for wells completed in low-permeability reservoirs, several days—or even years—are required to reach stabilized flow, while wells completed in high-permeability reservoirs stabilize in a short time. | + | As long as the radius of investigation is less than the distance to the no-flow boundary, stabilization has not been attained and the pressure behavior is transient. To illustrate the importance of stabilization times in deliverability testing, stabilization times were calculated as a function of permeability and drainage area for a well producing a gas with a specific gravity of 0.6 from a formation at 210°F and an average pressure of 3,500 psia [[File:Vol5 page 0843 inline 001.png|RTENOTITLE]], with a porosity of 10%. Table 1 shows that, for wells completed in low-permeability reservoirs, several days—or even years—are required to reach stabilized flow, while wells completed in high-permeability reservoirs stabilize in a short time. |

− | <gallery widths=300px heights=200px> | + | <gallery widths="300px" heights="200px"> |

File:Vol5 Page 0844 Image 0001.png|'''Table 1''' | File:Vol5 Page 0844 Image 0001.png|'''Table 1''' | ||

</gallery> | </gallery> | ||

Line 104: | Line 109: | ||

A more general equation for calculating stabilization time is | A more general equation for calculating stabilization time is | ||

− | [[File:Vol5 page 0844 eq 001.png]]....................(17) | + | [[File:Vol5 page 0844 eq 001.png|RTENOTITLE]]....................(17) |

− | where ''t''<sub>''DA''</sub> is dimensionless time for the beginning of pseudosteady-state flow. Values for ''t''<sub>''DA''</sub> are given in '''Table 2''' for various reservoir shapes and well locations. <ref name="r5" /> The time required for the pseudosteady-state equation to be exact is found from the entry in the column "Exact for ''t''<sub>''DA''</sub> >." | + | where ''t''<sub>''DA''</sub> is dimensionless time for the beginning of pseudosteady-state flow. Values for ''t''<sub>''DA''</sub> are given in '''Table 2''' for various reservoir shapes and well locations. <ref name="r5">Earlougher, R.C. Jr. 1977. Advances in Well Test Analysis, Vol. 5. Richardson, Texas: Monograph Series, SPE.</ref> The time required for the pseudosteady-state equation to be exact is found from the entry in the column "Exact for ''t''<sub>''DA''</sub> >." |

− | <gallery widths=300px heights=200px> | + | <gallery widths="300px" heights="200px"> |

File:Vol5 Page 0889 Image 0001.png|'''Table 2''' | File:Vol5 Page 0889 Image 0001.png|'''Table 2''' | ||

Line 116: | Line 121: | ||

</gallery> | </gallery> | ||

− | The Rawlins-Schellhardt and Houpeurt deliverability equations assume radial flow. If pseudoradial flow has been achieved, however, these analysis techniques can be used for hydraulically fractured wells. The time to reach the pseudoradial flow regime, ''t''<sub>prf</sub>, occurs<ref name="r6" /> at [[File:Vol5 page 0844 inline 001.png]] and is estimated with | + | The Rawlins-Schellhardt and Houpeurt deliverability equations assume radial flow. If pseudoradial flow has been achieved, however, these analysis techniques can be used for hydraulically fractured wells. The time to reach the pseudoradial flow regime, ''t''<sub>prf</sub>, occurs<ref name="r6">Lee, W.J. 1989. Postfracture Formation Evaluation. In Recent Advances in Hydraulic Fracturing, J.L. Gidley, S.A. Holditch, D.E. Nierode, and R.W. Veatch Jr. eds., Vol. 12. Richardson, Texas: Monograph Series, SPE.</ref> at [[File:Vol5 page 0844 inline 001.png|RTENOTITLE]] and is estimated with |

− | [[File:Vol5 page 0844 eq 002.png]]....................(18) | + | [[File:Vol5 page 0844 eq 002.png|RTENOTITLE]]....................(18) |

− | To illustrate the importance of achieving pseudoradial flow during a deliverability test, values of t<sub>prf</sub> were calculated for a hydraulically fractured well completed in a reservoir with ''ϕ'' = 0.15, [[File:Vol5 page 0844 inline 002.png]] = 0.03 cp, and [[File:Vol5 page 0789 inline 002.png]] = 1 × 10<sup>−4</sup> psia<sup>−1</sup> and with the range of permeabilities and hydraulic fracture half-lengths in '''Table 3'''. The results illustrate that a well with a long fracture in a low-permeability formation will take far too long to stabilize for conventional deliverability testing. | + | To illustrate the importance of achieving pseudoradial flow during a deliverability test, values of t<sub>prf</sub> were calculated for a hydraulically fractured well completed in a reservoir with ''ϕ'' = 0.15, [[File:Vol5 page 0844 inline 002.png|RTENOTITLE]] = 0.03 cp, and [[File:Vol5 page 0789 inline 002.png|RTENOTITLE]] = 1 × 10<sup>−4</sup> psia<sup>−1</sup> and with the range of permeabilities and hydraulic fracture half-lengths in '''Table 3'''. The results illustrate that a well with a long fracture in a low-permeability formation will take far too long to stabilize for conventional deliverability testing. |

− | <gallery widths=300px heights=200px> | + | <gallery widths="300px" heights="200px"> |

File:Vol5 Page 0844 Image 0002.png|'''Table 3''' | File:Vol5 Page 0844 Image 0002.png|'''Table 3''' | ||

</gallery> | </gallery> | ||

== Nomenclature == | == Nomenclature == | ||

+ | |||

{| | {| | ||

− | |||

− | |||

− | |||

|- | |- | ||

− | |''a''<sub>''f''</sub> | + | | ''a'' |

− | |= | + | | = |

− | |[[File:Vol5 page 0879 inline 002.png]], depth of investigation along major axis in fractured well, ft | + | | [[File:Vol5 page 0879 inline 001.png|RTENOTITLE]], stabilized deliverability coefficient, psia<sup>2</sup>-cp/MMscf-D |

+ | |- | ||

+ | | ''a''<sub>''f''</sub> | ||

+ | | = | ||

+ | | [[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'' | + | | ''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> |

|- | |- | ||

− | |''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 |

|- | |- | ||

− | |[[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> |

|- | |- | ||

− | |''D'' | + | | ''D'' |

− | |= | + | | = |

− | |non-Darcy flow constant, D/Mscf | + | | non-Darcy flow constant, D/Mscf |

|- | |- | ||

− | |''h'' | + | | ''h'' |

− | |= | + | | = |

− | |net formation thickness, ft | + | | net formation thickness, ft |

|- | |- | ||

− | |''k''<sub>''g''</sub> | + | | ''k''<sub>''g''</sub> |

− | |= | + | | = |

− | |permeability to gas, md | + | | permeability to gas, md |

|- | |- | ||

− | |''L''<sub>''f''</sub> | + | | ''L''<sub>''f''</sub> |

− | |= | + | | = |

− | |fracture half length, ft | + | | fracture half length, ft |

|- | |- | ||

− | |''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 |

|- | |- | ||

− | |[[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>''p''</sub> | + | | ''p''<sub>''p''</sub> |

− | |= | + | | = |

− | |pseudopressure, psia<sup>2</sup>/cp | + | | pseudopressure, psia<sup>2</sup>/cp |

|- | |- | ||

− | |''p''<sub>''wf''</sub> | + | | ''p''<sub>''wf''</sub> |

− | |= | + | | = |

− | |flowing BHP, psi | + | | flowing BHP, psi |

|- | |- | ||

− | |''q'' | + | | ''q'' |

− | |= | + | | = |

− | |flow rate at surface, STB/D | + | | flow rate at surface, STB/D |

|- | |- | ||

− | |''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>''s''</sub> | + | | ''t''<sub>''s''</sub> |

− | |= | + | | = |

− | |time required for stabilization, hours | + | | time required for stabilization, hours |

|- | |- | ||

− | |''T'' | + | | ''T'' |

− | |= | + | | = |

− | |reservoir temperature, °R | + | | reservoir temperature, °R |

|- | |- | ||

− | |Δ''p'' | + | | Δ''p'' |

− | |= | + | | = |

− | |pressure change since start of transient test, psi | + | | pressure change since start of transient test, psi |

|- | |- | ||

− | |[[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 |

|- | |- | ||

− | |''ϕ'' | + | | ''ϕ'' |

− | |= | + | | = |

− | |porosity, dimensionless | + | | porosity, dimensionless |

|} | |} | ||

− | ==References== | + | == References == |

− | |||

− | |||

− | < | + | <references /> |

− | + | == Noteworthy papers in OnePetro == | |

− | + | Use this section to list papers in OnePetro that a reader who wants to learn more should definitely read | |

− | + | == External links == | |

− | + | Use this section to provide links to relevant material on websites other than PetroWiki and OnePetro | |

− | |||

− | == | + | == See also == |

− | |||

− | + | [[Flow_equations_for_gas_and_multiphase_flow|Flow equations for gas and multiphase flow]] | |

− | |||

− | + | [[Fluid_flow_through_permeable_media|Fluid flow through permeable media]] | |

− | [[ | ||

− | [[ | + | [[PEH:Fluid_Flow_Through_Permeable_Media]] |

− | + | ==Category== | |

− | [[Category:3.3.3 Downhole and wellsite flow metering]] | + | [[Category:3.3.3 Downhole and wellsite flow metering]] [[Category:YR]] |

## Latest revision as of 15:13, 2 July 2015

This article summarizes the fundamental gas-flow equations, both theoretical and empirical, used to analyze deliverability tests in terms of pseudopressure. The four most common types of gas-well deliverability tests are discussed in separate articles: flow-after-flow, single-point, isochronal, and modified isochronal tests.

## Types and purposes of deliverability tests

Deliverability testing refers to the testing of a gas well to measure its production capabilities under specific conditions of reservoir and bottomhole flowing pressures (BHFPs). A common productivity indicator obtained from these tests is the absolute open flow (AOF) potential. The AOF is the maximum rate at which a well could flow against a theoretical atmospheric backpressure at the sandface. Although in practice the well cannot produce at this rate, regulatory agencies sometimes use the AOF to allocate allowable production among wells or to set maximum production rates for individual wells.

Another application of deliverability testing is to generate a reservoir inflow performance relationship (IPR) or gas backpressure curve. The IPR curve describes the relationship between surface production rate and BHFP for a specific value of reservoir pressure (that is, either the original pressure or the current average value). The IPR curve can be used to evaluate gas-well current deliverability potential under a variety of surface conditions, such as production against a fixed backpressure. In addition, the IPR can be used to forecast future production at any stage in the reservoir’s life.

Several deliverability testing methods have been developed for gas wells. Flow-after-flow tests are conducted by producing the well at a series of different stabilized flow rates and measuring the stabilized BHP. Each flow rate is established in succession without an intermediate shut-in period. A single-point test is conducted by flowing the well at a single rate until the BHFP is stabilized. This type of test was developed to overcome the limitation of long testing times required to reach stabilization at each rate in the flow-after-flow test.

Isochronal and modified isochronal tests were developed to shorten tests times for wells that need long times to stabilize. An isochronal test consists of a series of single-point tests usually conducted by alternately producing at a slowly declining sandface rate without pressure stabilization and then shutting in and allowing the well to build to the average reservoir pressure before the next flow period. The modified isochronal test is conducted similarly, except the flow periods are of equal duration and the shut-in periods are of equal duration (but not necessarily the same as the flow periods).

## Theory of deliverability test analysis

The theoretical equations developed by Houpeurt^{[1]} are exact solutions to the generalized radial-flow diffusivity equation, while the Rawlins and Schellhardt^{[2]} equation was developed empirically. All basic equations presented here assume radial flow in a homogeneous, isotropic reservoir and therefore may not be applicable to the analysis of deliverability tests from reservoirs with heterogeneities, such as natural fractures or layered pay zones. These equations should not be used to analyze tests from hydraulically fractured wells during the fracture-dominated linear or bilinear flow periods. Finally, these equations assume that wellbore-storage effects have ceased. Unfortunately, wellbore-storage distortion may affect the entire test period in short tests, especially those conducted in low-permeability reservoirs.

### Theoretical deliverability equations

The early-time transient solution to the diffusivity equation for gases for constant-rate production from a well in a reservoir with closed outer boundaries, written in terms of pseudopressure, *p*_{p},^{[3]} is

where *p*_{s} is the stabilized shut-in BHP measured before the deliverability test. In new reservoirs with little or no pressure depletion, this shut-in pressure equals the initial reservoir pressure, *p*_{s} = *p*_{i}, while in developed reservoirs, *p*_{s} < *p*_{i}.

The late-time or pseudosteady-state solution is

where is current drainage-area pressure. Gas wells cannot reach true pseudosteady state because *μ*_{g}(*p*)*c*_{t}(*p*) changes as decreases. Note that, unlike , which decreases during pseudosteady-state flow, *p*_{s} is a constant.

**Eqs. 1** and **2** are quadratic in terms of the gas flow rate, *q*. For convenience, Houpeurt^{[1]} wrote the transient flow equation as

and the pseudosteady-state flow equation as

where

The coefficients of *q* (*a*_{t} for transient flow and a for pseudosteady-state flow) include the Darcy flow and skin effects and are measured in (psia^{2}/cp)/(MMscf/D) when *q* is in MMscf/D. The coefficient of *q*^{2} represents the inertial and turbulent flow effects and is measured in (psia^{2}/cp)/(MMscf/D)^{2} when *q* is in MMscf/D.

The Houpeurt equations also can be written in terms of pressure squared and are derived directly from the solutions to the gas-diffusivity equation, assuming that *μ*_{g}*z* is constant over the pressure range considered. For transient flow,

and for pseudosteady-state flow,

The flow coefficients are

When the Houpeurt equation is presented in terms of pressure squared, the coefficients of *q* are measured in psia^{2}/(MMscf/D) when *q* is in MMscf/D, while the coefficient of *q*^{2} is measured in units of psia^{2}/(MMscf/D)^{2} when *q* is in MMscf/D. For convenience, all equations and examples in this section are presented with *q* measured in MMScf/D.

The pressure-squared form of the equation should be used only for gas reservoirs at low pressures (less than 2,000 psia) and high temperatures. To eliminate doubt about which equations to choose, use of the pseudopressure equations, which are applicable at all pressures and temperatures, is recommended. Consequently, all the analysis procedures in this section are presented in terms of pseudopressure.

An advantage of the pseudopressure form of the theoretical deliverability equation is that the flow coefficients are independent of the average reservoir pressure and, therefore, do not change as decreases during a flow test conducted under pseudosteady-state flow unless *s*, *k*, or *A* changes. Because the non-Darcy flow coefficient is a function of *μ*_{g}(*p*_{wf} ), the coefficient b will change slightly if the BHFP is changed. In contrast, because of the pressure dependency of the gas properties on average reservoir pressure, the flow coefficients for the pressure-squared form of the deliverability equation must be recalculated for every new value. When *s*, *k*, or *A* changes with time, the only way to update the deliverability curve is to retest the well.

### Empirical deliverability equations

In 1935, Rawlins and Schellhardt^{[2]} presented an empirical relationship that is used frequently in deliverability test analysis. The original form of their relation, given by **Eq. 13** in terms of pressure squared, is applicable only at low pressures:

In terms of pseudopressure, **Eq. 13** becomes

which is applicable over all pressure ranges. In **Eqs. 13** and **14**, *C* is the stabilized performance coefficient and *n* is the inverse slope of the line on a log-log plot of the change in pressure squared or pseudopressure vs. gas flow rate. Depending on the flowing conditions, the theoretical value of *n* ranges from 0.5, indicating turbulent flow throughout a well’s drainage area, to 1.0, indicating laminar flow behavior modeled by Darcy’s law. The value of *C* changes depending on the units of flow rate and whether **Eq. 13** or **14** is used. All equations and examples in this section are presented with q measured in MMscf/D.

Houpeurt proved that neither **Eq. 13** nor **Eq. 14** can be derived from the generalized diffusivity equation for radial flow of real gas through porous media. Although the Rawlins and Shellhardt equation is not theoretically rigorous, it is still widely used in deliverability analysis and has worked well over the years, especially when the test rates approach the AOF potential of the well and the extrapolation from test rates to AOF is minimal.

## Stabilization time

Unlike pressure-transient tests, the analysis techniques for conventional flow-after-flow and single-point tests require data obtained under stabilized flowing conditions. Although isochronal and modified isochronal tests were developed to circumvent the requirement of stabilized flow, they may still require a single, stabilized flow period at the end of the test. Consequently, there is a need to understand the meaning of stabilization time and have a method to estimate its value.

Stabilization time is defined as the time when the flowing pressure is no longer changing or is no longer changing significantly. Physically, stabilized flow can be interpreted as the time when the pressure transient is affected by the no-flow boundaries, either natural reservoir boundaries or an artificial boundaries created by active wells surrounding the tested well. Consider a graph of pressure as a function of radius for constant-rate flow at various times since the beginning of flow. As **Fig. 1** shows, the pressure in the wellbore continues to decrease as flow time increases. Simultaneously, the area from which fluid is drained increases, and the pressure transient moves farther out into the reservoir.

The radius of investigation, the point in the formation beyond which the pressure drawdown is negligible, is a measure of how far a transient has moved into a formation following any rate change in a well. The approximate position of the radius of investigation at any time for a gas well is estimated by **Eq. 15**^{[4]}:

Stabilized flowing conditions occur when the calculated radius of investigation equals or exceeds the distance to the drainage boundaries of the well (i.e., *r*_{i} ≥ *r*_{e}). Substituting r e and rearranging **Eq. 15**, yields an equation for estimating the stabilization time, *t*_{s}, for a gas well centered in a circular drainage area:

As long as the radius of investigation is less than the distance to the no-flow boundary, stabilization has not been attained and the pressure behavior is transient. To illustrate the importance of stabilization times in deliverability testing, stabilization times were calculated as a function of permeability and drainage area for a well producing a gas with a specific gravity of 0.6 from a formation at 210°F and an average pressure of 3,500 psia , with a porosity of 10%. Table 1 shows that, for wells completed in low-permeability reservoirs, several days—or even years—are required to reach stabilized flow, while wells completed in high-permeability reservoirs stabilize in a short time.

A more general equation for calculating stabilization time is

where *t*_{DA} is dimensionless time for the beginning of pseudosteady-state flow. Values for *t*_{DA} are given in **Table 2** for various reservoir shapes and well locations. ^{[5]} The time required for the pseudosteady-state equation to be exact is found from the entry in the column "Exact for *t*_{DA} >."

The Rawlins-Schellhardt and Houpeurt deliverability equations assume radial flow. If pseudoradial flow has been achieved, however, these analysis techniques can be used for hydraulically fractured wells. The time to reach the pseudoradial flow regime, *t*_{prf}, occurs^{[6]} at and is estimated with

To illustrate the importance of achieving pseudoradial flow during a deliverability test, values of t_{prf} were calculated for a hydraulically fractured well completed in a reservoir with *ϕ* = 0.15, = 0.03 cp, and = 1 × 10^{−4} psia^{−1} and with the range of permeabilities and hydraulic fracture half-lengths in **Table 3**. The results illustrate that a well with a long fracture in a low-permeability formation will take far too long to stabilize for conventional deliverability testing.

## Nomenclature

## References

- ↑
^{1.0}^{1.1}Houpeurt, A. 1959. On the Flow of Gases in Porous Media. Revue de L’lnstitut Francais du Petrole 14 (11): 1468. - ↑
^{2.0}^{2.1}Rawlins, E.L. and Schellhardt, M.A. 1935. Backpressure Data on Natural Gas Wells and Their Application to Production Practices, Vol. 7. Monograph Series, USBM. - ↑ Al-Hussainy, R., Jr., H.J.R., and Crawford, P.B. 1966. The Flow of Real Gases Through Porous Media. J Pet Technol 18 (5): 624-636. http://dx.doi.org/10.2118/1243-A-PA.
- ↑ Lee, W.J. 1977. Well Testing, Vol. 1. Richardson, Texas: Textbook Series, SPE.
- ↑ Earlougher, R.C. Jr. 1977. Advances in Well Test Analysis, Vol. 5. Richardson, Texas: Monograph Series, SPE.
- ↑ Lee, W.J. 1989. Postfracture Formation Evaluation. In Recent Advances in Hydraulic Fracturing, J.L. Gidley, S.A. Holditch, D.E. Nierode, and R.W. Veatch Jr. eds., Vol. 12. Richardson, Texas: Monograph Series, SPE.

## 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

Flow equations for gas and multiphase flow

Fluid flow through permeable media

PEH:Fluid_Flow_Through_Permeable_Media