Book summary “Petroleum production systems”: summary

Introduction

You face the complex challenge of managing continuous fluid flow from deep underground reservoirs to surface processing facilities. Petroleum production engineering focuses intensely on maximizing well deliverability and short-term recovery, contrasting sharply with broader reservoir engineering priorities.

Petroleum production systems integrate the subterranean reservoir with surface facilities, requiring coordinated inflow performance, well deliverability optimization, and targeted stimulation to maximize hydrocarbon recovery.

Michael J. Economides details the complete production system in this authoritative text. You will explore inflow behavior across various reservoir types, wellbore flow mechanics, artificial lift methods, well performance evaluation, stimulation techniques, and sand management.

To master these production operations, you must answer several critical questions:

• How do you quantify fluid flow through porous media and wellbores?
• What mechanisms cause near-wellbore formation damage, and how do you remove it?
• When should you apply gas lift versus pump-assisted lift?
• How do you design hydraulic fracturing and matrix acidizing treatments to boost deliverability?

1. The Role of Petroleum Production Engineering

Production engineering optimizes the continuous flow path from the subterranean reservoir through artificial structures to surface separation facilities. While reservoir engineering prioritizes long-term recovery and resource management, production engineering focuses on maximizing well deliverability in a cost-effective manner. This divergence in priorities can lead to organizational tension when international companies focused on rapid production must collaborate with national oil companies focused on long-term reserve management.

A petroleum production system consists of two intimately connected domains: the subterranean reservoir and the artificial surface and subsurface structures.

The production system integrates the porous reservoir medium with artificial structures including the wellbore, bottomhole and wellhead assemblies, and surface gathering, separation, and storage facilities. Production engineering addresses the entire continuous flow path from the reservoir up to the surface separation facility.

During management at Yukos E&P in Russia, appropriate production enhancement actions successfully increased field production by more than 15 percent, even after drilling operations were entirely suspended for up to a year.

System Components

Porosity, net reservoir thickness, and fluid saturations are the three critical variables required to evaluate the viability of a hydrocarbon prospect. Porosity measures the pore volume relative to bulk volume, net thickness defines the vertical extent of the productive porous medium, and fluid saturations determine the fractions of pore space occupied by water and hydrocarbons. Porosity values vary from over 0.3 to less than 0.1, and standard conditions are defined as usually 14.7 psi and 60 deg F. Electrical resistivity logs combined with spontaneous potential logs are standard tools used alongside core and well test data to identify sandstones, shales, and hydrocarbon-bearing zones.

Hydrocarbon reservoirs are classified thermodynamically using phase diagrams that map temperature and pressure boundaries. Reservoirs are categorized based on whether initial pressures sit above the bubble-point pressure in undersaturated oil systems, below it in saturated or two-phase systems, or beyond the dew-point and cricondentherm in gas and retrograde condensate systems. PVT properties determined from laboratory fluid samples dictate these phase behaviors.

Reservoir Classification and Flow Fundamentals

Fluid flow through a porous medium is governed by permeability, which quantifies the medium's capacity to transmit fluids under a pressure gradient. Originating from Henry Darcy's historical experiments with water flowing through sand packs, Darcy's law establishes that fluid velocity is linearly proportional to the pressure drop and inversely proportional to fluid viscosity. The ratio of permeability to viscosity is termed mobility.

The near-wellbore zone experiences disproportionate pressure drops and severe damage from drilling and completion activities, making its restoration critical. Due to radial convergence, up to 90 percent of the total pressure drop in the reservoir may be consumed in a zone just a few feet away from the well. Matrix stimulation, hydraulic fracturing, casing, cementing, perforation, and gravel packing are employed to mitigate this damage, provide zonal isolation, and prevent sand production.

Well deliverability is determined by the intersection of the inflow performance relationship and the vertical flow performance. The inflow performance relationship describes what the reservoir can deliver, while vertical flow performance describes what the well can lift to the surface. Combining these two relationships yields the operating well deliverability, which can be optimized by lowering bottomhole flowing pressure or increasing the dimensionless productivity index through stimulation and modern well architecture.

