10.1 Petroleum geophysics and seismic exploration

Petroleum geophysics and seismic exploration are the primary tools used to find and develop oil and gas reservoirs. These methods use controlled sound waves to create detailed images of underground rock layers, allowing geophysicists to identify structures that may trap hydrocarbons.

Seismic data interpretation, combined with well log measurements, forms the backbone of modern exploration workflows. Understanding how these datasets work together is essential for mapping reservoirs, estimating reserves, and planning drilling operations.

Seismic Methods for Petroleum Exploration

Principles of Seismic Reflection and Refraction

Seismic reflection and refraction are the two foundational methods for imaging subsurface geology. Both rely on generating seismic waves with controlled sources (vibroseis trucks on land, air guns offshore, or explosives) and recording how those waves interact with rock layers below.

Seismic reflection occurs when waves hit an interface between two layers with different acoustic impedances and bounce back to the surface. Acoustic impedance is the product of a layer's density and seismic velocity ($Z = \rho \cdot v$). The greater the contrast in acoustic impedance between two layers, the stronger the reflection. This is what produces the visible horizons on a seismic section.

Seismic refraction occurs when waves encounter a layer with higher velocity, causing them to bend and travel along the interface before returning to the surface. The bending follows Snell's Law, which relates the angles of incidence and refraction to the velocities of the two layers:

$\frac{\sin \theta_1}{v_1} = \frac{\sin \theta_2}{v_2}$

Refraction is particularly useful for determining the velocity structure of the subsurface, which feeds into depth calculations.

Two-way travel time (TWT) is the time it takes for a seismic wave to travel from the source down to a reflector and back to the surface. Combined with a velocity model (built from refraction data and well logs), TWT is converted to actual depth, giving you the geometry of geological structures.

Applications and Advantages of Seismic Methods

Reflection seismic dominates petroleum exploration because it offers higher resolution and can image complex structures like faults, folds, and stratigraphic traps. High-resolution surveys can detect thin beds and subtle features down to several kilometers depth.

Refraction seismic plays a complementary role. It's most valuable for characterizing complex near-surface geology and providing velocity information used for static corrections and depth conversion of reflection data. It also sees use in deep crustal studies.

Both methods are non-invasive and can cover large areas efficiently, making them cost-effective compared to drilling exploratory wells. Advances in acquisition hardware, processing algorithms, and interpretation software continue to push the resolution and accuracy of subsurface imaging.

Interpreting Seismic Data for Reservoirs

Seismic Data Interpretation Techniques

Seismic sections are visual cross-sections of the subsurface, displaying two-way travel time on the vertical axis and distance along a survey line on the horizontal axis. They can be shown as wiggle traces, variable density plots, or color-coded displays.

The primary goal of interpretation is identifying hydrocarbon traps, which are geological configurations that accumulate and seal oil or gas. The three main trap types are:

• Anticlines: Upward-folded structures where hydrocarbons migrate to the crest and become trapped beneath a seal rock.
• Fault traps: Faults juxtapose permeable reservoir rock against impermeable rock, creating a lateral seal.
• Stratigraphic traps: Lateral changes in rock type or pinchouts of permeable layers create traps without structural deformation.

Certain patterns on seismic sections can directly indicate hydrocarbons:

• Bright spots: High-amplitude anomalies that may indicate gas or light oil in the pore space (the gas lowers acoustic impedance, increasing the reflection contrast).
• Flat spots: Horizontal reflections that cut across dipping structure, representing fluid contacts (gas-oil or oil-water boundaries).
• Dim spots / gas chimneys: Zones of reduced amplitude above a reservoir, caused by gas leaking upward and absorbing seismic energy.

Seismic attributes extract additional information from the data beyond simple reflection patterns:

• Amplitude attributes (e.g., RMS amplitude) highlight lithology or fluid changes.
• Frequency attributes (e.g., instantaneous frequency) can indicate variations in bed thickness or fluid type.
• Phase attributes help identify stratigraphic features and discontinuities.

Integrated Interpretation and Uncertainty Reduction

Seismic facies analysis classifies regions of a seismic section by their distinct reflection patterns, amplitude, frequency, and continuity. These facies can then be tied to depositional environments (channels, submarine fans, reefs) and lithologies (sandstone, shale, carbonate).

No single data type tells the full story. Interpretation improves significantly when seismic data is integrated with:

• Well logs: Detailed rock property and fluid measurements at specific borehole locations.
• Core analysis: Direct laboratory measurements of porosity, permeability, and fluid saturation from rock samples.
• Regional geological models: Broader structural and depositional context that guides and validates seismic picks.

Two advanced techniques are particularly important for reducing uncertainty:

1. Seismic inversion converts reflection data into quantitative rock properties (acoustic impedance or velocity), which can be directly compared with well log measurements.
2. AVO (Amplitude Versus Offset) analysis examines how reflection amplitude changes with source-receiver distance. Different fluids and lithologies produce characteristic AVO responses, making it a powerful tool for distinguishing gas-filled rock from brine-filled rock.

Reservoir Characterization with Logging and Attributes

