Understanding Horizontal Well Critical Rate Correlation: Chaperon Formula in Excel

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:

    \[ Q_{cr} = \frac{7758 \cdot h \cdot c \cdot (\sigma_w - \sigma_o)}{\phi \cdot \mu \cdot B_{oi} \cdot (P_i - P_{wf})} \]

Where:

  • Q_{cr} is the critical rate (STB/D),
  • h is the reservoir thickness (ft),
  • c is the reservoir compressibility (1/psi),
  • \sigma_w and \sigma_o are the water and oil saturations at critical conditions,
  • \phi is the porosity of the reservoir,
  • \mu is the viscosity of the oil (cp),
  • B_{oi} is the initial formation volume factor of the oil (bbl/STB),
  • P_i is the initial reservoir pressure (psi),
  • P_{wf} is the bottomhole flowing pressure (psi).

Procedures

  1. Collect Reservoir Data: Gather data on reservoir thickness (h), compressibility (c), porosity (\phi), initial oil saturation (\sigma_o), initial reservoir pressure (P_i), and oil viscosity (\mu).
  2. Determine Critical Gas Saturation: Calculate the critical gas saturation (\sigma_w) using reservoir data and empirical correlations.
  3. Input Data into Chaperon Formula: Substitute the values into the Chaperon formula to find the critical rate (Q_{cr}).
  4. 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.
  5. Scenario: Assume a reservoir with h = 50 ft, c = 10^{-6} 1/psi, \phi = 0.20, \sigma_o = 0.30, \mu = 5 cp, B_{oi} = 1.5 bbl/STB, P_i = 3000 psi, and P_{wf} = 2500 psi.
  6. Perform Calculation in Excel: Input the scenario data into your Excel spreadsheet, apply the Chaperon formula, and calculate the critical rate.
  7. 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 (h) 50 ft
Compressibility (c) 10^{-6} 1/psi
Porosity (\phi) 0.20
Oil Saturation (\sigma_o) 0.30
Oil Viscosity (\mu) 5 cp
Formation Factor (B_{oi}) 1.5 bbl/STB
Initial Pressure (P_i) 3000 psi
Bottomhole Pressure (P_{wf}) 2500 psi

Using the Chaperon formula in Excel, the calculated critical rate (Q_{cr}) 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.

Comments

No comments yet. Why don’t you start the discussion?

Leave a Reply

Your email address will not be published. Required fields are marked *