Production engineers utilize rate-versus-pressure relationships, well testing, and production logging to diagnose and resolve well performance issues. Diagnostic procedures involve analyzing flow paths from reservoir to separator, conducting well testing to evaluate flow potential and skin effects, and using permanent downhole instruments or production logging to map fluid distribution and completion problems. Remedial actions range from hydraulic fracturing to resizing surface flow lines.

Petroleum production engineering calculations predominantly employ oilfield units due to industry convention, requiring specific conversions for consistent SI application. In Example 1-1, steady-state radial flow calculations are demonstrated by converting SI data to oilfield units ($q = 0.001\text{ m}^3/\text{s}$, $B = 1.1\text{ res m}^3/\text{ST m}^3$, $\mu = 2 \times 10^{-3}\text{ Pa-s}$, $k = 10^{-14}\text{ m}^2$, $h = 10\text{ m}$, $r_e = 575\text{ m}$, $r_w = 0.1\text{ m}$, $s = 0$) and alternatively by deriving the constant $0.159$ for SI units, where the calculated pressure drawdown equals $3.043 \times 10^7\text{ Pa}$. To work with these systems:

• Convert all input parameters from SI units to oilfield units using standard conversion factors.
• Calculate the desired result using standard oilfield unit equations containing established constants.
• Convert the resulting output values back into SI units as required.

2. Production from Undersaturated and Saturated Reservoirs

Radial Darcy's law governs single-phase fluid flow from the reservoir into the wellbore, driven by pressure gradients and rock permeability. In 1856, Darcy formulated this fundamental relationship for porous media.

Radial Darcy's law dictates that flow rate increases with higher pressure gradients, greater permeability, and larger reservoir height.

The flow rate decreases with fluid viscosity increases under single-phase conditions. Steady-state well performance requires constant boundary and bottomhole pressures, approximated primarily in waterflooded reservoirs.

Example 2-1 calculates a steady-state production rate for a 640-acre drainage area. The system features a 2980 ft drainage radius, a 5651 psi outer boundary pressure, and a 4500 psi bottomhole flowing pressure.

Steady-State and Pseudosteady-State Flow

When production wells feel their boundaries under no-flow conditions, the reservoir enters a pseudosteady state where pressure declines uniformly everywhere. The radial diffusivity equation solution during pseudosteady state utilizes the volumetrically weighted average reservoir pressure determined from pressure buildup tests.

Example 2-4 calculates an average reservoir pressure of 4114 psi after a specific pressure drawdown in a 640-acre reservoir. Dietz shape factors account for irregular drainage shapes and asymmetrical well positions within drainage volumes.

Example 2-5 demonstrates that placing a well at the corner quadrant of a square drainage shape reduces its flow rate by 10% compared to a centrally located well. The shape factor drops from 30.9 for a well at the center of a square to 4.5 for a well at the center of a square quadrant.

The skin effect quantifies near-wellbore formation damage or stimulation, creating an additional steady-state pressure drop. Introduced by Van Everdingen and Hurst in 1949, the skin factor modifies the pressure drop equations. Positive values indicate damage and a reduced effective wellbore radius, whereas negative values indicate stimulation or hydraulic fracturing.

In a damaged well with $s = 10$, the effective radius drops to $4.5 \times 10^{-5}$ of the well radius. A fractured well with $s = -6$ expands the effective radius to 402 times the well radius.

Two-Phase Flow Correlations

Production from saturated oil reservoirs involves free gas release below the bubble-point pressure, governed by complex PVT property variations. Expanding free or solution gas is a much more efficient recovery mechanism than undersaturated oil expansion.

Example 3-1 computes a total formation volume factor of $1.61\text{ res bbl/STB}$ at 3000 psi from an initial bubble-point pressure of 4336 psi. This reflects a 30 MMSTB oil volume reduction as gas evolves from a solution gas-oil ratio of 800 SCF/STB at the bubble point.

Vogel introduced an empirical dimensionless inflow performance relationship in 1968 for solution-gas drive reservoirs operating below the bubble point. Normalized against absolute open flow potential, Vogel's equation eliminates the need for explicit relative permeability measurements by using an empirical curve shape for two-phase flow.