Well Logging for Reservoir Properties

Well logging measures physical properties of rock formations along the length of a borehole, providing high-resolution vertical profiles of the reservoir. Logging tools are lowered into the well on a wireline cable, though modern techniques also include logging while drilling (LWD) and measurement while drilling (MWD), which record data in real time as the well is drilled.

The most common log types and what they measure:

These logs are combined to determine the four key reservoir parameters:

• Lithology: Identified by cross-plotting gamma ray, density, and neutron logs.
• Porosity: Estimated from density, neutron, and sonic logs (often using two or more for cross-checking).
• Permeability: Estimated from porosity using empirical relationships (e.g., the Kozeny-Carman equation or core-calibrated transforms).
• Fluid saturation: Inferred primarily from resistivity logs using Archie's equation: $S_w = \left(\frac{a \cdot R_w}{\phi^m \cdot R_t}\right)^{1/n}$, where $S_w$ is water saturation, $R_w$ is formation water resistivity, $R_t$ is true formation resistivity, $\phi$ is porosity, and $a$, $m$, $n$ are empirical constants.

Seismic Attributes and Reservoir Characterization

Seismic attributes are quantitative values extracted from seismic trace data using mathematical algorithms. They can be computed along interpreted horizons, on time slices, or across entire 3D volumes.

Key attribute categories:

• Amplitude attributes (RMS amplitude, instantaneous amplitude): Highlight acoustic impedance contrasts tied to lithology or fluid changes.
• Frequency attributes (instantaneous frequency, dominant frequency): Sensitive to bed thickness variations and fluid type.
• Phase attributes (instantaneous phase, cosine of phase): Useful for identifying stratigraphic features and lateral discontinuities.
• Coherence attributes: Measure trace-to-trace similarity, making faults, fractures, and channel edges stand out clearly.
• Curvature attributes: Quantify the degree of folding or bending of reflectors, highlighting structural and stratigraphic features.

Seismic inversion deserves special attention because it bridges the gap between seismic data and rock properties. Two main approaches exist:

• Deterministic inversion: Uses a single starting model and produces one "best-fit" output of acoustic impedance or velocity.
• Stochastic inversion: Uses a range of input models and generates multiple equally probable realizations, capturing uncertainty in the result.

The real power comes from integrating well logs and seismic attributes together. Well logs give you high-resolution vertical detail at discrete points; seismic attributes give you spatially continuous coverage across the entire field. Geostatistical methods like kriging and co-kriging combine both data types to build 3D reservoir models that honor the well data while using seismic trends to interpolate between wells.

3D and 4D Seismic in Field Development

3D Seismic Surveys and Interpretation

3D seismic surveys acquire data across a dense grid of source and receiver lines, producing a full three-dimensional volume of the subsurface rather than isolated 2D cross-sections. The closely spaced lines provide high fold coverage (many overlapping recordings of the same subsurface point), which improves signal-to-noise ratio and image quality.

3D data enables far more accurate imaging of complex geology. 3D migration algorithms (Kirchhoff migration, wave-equation migration) reposition reflections to their true subsurface locations, correcting for the distortions that occur when structures are steeply dipping or laterally variable.

Three key interpretation techniques take advantage of the 3D volume:

1. Volume rendering: Displays the entire seismic cube as a semi-transparent volume, letting you visualize spatial relationships between faults, channels, and other features from any angle.
2. Horizon slicing: Extracts seismic attributes along a picked horizon surface, mapping lateral variations in reservoir properties across the field.
3. Attribute analysis on volumes: Applies coherence, curvature, or amplitude attributes across the full 3D dataset to highlight faults, channel systems, or fluid contacts that might not be visible on individual sections.

4D Seismic Monitoring and Reservoir Management

4D seismic (time-lapse seismic) repeats a 3D survey over the same area at different times during production, typically at intervals of several years. The surveys must be processed carefully to ensure consistency, so that any differences between them reflect actual reservoir changes rather than acquisition artifacts.

4D seismic detects production-related changes in the reservoir:

• Fluid substitution: As water replaces oil or gas during production, acoustic impedance changes, producing a detectable shift in seismic response.
• Pressure depletion: Reduced pore pressure can cause compaction and velocity changes, altering travel times.
• Bypassed reserves: Areas that haven't been swept by the production process show up as unchanged zones, guiding decisions about where to drill infill wells.

Two primary interpretation techniques are used with 4D data:

1. Difference volumes: Subtract the baseline (pre-production) survey from the monitor (later) survey. Any non-zero signal highlights where the reservoir has changed.
2. Time-shift analysis: Measures travel-time differences between baseline and monitor surveys. These shifts indicate velocity changes caused by pressure depletion or fluid substitution.

Integrating 3D and 4D seismic with reservoir simulation models closes the loop between observation and prediction. Simulation models use seismic-derived properties (porosity, permeability) to predict fluid flow and production. When 4D seismic reveals discrepancies between predicted and actual reservoir behavior, the model can be calibrated and updated. This integrated workflow, combining seismic, well, and production data, enables data-driven decisions about well placement, injection strategy, and overall field management.