Difference between revisions of "Vogel's IPR"

From wiki.pengtools.com
Jump to: navigation, search
(1. Calculate the Productivity Index, J:)
 
(11 intermediate revisions by the same user not shown)
Line 1: Line 1:
 
__TOC__
 
__TOC__
 +
<div style='text-align: right;'>By Mikhail Tuzovskiy on {{REVISIONTIMESTAMP}}</div>
 
==Vogel's Inflow Performance Relationship==
 
==Vogel's Inflow Performance Relationship==
 
[[File:Vogel's.png|thumb|right|300px| Vogel's IPR<ref name=Vogel />]]
 
[[File:Vogel's.png|thumb|right|300px| Vogel's IPR<ref name=Vogel />]]
Line 42: Line 43:
 
:: <math>J = \frac{kh\ J_D}{141.2 B \mu}</math>
 
:: <math>J = \frac{kh\ J_D}{141.2 B \mu}</math>
 
:1.3 J from kh  and skin:
 
:1.3 J from kh  and skin:
:: <math>J = \frac{kh\ \frac{1}{1 / 0.13 + s}}{141.2 B \mu}</math>
+
:: <math>J = \frac{kh\ \frac{1}{1 / 0.13 + S}}{141.2 B \mu}</math>
  
 
===2. Calculate the flowing rates, q<sub>o</sub>:===
 
===2. Calculate the flowing rates, q<sub>o</sub>:===
Line 78: Line 79:
 
== Nomenclature  ==
 
== Nomenclature  ==
 
:<math> B </math> = formation volume factor, bbl/stb
 
:<math> B </math> = formation volume factor, bbl/stb
 +
:<math> J </math> = productivity index, stb/d/psia
 
:<math> J_D </math> = dimensionless productivity index, dimensionless
 
:<math> J_D </math> = dimensionless productivity index, dimensionless
 
:<math> kh</math> = permeability times thickness, md*ft
 
:<math> kh</math> = permeability times thickness, md*ft
 +
:<math> P </math> = pressure, psia
 
:<math> \bar{P} </math> = average reservoir pressure, psia
 
:<math> \bar{P} </math> = average reservoir pressure, psia
:<math> P_{\bar{P}} </math> = average reservoir pseudopressure, psia<sup>2</sup>/cP
 
:<math> P_{wf} </math> = well flowing pressure, psia
 
:<math> P_{P_{wf}} </math> = average well flowing pseudopressure, psia<sup>2</sup>/cP
 
 
:<math> q </math> = flowing rate, stb/d
 
:<math> q </math> = flowing rate, stb/d
:<math> q_g </math> = gas rate, MMscfd
+
:<math> S </math> = skin factor, dimensionless
:<math> T </math> = temperature, °R
 
  
 
===Greek symbols===
 
===Greek symbols===
Line 94: Line 93:
 
===Subscripts===
 
===Subscripts===
  
:g = gas<BR/>
+
:b = at bubble point pressure<BR/>
:K = °K<BR/>
+
:max = maximum<BR/>
:L = liquid<BR/>
 
:m = gas/liquid mixture<BR/>
 
 
:o = oil<BR/>
 
:o = oil<BR/>
:R = °R<BR/>
+
:test = well test<BR/>
:SL = superficial liquid<BR/>
+
:wf = well flowing bottomhole pressure<BR/>
:SG = superficial gas<BR/>
 
:w = water<BR/>
 
  
=== References ===
+
== References ==
 
<references>
 
<references>
 