Fetkovich proposed an exponential well test approximation in 1973 to model two-phase flow deliverability using stabilized flow rate data. Example 3-8 applies Fetkovich's approximation across three test rates of 1875, 1427, and 860 STB/d at respective bottomhole pressures to determine deliverability parameters.

3. Production from Natural Gas Reservoirs

Gas gravity links almost all gas properties because physical conditions vary significantly. The molecular weight of air is taken as 28.97, reflecting approximately 79 percent nitrogen and 21 percent oxygen. Pure methane has a gravity of 0.55, whereas rich or heavy gas reservoirs can have a gravity of 0.75 or higher.

Natural gas production calculations require real gas law compressibility factors, non-Darcy turbulence coefficients, and real gas pseudopressures.

Example 4-1 calculates the gas gravity for a natural gas mixture containing 88 percent C1, 8.2 percent C2, 2.1 percent C3, and 1.7 percent CO2, resulting in a molecular weight contribution of 18.225 and a gas gravity of 0.63.

Gas Properties and PVT Behavior

The real gas law incorporates a compressibility factor, or gas deviation factor, to account for deviations of real gas mixtures from ideal gas behavior at elevated pressures and temperatures. The universal gas constant R is 10.73 psi ft3/lb-mol-deg-R, and the compressibility factor Z_sc equals 1 at standard conditions of 14.7 psi and 60 degrees Fahrenheit. The Dranchuk and Abou-Kassem correlation reproduces Standing-Katz charts with an accuracy better than 1 percent.

Example 4-2 calculates the volume of 1 lb-mol of a natural gas mixture at 180 degrees Fahrenheit and 4000 psi, yielding a compressibility factor Z of 0.855 from the Standing-Katz chart and 0.89 from the Dranchuk-Abou-Kassem correlation.

Gas Well Deliverability

High-rate gas wells experience turbulence skin effects that require incorporating a non-Darcy coefficient into the deliverability equation to prevent severe overestimation of production capacity. Example 4-7 develops a deliverability curve showing that taking into account non-Darcy effects yields an absolute open-flow potential of 460 MSCF/d, whereas considering only Darcy flow incorrectly predicts 1420 MSCF/d. The non-Darcy coefficient used is 4.9 times 10 to the power of minus 2.

To evaluate non-Darcy flow deliverability, you must execute these steps:

• Calculate the time required for stabilized flow using drainage radius and formation properties.
• Calculate coefficients a and b for the deliverability equation.
• Graph the Darcy flow result and the correct curve incorporating non-Darcy effects on a log-log plot.

4. Production from Horizontal Wells

Horizontal wells extend contact with the reservoir, changing the drainage pattern and making them successful in thin, heavy oil, tight, or heterogeneous formations. Since the 1980s, horizontal wells have captured an ever-increasing share of hydrocarbon production.

Horizontal wells extend reservoir contact and alter drainage patterns, requiring specialized steady-state and pseudosteady-state inflow models.

Horizontal wells change flow from radial in vertical wells to a combination of radial, linear, and elliptical flows. They help reduce water or gas coning. You must evaluate thin reservoirs with a thickness threshold of 50 ft and tight formations with a permeability threshold of 0.1 md. Sufficient vertical permeability is required, and you must evaluate horizontal-to-vertical permeability anisotropy because drilling normal to the larger horizontal permeability direction improves production.

Steady-State Models

The steady-state model by Joshi derives horizontal well inflow by combining flow resistance in the horizontal plane and the vertical plane while incorporating anisotropy. Joshi treated horizontal well inflow as a combination of two-dimensional flow in the x-y plane with elliptical isobars and radial flow in the y-z plane from a vertical boundary located at half the reservoir thickness from the well.

A 2000-ft long horizontal lateral in a 100-ft thick reservoir with 10 md horizontal permeability and 1 md vertical permeability yields 1149 STB/d at a 2000 psi bottomhole pressure. The well has a diameter of 6 inches, a pressure at the drainage boundary of 4000 psia, an oil viscosity of 5 cp, and a formation volume factor of 1.1.

