Nodal Analysis (IPR + VLP)
Spreadsheet
29 rows x 4 columns
| A | B | C | D | |
|---|---|---|---|---|
| 1 | Nodal Analysis | |||
| 2 | Reservoir | |||
| 3 | Productivity Index (J) | 5 | bbl/(D.psi) | |
| 4 | Reservoir Pressure (Pr) | 3500 | psi | |
| 5 | Bubble Point (Pb) | 3000 | psi | |
| 6 | GOR | 500 | scf/STB | |
| 7 | Gas Specific Gravity | 0.75 | dimensionless | |
| 8 | Wellhead Pressure (Pwh) | 200 | psi | |
| 9 | Temperature | 200 | degF | |
| 10 | Liquid Density | 50 | lb/ft3 | |
| 11 | Liquid Viscosity | 2 | cP | |
| 12 | Interfacial Tension | 20 | dynes/cm | |
| 13 | Tubing ID | 2.992 | in | |
| 14 | TVD | 8000 | ft | |
| 15 | Pipe Roughness | 0.0006 | dimensionless | |
| 16 | ||||
| 17 | Nodal Analysis Table | |||
| 18 | Pwf (psi) | q IPR (bbl/D) | Qg (mmscf/D) | Pwf VLP (psi) |
| 19 | 3000 | 2500 | 1.25 | 1678.526921 |
| 20 | 2500 | 4814.814815 | 2.407407407 | 2232.356748 |
| 21 | 2000 | 6759.259259 | 3.37962963 | 2809.781682 |
| 22 | 1500 | 8333.333333 | 4.166666667 | 3311.32176 |
| 23 | 1000 | 9537.037037 | 4.768518519 | 3714.089468 |
| 24 | 500 | 10370.37037 | 5.185185185 | 4003.091664 |
| 25 | ||||
| 26 | Summary | |||
| 27 | Operating Rate | 5297.834483 | bbl/D | |
| 28 | Operating Pwf | 2375.794942 | psi | |
| 29 | AOF (at Pwf = 0) | 10833.33333 | bbl/D |
| Name | Refers To | Comment |
|---|---|---|
| J | $B$3 | Productivity index, bbl/(D.psi) |
| Pr | $B$4 | Average reservoir pressure, psia |
| Pb | $B$5 | Bubble point pressure, psia |
| GOR | $B$6 | Producing GOR, scf/STB |
| SG | $B$7 | Gas specific gravity |
| Pwh | $B$8 | Wellhead pressure, psia |
| Tf | $B$9 | Flowing temperature, degF |
| Rho_l | $B$10 | Liquid density, lb/ft3 |
| Ul | $B$11 | Liquid viscosity, cP |
| IFT | $B$12 | Interfacial tension, dynes/cm |
| ID | $B$13 | Tubing ID, in |
| TVD | $B$14 | True vertical depth, ft |
| Rough | $B$15 | Pipe roughness, dimensionless |
Description
Determine well operating point by intersecting the IPR curve (Vogel) with the VLP curve (Beggs-Brill multiphase tubing correlation). The node is at bottomhole — sweep flowing pressure to compute IPR rates and corresponding VLP bottomhole pressures. The operating point occurs where IPR Pwf equals VLP Pwf at the same rate.
Operating point: The intersection of IPR and VLP curves gives the sustainable production rate and flowing pressure for the current reservoir and completion conditions.
Workflow
- Inputs: Reservoir (J, Pr, Pb), completion (tubing ID, TVD), fluid (GOR, SG_gas, ρl, μl, σ), surface (Pwh)
- Step 1: Sweep Pwf and compute IPR rate using Vogel
- Step 2: For each rate, compute VLP bottomhole pressure using Beggs-Brill Pin
- Step 3: Find intersection — operating rate and flowing pressure
- Output: Operating rate, flowing pressure, AOF (absolute open flow potential)
How to use this blueprint
- In Excel, go to the Petroleum Office ribbon tab and click Blueprint Manager
- Search for Nodal Analysis (IPR + VLP)
- Click on the blueprint to preview
- Click Insert to place it into your worksheet. Modify the input values to match your data.