The Hammerlindl method is used to estimate the original gas in place (OGIP) of abnormally pressured gas reservoirs. The method was developed by Hammerlindl and involves plotting the cumulative gas production against the pressure depletion. The slope of the straight line obtained from the plot is used to estimate the OGIP. The method is based on the assumption that the reservoir pressure is proportional to the gas saturation pressure.
Basic Theory
The Hammerlindl method is based on the principle of material balance. It assumes that the reservoir behaves as a single-phase system and considers the gas expansion and compression in the reservoir…
The formula for gas in place (GIP) using the Hammerlindl method is given by:
Procedures
- Input Data: Collect reservoir data including net pay thickness (), porosity (), formation volume factor ratio (), initial reservoir pressure (), wellbore flowing pressure (), and initial gas formation volume factor ().
- Formula Application: Apply the Hammerlindl formula to calculate gas in place ().
- Excel Implementation:
- Create an Excel table with columns for each input parameter and the formula for .
- Input the collected data into the corresponding columns.
- Use Excel formulas to calculate based on the Hammerlindl method.
- Scenario Example:
- Let’s consider a reservoir with the following data:
- Net pay thickness (): 50 feet
- Porosity (): 0.20
- Formation volume factor ratio (): 1.5
- Initial reservoir pressure (): 3000 psia
- Wellbore flowing pressure (): 500 psia
- Initial gas formation volume factor (): 0.005 reservoir cubic feet per standard cubic foot.
- Let’s consider a reservoir with the following data:
Excel Implementation
Parameter | Value |
---|---|
Net Pay Thickness () | 50 |
Porosity () | 0.20 |
Formation Volume Factor Ratio () | 1.5 |
Initial Reservoir Pressure () | 3000 |
Wellbore Flowing Pressure () | 500 |
Initial Gas Formation Volume Factor () | 0.005 |
MATLAB Comparison
% MATLAB code for Hammerlindl method
h = 50; % Net pay thickness (feet)
phi = 0.20; % Porosity
c = 1.5; % Formation volume factor ratio
pi = 3000; % Initial reservoir pressure (psia)
pwf = 500; % Wellbore flowing pressure (psia)
Bgi = 0.005; % Initial gas formation volume factor (reservoir scf/scf)
GIP_MATLAB = (7758 * h * phi * c * (pi - pwf)) / Bgi;
disp(['GIP_MATLAB: ', num2str(GIP_MATLAB), ' scf']);