Furui et al. modeled horizontal well inflow assuming radial flow near the wellbore that transitions to linear flow farther away in a box-shaped drainage boundary. Total pressure drop is calculated by summing the pressure drop in the radial flow region near the well and the pressure drop in the linear flow region farther out.

Pseudosteady-State Flow and Multilateral Wells

The Babu and Odeh model uses a box-shaped drainage area to calculate pseudosteady-state horizontal well performance using shape factors and partial penetration skin factors. The model applies radial flow in the y-z plane adjusted by a geometry shape factor, while flow from beyond the wellbore in the x direction is accounted for through a partial penetration skin factor.

For a 4000-ft long reservoir with a 2000-ft long horizontal well, a width of 1414 ft, a height of 100 ft, and a bottomhole pressure of 2000 psi, the model predicts a flow rate of 1527 STB/d. The productivity index ratio of the horizontal well to the vertical well is 4.6.

Economides et al. developed a general model for arbitrary well trajectories by integrating unit length point sources within no-flow boundary boxes. The model decomposes three-dimensional dimensionless pressure into one two-dimensional part using horizontal shape factors and one one-dimensional part for vertical flow effects.

Horizontal gas well inflow equations are derived by substituting the oil formation volume factor with the gas formation volume factor and applying real gas pseudopressure. Equations for gas wells adapt oil models by using the real gas law or the real gas pseudopressure function introduced by Al-Hussainy and Ramey to handle pressure-dependent gas properties.

Multilateral wells feature two or more laterals contacting different parts of a reservoir, helping reduce wellbore frictional pressure drop compared to long single laterals. Multilateral wells are applied in coal production, heavy oil, unconventional low-permeability gas formations, and offshore fields. Flow rates from each lateral cannot simply be added together because junction pressures must balance across all linked laterals. Complex multilateral wells are applied in heavy oil recovery in the Orinoco Heavy-Oil-Belt in Venezuela.

5. Skin Effects and Formation Damage

Near-wellbore formation damage creates severe pressure drops quantified by skin factors and described by Hawkins formula.

Introduced in 1949, the Van Everdingen and Hurst skin effect uses a dimensionless skin factor $s$ analogous to a heat transfer film coefficient. A positive skin reduces well productivity, whereas a negative skin increases it.

Derived by Hawkins in 1956, Hawkins formula relates the damage skin effect to the radius of damage, undamaged permeability, and damaged zone permeability. Calculations show that permeability impairment has a much larger effect on the skin value than the depth of damage penetration. Example 6-1 demonstrates that a damage penetration of 3 ft with $k/k_s = 5$ yields a skin of 9.3, whereas $k/k_s = 10$ yields a skin of 20.9.

Skin Components

Apparent skin effects can be driven by high flow rates causing turbulence or by near-wellbore phase changes. Rate-dependent skin equals $Dq$ where $D$ is the non-Darcy coefficient, which you can isolate using multi-rate well tests. Gas condensate reservoirs suffer permanent permeability impairment from liquid condensate deposition that is not eliminated following shut-in. Partial completion creates a positive skin effect due to flow convergence, whereas well deviation introduces a negative skin effect.

Horizontal wells experience anisotropic flow and varying exposure times that alter damage skin distributions. Standard Hawkins formula fails for horizontal wells because radial flow is subject to horizontal-to-vertical permeability anisotropy. Cased and perforated completions generate mechanical skins that depend on perforation geometry, density, and phase angle. Gravel pack completions prevent sand production in unconsolidated formations but introduce series flow resistances.

Damage Mechanisms and Operational Sources

Formation damage arises from particle plugging, fines migration, chemical precipitation, fluid effects, and mechanical crushing. The Prudhoe Bay field experiences calcium carbonate precipitation from bicarbonate-rich connate water as pressure declines.

Drilling, completion, production, and injection operations are the primary operational sources of formation damage.

• Ensure aqueous fluids contacting formations contain at least 2 wt percent of KCl or at least one-tenth divalent cations to prevent clay dispersion.
• Filter completion fluids so they contain no more than 2 ppm of solids smaller than 2 microns.
• Use underbalanced or extreme overbalanced perforating techniques to keep debris away from perforation tunnels.