<ref name=Vogel>{{cite journal
 
<ref name=Vogel>{{cite journal
Line 136: Line 131:
 
==See also==
 
==See also==
 
:[[IPR]]<BR/>
 
:[[IPR]]<BR/>
:[[141.2 derivation]]<BR/>
+
:[[3 Phase IPR]]<BR/>
 
:[[Darcy's law]]<BR/>
 
:[[Darcy's law]]<BR/>
 
:[[JD]]<BR/>
 
:[[JD]]<BR/>
:[[Production Potential]]
+
:[[141.2 derivation]]<BR/>
  
 
{{#seo:
 
{{#seo:

Latest revision as of 05:29, 3 January 2023

By Mikhail Tuzovskiy on 20230103052908

Vogel's Inflow Performance Relationship

Vogel's IPR[1]

Vogel's IPR is an empirical two-phase (oil + gas) inflow performance relationship correlation published in 1968 [1].

Vogel's IPR is based on computer simulations to several solution gas drive reservoirs for different fluid and reservoir relative permeability properties.

Vogel's IPR is the default IPR correlation to calculate oil wells performance in the PQplot nodal analysis software which is available online at petroleum engineering site www.pengtools.com.

Math and Physics

Vogel's IPR equation

 \frac{q_o}{q_{o_{max}}} = 1-0.2 \frac{P_{wf}}{\bar{P}} - 0.8 \left ( \frac{P_{wf}}{\bar{P}} \right )^2 [1]

Darcy's law + Vogel's IPR

Combination single phase liquid constant PI equation and Vogel's IPR:

Combination Constant PI and Vogel's IPR[2]

 q_{ob} = J (\bar{P} - P_b) [2] , oil flow rate at the bubble point.
 q_{o_{max}} = q_{ob} + \frac{J P_b}{1.8} [3] , maximum oil rate or absolute open flow (AOF).
 q_o = q_{ob} + (q_{o_{max}} - q_{ob})  \left (1-0.2 \frac{P_{wf}}{\bar{P}} - 0.8 \left ( \frac{P_{wf}}{\bar{P}} \right )^2 \right )[2] , oil rate at given flowing bottomhole pressure.
 J = \frac{q_o}{\bar{P}-P_b + \frac{P_b}{1.8} \left (1-0.2 \frac{P_{wf}}{\bar{P}} - 0.8 \left ( \frac{P_{wf}}{\bar{P}} \right )^2 \right ) } [2] , productivity index for test below the bubble point pressure.

Discussion

Why Vogel's IPR?

Vogel's IPR solution has been found to be very good and is widely used in prediction of IPR curves.
— Kermit Brown et al[2]

Vogel's IPR calculation workflow

1. Calculate the Productivity Index, J:

1.1 J from the flow test:
  • Test is above the bubble point:
J=\frac{q_{o_{test}}}{\bar{P}-P_{wf}}
  • Test is below the bubble point:
 J = \frac{q_{o_{test}}}{\bar{P}-P_b + \frac{P_b}{1.8} \left (1-0.2 \frac{P_{wf}}{\bar{P}} - 0.8 \left ( \frac{P_{wf}}{\bar{P}} \right )^2 \right ) }
1.2 J from kh and JD:
J = \frac{kh\ J_D}{141.2 B \mu}
1.3 J from kh and skin:
J = \frac{kh\ \frac{1}{1 / 0.13 + S}}{141.2 B \mu}

2. Calculate the flowing rates, qo:

For each Pwf from \bar{P} to 0:

2.1 Calculate oil flow rate at the bubble point:
 q_{ob} = J (\bar{P} - P_b)
2.2 Calculate maximum oil rate:
 q_{o_{max}} = q_{ob} + \frac{J P_b}{1.8}
2.3 Calculate oil rate at given flowing bottomhole pressure, Pwf:
 q_o = q_{ob} + (q_{o_{max}} - q_{ob})  \left (1-0.2 \frac{P_{wf}}{\bar{P}} - 0.8 \left ( \frac{P_{wf}}{\bar{P}} \right )^2 \right )

Vogel's IPR calculation example

Following the example problem #6, page 15 [2]:

Given:

\bar{P}=4200 psi
J=2 b/d/psi
Pb=3000 psi

Calculate:

(1) qob
(2) qmax
(3) q for Pwf=1500 psi

Solution:

(1)  q_{ob} = J (\bar{P} - P_b) = 2 \times (4200 - 3000) = 2400\ b/d
(2)  q_{o_{max}} = q_{ob} + \frac{J P_b}{1.8} = 2400 + \frac{2 \times 3000}{1.8} = 2400 + 3333 = 5733\ b/d
(3 ) q_o = q_{ob} + (q_{o_{max}} - q_{ob})  \left (1-0.2 \frac{P_{wf}}{\bar{P}} - 0.8 \left ( \frac{P_{wf}}{\bar{P}} \right )^2 \right )
 q_o = 2400 + (5733 - 2400)  \left (1-0.2 \frac{1500}{4200} - 0.8 \left ( \frac{1500}{4200} \right )^2 \right )
 q_o = 2400 + 3333 \times 0.7 = 4733\ b/d


With the help of the PQplot software, other values of flowing pressure has been assumed, and the corresponding values of flow rate were determined. From these, the IPR cure was plotted (Fig.1).

The PQplot model from this example is available online by the following link: Vogel's IPR calculation example

Nomenclature

 B = formation volume factor, bbl/stb
 J = productivity index, stb/d/psia
 J_D = dimensionless productivity index, dimensionless
 kh = permeability times thickness, md*ft
 P = pressure, psia
 \bar{P} = average reservoir pressure, psia
 q = flowing rate, stb/d
 S = skin factor, dimensionless

Greek symbols

 \mu = viscosity, cp

Subscripts

b = at bubble point pressure
max = maximum
o = oil
test = well test
wf = well flowing bottomhole pressure

References

  1. 1.0 1.1 1.2 Vogel, J. V. (1968). "Inflow Performance Relationships for Solution-Gas Drive Wells". Journal of Petroleum Technology. 20 (SPE-1476-PA). 
  2. 2.0 2.1 2.2 2.3 2.4 2.5 Brown, Kermit (1984). The Technology of Artificial Lift Methods. Volume 4. Production Optimization of Oil and Gas Wells by Nodal System Analysis. Tulsa, Oklahoma: PennWellBookss. 
  3. Neely, A.B. (1976). Use of IPR Curves. Houston, Texas: Shell Oil Co. 

See also

IPR
3 Phase IPR
Darcy's law
JD
141.2 derivation