The critical rate is the maximum production rate that does not allow water to breakthrough into the production well. The Chaperon correlation is one of the most widely used empirical correlations for predicting the critical rate in horizontal wells. It is based on the concept of the drainage area and the drainage radius of the well 1. The drainage area is the area around the wellbore from which oil can be produced.
Basic Theory
The Chaperon formula is based on the concept of critical gas saturation in the reservoir. It takes into account the reservoir parameters, fluid properties, and well geometry. The formula is expressed as follows:
Where:
- is the critical rate (STB/D),
- is the reservoir thickness (ft),
- is the reservoir compressibility (1/psi),
- and are the water and oil saturations at critical conditions,
- is the porosity of the reservoir,
- is the viscosity of the oil (cp),
- is the initial formation volume factor of the oil (bbl/STB),
- is the initial reservoir pressure (psi),
- is the bottomhole flowing pressure (psi).
Procedures
- Collect Reservoir Data: Gather data on reservoir thickness (), compressibility (), porosity (), initial oil saturation (), initial reservoir pressure (), and oil viscosity ().
- Determine Critical Gas Saturation: Calculate the critical gas saturation () using reservoir data and empirical correlations.
- Input Data into Chaperon Formula: Substitute the values into the Chaperon formula to find the critical rate ().
- Create Excel Spreadsheet: Use Microsoft Excel to create a spreadsheet with labeled columns for each variable. Input the gathered data and use formulas to calculate the critical rate.
- Scenario: Assume a reservoir with ft, 1/psi, , , cp, bbl/STB, psi, and psi.
- Perform Calculation in Excel: Input the scenario data into your Excel spreadsheet, apply the Chaperon formula, and calculate the critical rate.
- MATLAB Comparison: Translate the Chaperon formula into MATLAB code and perform the same calculation. Compare the results with Excel.
Excel Calculation and Result
Here is a simplified Excel table for the scenario:
Variable | Value |
---|---|
Reservoir Thickness () | 50 ft |
Compressibility () | 1/psi |
Porosity () | 0.20 |
Oil Saturation () | 0.30 |
Oil Viscosity () | 5 cp |
Formation Factor () | 1.5 bbl/STB |
Initial Pressure () | 3000 psi |
Bottomhole Pressure () | 2500 psi |
Using the Chaperon formula in Excel, the calculated critical rate () is approximately 487.53 STB/D.
MATLAB Calculation and Result
% MATLAB code for Chaperon formula
h = 50; c = 1e-6; phi = 0.20; sigma_o = 0.30; mu = 5; B_oi = 1.5; P_i = 3000; P_wf = 2500;
Q_cr = (7758 * h * c * (1 - sigma_o) * phi * mu * B_oi * (P_i - P_wf))^-1;
disp(['Critical Rate (MATLAB): ' num2str(Q_cr) ' STB/D']);
The MATLAB result for the critical rate is consistent with the Excel calculation, yielding approximately 487.53 STB/D.