6. Wellbore Flow Performance

Production engineering optimizes the continuous flow path from the reservoir to surface facilities. Wellbore flow calculations combine mechanical energy balances for single-phase fluids with empirical multiphase holdup and pressure gradient correlations.

Single-Phase Flow

Single-phase pipe flow is classified as laminar or turbulent based on the dimensionless Reynolds number. The general Reynolds number value where transition from laminar to turbulent flow occurs in circular pipes is 2100.

Example 7-1 constructs a graph of Reynolds number versus volumetric flow rate for 1.03-specific gravity water in a 7-in., 32-lb/ft casing. You find that the laminar-to-turbulent transition occurs at a flow rate of 81 bbl/d.

Overall pressure drop in single-phase pipe flow is the sum of potential energy, kinetic energy, and frictional contributions. Example 7-2 calculates a potential energy pressure drop of -292 psi for 1000 bbl/d of brine injected through tubing deviated 50 degrees from vertical over a distance of 1000 ft. Calculating pressure drop in gas wells requires integrating the mechanical energy balance while accounting for fluid compressibility. Example 7-5 uses a tubing length of 10,000 ft and a gas production rate of 2 MMSCF/d to calculate a bottomhole pressure of 1078 psia.

Multiphase Flow Correlations

In multiphase wellbore flow, the lighter phase typically moves faster than the denser phase, causing the denser phase to be held up. Example 7-6 determines that for a gas-liquid flow with a slip velocity of 60 ft/min and superficial velocities of 60 ft/min for each phase, the liquid holdup is 0.62 and the gas holdup is 0.38.

Gas-liquid upward flow exhibits distinct flow regimes depending on phase flow rates and fluid properties. Example 7-7 places flow rates of 200 bbl/d of water and 10,000 ft3/day of air in high-velocity slug or churn flow with a liquid velocity number of 1.11 and a gas velocity number of 9.8.

The modified Hagedorn and Brown method is an empirical two-phase flow correlation. Example 7-8 calculates a liquid holdup of 0.46 for 2000 bbl/d of oil and 1 MMSCF/d of gas. Stepwise pressure traverse calculations use a fixed-length increment, such as 200 ft for tubing flow, across a depth of 10,000 ft to determine overall pressure drops.

1. Select a length increment, such as 200 ft for tubing flow.
2. Estimate the pressure drop using the no-slip average density and potential energy pressure gradient.
3. Calculate all fluid properties at the average pressure and average temperature.
4. Calculate the pressure gradient using a two-phase flow correlation.
5. Obtain a new estimate of pressure drop from the calculated pressure gradient and length increment.
6. Repeat the procedure from step 3 if the new and old pressure drop estimates do not match within tolerance.

7. Flow in Horizontal Wellbores, Wellheads, and Gathering Systems

Surface gathering networks, wellhead chokes, and horizontal wellbore friction losses dictate operating wellhead and bottomhole pressures.

Horizontal Pipe Flow

Single-phase liquid flow in horizontal pipes is calculated using mechanical energy balance equations where potential energy pressure drop is zero. For an incompressible fluid in a constant diameter horizontal pipe, both potential and kinetic energy pressure drops are zero, simplifying the mechanical energy balance. Friction factors are determined using the Chen equation or the Moody diagram based on the Reynolds number.

Calculating the pressure drop for 1000 bbl/d of injection water through a 3000-ft long, 1 1/2-in. I.D. galvanized iron pipe with relative roughness 0.004 yielded a significant pressure loss of over 3000 ft. High-rate, low-pressure horizontal gas lines require incorporating a kinetic energy term into the compressible fluid pressure drop calculations. Integrating the mechanical energy balance for real gas yields an implicit equation containing an additional natural logarithm term to account for the kinetic energy pressure drop.

Horizontal two-phase gas-liquid flow regimes are divided into segregated, intermittent, and distributive categories, which are predicted using specific flow regime maps. Maps developed by Baker, Mandhane et al., Beggs and Brill, and Taitel and Dukler are used to predict these regimes. Pressure gradients in horizontal gas-liquid flow are commonly calculated using empirical correlations developed by Beggs and Brill, Eaton et al., and Dukler.

