CU PET 323 - Reservoir Engineering I - Full #118
CU PET 323 - Reservoir Engineering I - Full #118
• Petroleum Reservoir Fluid Property Correlations, McCain. Jr., Spivey and Lenn
• Outline:
– Essence of Reservoir Engineering
– Roles of Reservoir Engineering (in a multi-disciplinary field
development team)
– Specific Responsibilities of Reservoir Engineers
• Stages in Reservoir Engineering Operations
– Methods in Reservoir Engineering
– Physical Principles underlying Reservoir Engineering
(Don‟t tell other guys I told you! They will be green with envy)
Formation properties:
porosity, permeability,
payzone thickness, fluid
saturations.
Reservoir Engineering
Constructing models to predict the performance of
the field under a variety of development scenarios,
in conjunction with Production/Process
Engineering
oil/gas/water
of
Project Engineering
Constructing platforms
(offshore), surface
production/injection and
process facilities
PART 1: INTRODUCTION TO RESERVOIR
ENGINEERING
Specific Responsibilities of Reservoir Engineers
• Estimation of amount of hydrocarbon-in-place, in conjunction with
geologists and petrophysicist.
• Conservation of mass:
Rate of flow IN – Rate of flow OUT = Rate of Accumulation.
It is the principle underlying the MBE, a prominent analytical tool in RE
• Ideally, lab procedures are the most accurate way to determine values of
properties of a particular fluid.
• At the end of this chapter, the class shall be divided into groups; with
each group assigned to develop spreadsheets for given correlations.
𝒏𝑹𝑻
For ideal gas, 𝑷𝑽 = 𝒏𝑹𝑻, 𝑽 = i.e.
𝑷
𝑽 𝟏 𝒏𝑹𝑻
Dividing through by mass; = =
𝒎 𝝆 𝒎𝑷
𝒏 𝟏 𝟏
But = =
𝒎 𝑴𝒐𝒍𝒂𝒓 𝒘𝒆𝒊𝒈𝒉𝒕 𝑴
𝑷𝑴
Therefore, 𝝆 = … … 𝟐. 𝟏𝟎
𝑹𝑻
With gas composition known, the Apparent Molecular weight, Mg of a gas mixture can be
calculated using the simple mixing rule below:
𝒏𝒊
𝑴𝒈 = 𝒚𝒊 𝑴𝒊 … … 𝟐. 𝟏𝟐
𝒊<𝟏
yi is the mole fraction of the ith component and Mi is the molecular weight of the ith
component
𝑻
𝑻𝒑𝒓 = … … 𝟐. 𝟖
𝑻𝒑𝒄
TTOWG!
PART 2: PETROLEUM RESERVOIR FLUID
PROPERTIES
Calculation of Pseudo-critical Gas Properties
As discussed earlier, there are two scenarios: gas composition known;
gas gravity, γg known (composition unknown). The two scenarios
are handled separately by the following correlations:
• Gas gravity and Gas Composition known; use mixing rule.
• Gas Gravity, γg known; use Sutton‟s correlation
Brown‟s chart
• At 3000psia and 180
deg F, find the z factor
of the mixture whose
composition is shown
in the table to the right:
• Also find the critical P
and T using Standing‟s
correlation
𝑻′𝒑𝒄 = 𝑻𝒑𝒄 − 𝜺, … . . 𝟐. 𝟏𝟕
𝑷 𝑻′
𝒑𝒄 𝒑𝒄
𝑷′𝒑𝒄 = … . . 𝟐. 𝟏𝟖,
𝑻𝒑𝒄 + 𝑩 𝟏 − 𝑩 𝜺
𝑩 = 𝒎𝒐𝒍𝒆 𝒇𝒓𝒂𝒄𝒕𝒊𝒐𝒏 𝒐𝒇 𝑯𝟐 𝑺 𝒊𝒏 𝒕𝒉𝒆 𝒈𝒂𝒔 𝒎𝒊𝒙𝒕𝒖𝒓𝒆,
𝜺 = 𝟏𝟐𝟎 𝑨𝟎.𝟗 − 𝑨𝟏.𝟔 + 𝟏𝟓 𝑩𝟎.𝟓 − 𝑩𝟒.𝟎 … … 𝟐. 𝟏𝟗
𝑨 = 𝒚𝑯𝟐𝑺 + 𝒚𝑪𝑶𝟐 … … 𝟐. 𝟐𝟎
𝟏.𝟐𝟎𝟒𝟖
𝑷
𝑹𝒔 = 𝜸𝒈 + 𝟏. 𝟒 𝟏𝟎𝒙
𝟏𝟖. 𝟐
Where:
𝒙 = 𝟎. 𝟎𝟏𝟐𝟓𝑨𝑷𝑰 − 𝟎. 𝟎𝟎𝟎𝟗𝟏 𝑻 − 𝟒𝟔𝟎
𝑻 = 𝒕𝒆𝒎𝒑𝒆𝒓𝒂𝒕𝒖𝒓𝒆 𝒊𝒏 °𝑹,
𝒑 = 𝒔𝒚𝒔𝒕𝒆𝒎 𝒑𝒓𝒆𝒔𝒔𝒖𝒓𝒆 𝒊𝒏 𝒑𝒔𝒊𝒂,
𝜸𝒈 = 𝒔𝒐𝒍𝒖𝒕𝒊𝒐𝒏 𝒈𝒂𝒔 𝒔𝒑𝒆𝒄𝒊𝒇𝒊𝒄 𝒈𝒓𝒂𝒗𝒊𝒕𝒚,
𝑵𝑩, 𝑷 ≤ 𝑷𝒃
𝟏.𝟐𝟎𝟒𝟖
𝑷
𝑹𝒔 = 𝜸𝒈 + 𝟏. 𝟒 𝟏𝟎𝒙
𝟏𝟖. 𝟐
Where:
𝒙 = 𝟎. 𝟎𝟏𝟐𝟓𝑨𝑷𝑰 − 𝟎. 𝟎𝟎𝟎𝟗𝟏 𝑻 − 𝟒𝟔𝟎
𝟎.𝟖𝟑
𝑹𝒔 𝒂
𝑷𝒃 = 𝟏𝟖. 𝟐 𝟏𝟎 − 𝟏. 𝟒 ,
𝜸𝒈
𝟎.𝟖𝟑
𝑹𝒔 𝒂
𝑷𝒃 = 𝟏𝟖. 𝟐 𝟏𝟎 − 𝟏. 𝟒 ,
𝜸𝒈
𝑽𝟎 𝒑,𝑻
𝑩𝒐 = , 𝑩𝟎 ≥ 𝟏
𝑽𝒐 𝒔𝒄
Bo – Oil formation volume factor, bbl/STB
𝑩𝒐 = 𝒇(𝑹𝒔 , 𝜸𝒈 , 𝜸𝒐 , 𝑻)
Standing‟s Correlation for Bo
𝟎.𝟓 𝟏.𝟐
𝜸𝒈
𝑩𝒐 = 𝟎. 𝟗𝟕𝟓𝟗 + 𝟎. 𝟎𝟎𝟎𝟏𝟐𝟎 𝑹𝒔 + 𝟏. 𝟐𝟓 𝑻 − 𝟒𝟔𝟎
𝜸𝒐
T = Temperature in Deg. Rankine
Characteristics Trend
Characteristics Trends
Reservoir Pressure Remains High
Surface GOR Remains low
Starts Early and Increases to
Water Production
appreciable amounts
Flows until water production gets
Well behavior
excessive (liquid loading of wells)
Similarly, the reduction in the pore volume due to the expansion of the
reservoir rock is given by:
𝑁𝐵𝑜𝑖 (1 + 𝑚)
𝑐𝑎𝑛𝑔𝑒 𝑖𝑛 𝑝𝑜𝑟𝑒 𝑣𝑜𝑙𝑢𝑚𝑒 = 𝑐𝑓 ∆𝑝
1 − 𝑆𝑤𝑖
Combining equations from the volume occupied by initial oil, through total
pore volume with the equation of MBE written in a generalized form and
rearranging gives:
𝑁𝑝 𝐵𝑜 + (𝐺𝑝 −𝑁𝑝 𝑅𝑠 )𝐵𝑔 − (𝑊𝑒 −𝑊𝑝 𝐵𝑤 ) − 𝐺𝑖𝑛𝑗 𝐵𝑔𝑖𝑛𝑗 − 𝑊𝑖𝑛𝑗 𝐵𝑤
𝑁=
𝐵𝑔 𝑆𝑤𝑖 𝑐𝑤 + 𝑐𝑓
(𝐵𝑜 −𝐵𝑜𝑖 ) + (𝑅𝑠𝑖 −𝑅𝑠 )𝐵𝑔 + 𝑚𝐵𝑜𝑖 − 1 + 𝐵𝑜𝑖 (1 + 𝑚)[ ]∆𝑝
𝐵𝑔𝑖 1 − 𝑆𝑤𝑖
𝐺𝑝 = 𝑅𝑝 𝑁𝑝
𝑁𝑝 [𝐵𝑜 +(𝑅𝑝 −𝑅𝑠 )𝐵𝑔 ] − (𝑊𝑒 −𝑊𝑝 𝐵𝑤 ) − 𝐺𝑖𝑛𝑗 𝐵𝑔𝑖𝑛𝑗 − 𝑊𝑖𝑛𝑗 𝐵𝑤𝑖
𝑁=
𝐵𝑔 𝑆𝑤𝑖 𝑐𝑤 + 𝑐𝑓
(𝐵𝑜 −𝐵𝑜𝑖 ) + (𝑅𝑠𝑖 −𝑅𝑠 )𝐵𝑔 + 𝑚𝐵𝑜𝑖 − 1 + 𝐵𝑜𝑖 1 + 𝑚 [ ] ∆𝑝
𝐵𝑔𝑖 1 − 𝑆𝑤𝑖
Solution
Because the reservoir contains a gas cap, the rock and fluid expansion can be
neglected, i.e., set cf and cw = 0. For illustration purposes, however, the rock
and fluid expansion term will be included in the calculations.
Part A. Cumulative water influx
Step 1. Calculate cumulative gas-oil ration Rp:
1100 × 106
𝑅𝑝 = = 1100𝑠𝑐𝑓/𝑆𝑇𝐵
1 × 106
Step 2. Arrange Equation 11-17 to solve for We:
𝑁 = 𝑁𝑝 [𝐵𝑡 +(𝑅𝑝 −𝑅𝑠𝑖 )𝐵𝑔 ]
𝐵𝑔 𝑆𝑤𝑖 𝑐𝑤 + 𝑐𝑓
−𝑁 𝐵𝑔 − 𝐵𝑡𝑖 + 𝑚𝐵𝑡𝑖 − 1 + 𝐵𝑡𝑖 1 + 𝑚 ∆𝑝 + 𝑊𝑝 𝐵𝑤𝑝
𝐵𝑔𝑖 1 − 𝑆𝑤𝑖
0.00092
We = 106 [1.655 + (1100-1040)0.00092] -107 [(1.655 – 1.58) + 0.25(1.58) − 1
0.00080
0.2(1.5×10−6
+ 1.58(1 + 0.25) 1;0.2
(3000 − 2800)] + 50,000 = 411,281𝑏𝑏𝑙
Solution
Step 1. Assuming the same porosity and connate water for the oil
and gas
zones, calculate m:
20,000
𝑚= = 0.2
100,000
𝟕𝟕𝟓𝟖𝑨𝒉∅ 𝟏;𝑺𝒘𝒊
𝑰𝒏𝒊𝒕𝒊𝒂𝒍 𝑶𝑰𝑷, 𝑵 =
𝑩𝒐𝒊
𝟕𝟕𝟓𝟖𝑨𝒉∅ 𝟏;𝑺𝒈 ;𝑺𝒘𝒊
𝑶𝒊𝒍 𝑹𝒆𝒎𝒂𝒊𝒏𝒊𝒏𝒈 =
𝑩𝒐
𝟕𝟕𝟓𝟖𝑨𝒉∅ 𝟏;𝑺𝒘𝒊 𝟕𝟕𝟓𝟖𝑨𝒉∅ 𝟏;𝑺𝒈 ;𝑺𝒘𝒊
𝑹𝒆𝒔𝒆𝒓𝒗𝒆𝒔 = −
𝑩𝒐𝒊 𝑩𝒐
𝟏;𝑺𝒘𝒊 𝟏;𝑺𝒈 ;𝑺𝒘𝒊
𝑹𝒆𝒔𝒆𝒓𝒗𝒆𝒔 = 𝟕𝟕𝟓𝟖𝑨𝒉∅ −
𝑩𝒐𝒊 𝑩𝒐
𝟕𝟕𝟓𝟖𝑨𝒉∅ 𝟏;𝑺𝒘𝒊
𝑰𝒏𝒊𝒕𝒊𝒂𝒍 𝑶𝑰𝑷, 𝑵 =
𝑩𝒐𝒊
𝟕𝟕𝟓𝟖𝑨𝒉∅𝑺𝒐𝒓
𝑶𝒊𝒍 𝑹𝒆𝒎𝒂𝒊𝒏𝒊𝒏𝒈𝒕 = ,
𝑩𝒐𝒊
𝑷 ≅ 𝑷𝒊 𝒅𝒖𝒆 𝒕𝒐 𝒕𝒉𝒆 𝒂𝒒𝒖𝒊𝒇𝒆𝒓 𝒔𝒖𝒑𝒑𝒐𝒓𝒕
𝟕𝟕𝟓𝟖𝑨𝒉∅ 𝟏;𝑺𝒐𝒓 ;𝑺𝒘𝒊
𝑹𝒆𝒔𝒆𝒓𝒗𝒆𝒔𝒕 =
𝑩𝒐𝒊
CU PET 323, 2012/2013 TTOWG! 109
PET 323– Reservoir Engineering I – I. Steteyeobot and D. Alaigba
PART 8: VOLUMETRIC ESTIMATION OF HIP
Volumetric Estimation of Reserves and Recovery Factors for
Initially Under saturated Oil Reservoirs under Active Water
Drive
𝑰𝒏𝒊𝒕𝒊𝒂𝒍 𝑶𝑰𝑷;𝑶𝒊𝒍 𝑹𝒆𝒎𝒂𝒊𝒏𝒊𝒏𝒈
𝑹𝒆𝒄𝒐𝒗𝒆𝒓𝒚 𝑭𝒂𝒄𝒕𝒐𝒓𝒚, 𝑹𝑭 =
𝑰𝒏𝒊𝒕𝒊𝒂𝒍 𝑶𝑰𝑷
𝟕𝟕𝟓𝟖𝑨𝒉∅ 𝟏−𝑺𝒐𝒓 −𝑺𝒘𝒊
𝑹𝒆𝒔𝒆𝒓𝒗𝒆𝒔 𝑩𝒐𝒊
𝑹𝑭 = = 𝟕𝟕𝟓𝟖𝑨𝒉∅ 𝟏−𝑺𝒘𝒊
𝑰𝒏𝒊𝒕𝒊𝒂𝒍 𝑶𝑰𝑷
𝑩𝒐𝒊
𝟏;𝑺𝒐𝒓 ;𝑺𝒘𝒊
𝑹𝑭 =
𝟏;𝑺𝒘𝒊
𝑹𝒆𝒔𝒆𝒓𝒗𝒆𝒔 = 𝑮 − 𝑮𝒂 ,
𝒘𝒉𝒆𝒓𝒆 𝒂 𝒊𝒔 𝒕𝒉𝒆 𝒂𝒃𝒂𝒏𝒅𝒐𝒏𝒎𝒆𝒏𝒕 𝒑𝒓𝒆𝒔𝒔𝒖𝒓𝒆
𝟒𝟑𝟓𝟔𝟎𝑨𝒉∅ 𝟏;𝑺𝒘𝒊
𝑮𝒂 =
𝑩𝒈𝒂
𝟏 𝟏
𝑹𝒆𝒔𝒆𝒓𝒗𝒆𝒔 = 𝟒𝟑𝟓𝟔𝟎𝑨𝒉∅ 𝟏 − 𝑺𝒘𝒊 −
𝑩𝒈𝒊 𝑩𝒈𝒂
𝟏 𝟏
𝟒𝟑𝟓𝟔𝟎𝑨𝒉∅ 𝟏;𝑺𝒘𝒊 ;
𝑩𝒈𝒊 𝑩𝒈𝒂
𝑹𝑭 = 𝟒𝟑𝟓𝟔𝟎𝑨𝒉∅ 𝟏−𝑺𝒘𝒊
𝑩𝒈𝒊
𝑩𝒈𝒊
𝑹𝑭 = 𝟏 −
𝑩𝒈𝒂
𝟏;𝑺𝒘𝒊 ;𝑺𝒈𝒓
𝑹𝑭 =
𝟏;𝑺𝒘𝒊
DDI 0.440971033
SGI 0.235013558
b
WDI 0.308154427
EDI 0.015860982
Sum 1