Chokes and Gathering Systems

Wellhead chokes restrict production rates, and when fluid velocity reaches sonic velocity in the throat, critical flow occurs where downstream pressure changes do not affect flow rate. Surface gathering systems collect flow from multiple wells into central processing facilities, where pressure drops through lines, valves, and fittings dictate wellhead operating pressures.

Pressure drop along a horizontal wellbore is negligible in low-permeability reservoirs but becomes a major productivity-limiting factor in high-permeability formations. Evaluating a 4000-ft long, 4-in. diameter horizontal well showed that for 100 md permeability the wellbore pressure drop is only 8 percent of drawdown, but for 1000 md permeability with 100 psi drawdown, it consumes 78 percent of the reservoir drawdown. Minor undulations in nominally horizontal wellbores cause severe changes in multiphase flow behavior, accumulating stagnant liquids in sumps and increasing potential energy pressure drops.

8. Well Deliverability and Production Forecasting

Combining inflow performance relationships with vertical flow performance curves yields operating deliverability matched with material balance forecasting.

You determine well deliverability graphically by finding the intersection point between the reservoir inflow performance relationship and the vertical flow performance curve. The inflow performance relationship curve shows what the reservoir delivers into the well at a specific time. The vertical flow performance curve relates wellhead flowing pressure to bottomhole flowing pressure based on hydrostatic and frictional pressure drops.

Example 9-1 combines a single-phase inflow performance relationship with a flat vertical flow performance curve for an 8000-ft well. You find an expected production rate of 430 STB/d at a bottomhole pressure of 3055 psi. Increasing the flowing gas-liquid ratio reduces bottomhole flowing pressure and shifts gradient curves.

Reservoir inflow performance is degraded by formation damage and by declining reservoir pressure over time. A zero skin effect yields 577 STB/d, but skin values of 10 and 50 reduce the flow rate to 273 and less than 73 STB/d respectively.

Material Balance Forecasting

Fluid recovery in undersaturated oil reservoirs depends entirely on fluid expansion driven by total compressibility. Total compressibility combines oil, water, and rock compressibilities weighted by their respective saturations.

Example 10-2 evaluates a 40-acre oil well under pseudosteady-state conditions. After 44 days and a 1650 psi pressure drop, total recovery reaches only 2.15% of original oil-in-place. Tarner's method provides a stepwise calculation procedure for solution gas-drive reservoirs to predict pressure decline, produced gas-oil ratios, incremental oil recovery, and time duration.

You follow a strict procedure to forecast production from a two-phase reservoir:

• Define pressure interval delta p for calculation.
• Calculate phi_n and phi_g using properties at average pressure in the interval.
• Assume a value of producing gas-oil ratio guess in the interval.
• Calculate incremental cumulative oil production.
• Calculate incremental gas cumulative production.
• Calculate the oil saturation.
• Obtain relative permeability ratio from a plot versus oil saturation.
• Calculate calculated producing gas-oil ratio.
• Compare guessed and calculated gas-oil ratios and repeat if necessary.

Example 10-4 forecasts performance for a 40-acre two-phase reservoir under solution gas drive. Recovery is approximately seven times higher than in a comparable undersaturated reservoir over a similar pressure drop.

9. Gas Lift and Pump-Assisted Lift

Artificial lift systems restore or enhance production by injecting gas to lighten fluid columns or employing downhole pumps to boost pressure.

Gas lift reduces the hydrostatic pressure component in the wellbore by injecting gas continuously or intermittently, lowering the density of the fluid mixture. Continuous gas-lift design balances the required bottomhole flowing pressure with the artificial gas-liquid ratio and the injection pressure down the annulus. Example 11-1 calculates the required working GLR of 400 SCF/STB for a well with 8000 ft depth and 1500 psi bottomhole pressure, using a natural GLR of 300 SCF/STB.

Newly drilled wells filled with high-density completion fluids are unloaded in steps using multiple gas-lift valves starting from the topmost valve. Example 11-5 designs an unloading procedure starting with a kickoff pressure of 1000 psi and a first valve location calculated at 2855 ft using a kill fluid gradient of 0.4 psi/ft.

Follow these steps for unloading wells:

• Inject gas through the first valve to lift the liquid column from the valve location to the surface.
• Close the first valve and open the next valve below to lift the column from the second valve location upward.
• Repeat the procedure until the entire wellbore is unloaded, using a 100-200 psi safety factor in design tubing pressure for subsequent valves.

Downhole Pumps and ESPs

Downhole pumps boost well productivity by increasing pressure at the bottom of the tubing to lift liquids to the surface. Electrical submersible pumps are multistage centrifugal pumps capable of producing high liquid volumes from deep wells by imparting a constant pressure head. Example 12-5 details the complete design of an ESP for a 10,000-ft-deep well producing 3000 STB/d of saturated oil, calculating a required pump pressure increase of 3004 psi and selecting a 100-stage pump.

10. Well Performance Evaluation, Stimulation, and Sand Management

Production optimization relies on open-hole wireline logs to evaluate static properties like lithology, porosity, and fluid saturations. When wells exhibit abnormally low productivity indices, production logging identifies damaged intervals for targeted stimulation.

Production logging, matrix acidizing, hydraulic fracturing, and sand control safeguard well productivity and maximize ultimate recovery.

Example 13-2 analyzes Well A-1 in Reservoir Alpha using temperature and spinner flowmeter logs to reveal a skin factor of 20, showing zone A was severely damaged while zone B produced 70 percent of the total rate.

Well Evaluation and Matrix Acidizing

Matrix acidizing dissolves near-wellbore minerals below parting pressure. Stoichiometric equations determine acid requirements.

For sandstone acid treatments:

• assume the HCl-carbonate reaction is very fast so the acid front is sharp;
• calculate the acid volume needed to consume the calcite divided by volumetric dissolving power;
• add the volume of pore space left in the region after carbonate removal.

Diversion techniques ensure untreated zones receive treatment. When placing diverting agents:

• use mechanical isolation with opposed cup packers, squeeze packers with retrievable bridge plugs, or inflatable straddle packers for zonal isolation;
• employ twice as many dense ball sealers as perforations, or a 50% excess when using buoyant ball sealers;
• add particulate diverting agents continuously or in batches to form a low-permeability filter cake easily removed after treatment.

Hydraulic Fracturing and Sand Management

Hydraulic fracturing creates conductive paths in tight formations. Unified Fracture Design optimizes fracture dimensions using proppant and drainage volumes.

Sand management requires exclusion or control. When designing sand exclusion:

• obtain a representative sample of formation material from rubber-sleeve cores, conventional cores, or sidewall cores;
• perform a sieve analysis using a series of standard sieve trays to measure formation particle size distribution;
• plot cumulative weight of formation material retained versus grain size on a semilog plot.

10 Key Ideas

1. A petroleum production system consists of two intimately connected domains: the subterranean reservoir and the artificial surface and subsurface structures.
2. Radial Darcy's law dictates that flow rate increases with higher pressure gradients, greater permeability, and larger reservoir height.
3. Natural gas production calculations require real gas law compressibility factors, non-Darcy turbulence coefficients, and real gas pseudopressures.
4. Horizontal wells extend reservoir contact and alter drainage patterns, requiring specialized steady-state and pseudosteady-state inflow models.
5. Near-wellbore formation damage creates severe pressure drops quantified by skin factors and described by Hawkins formula.
6. Wellbore flow calculations combine mechanical energy balances for single-phase fluids with empirical multiphase holdup and pressure gradient correlations.
7. Surface gathering networks, wellhead chokes, and horizontal wellbore friction losses dictate operating wellhead and bottomhole pressures.
8. Combining inflow performance relationships with vertical flow performance curves yields operating deliverability matched with material balance forecasting.
9. Artificial lift systems restore or enhance production by injecting gas to lighten fluid columns or employing downhole pumps to boost pressure.
10. Production logging, matrix acidizing, hydraulic fracturing, and sand control safeguard well productivity and maximize ultimate recovery.