Skip to content
Brent$117.6Gas TTF$1,192USD453.87 ₸
MunaiHub
GeoscienceCorrelation and maps Intermediate 23 min read

Well correlation

Well correlation is the matching of the sections of wells. It means finding which layers, members (groups of beds) and surfaces in one well correspond to which in another. Every line on a correlation panel (a drawing of well sections side by side) is a hypothesis. It says either that a layer continues unbroken between the wells (lithostratigraphic correlation), or that two points of the sections are of the same age (chronostratigraphic correlation). We correlate from logs (measurements made in the well), core (rock cylinders cut from the well), cuttings (rock chips brought up by the drilling fluid), fossils and marker beds. We check the result against seismic data, pressures and the behaviour of wells in production.

Why it matters

Almost everything that we build afterwards depends on correlation. This includes the picks for structure and thickness maps (the depths chosen for the tops and bases of layers), the identification of layers and development units (groups of layers developed as one), the layers of the geological model, the expected connection between injection and production wells, and the perforation intervals. A correlation error can turn one layer into two, or two into one. It can join sand bodies that do not communicate, or hide a fault. Then injection produces no response, water breaks through early, and the hydrocarbons in place are calculated wrongly.

Quick refresher

  • Lithological boundaries are diachronous (the same rock type is older in one place and younger in another). Time lines never cross one another. Examples of time lines are marker beds, biostratigraphic events, ash beds and flooding surfaces. “Sand-to-sand” correlation can join bodies of different ages that do not communicate.
  • A good marker bed was deposited quickly and everywhere, does not depend on facies (the type of depositional environment), and is easy to recognise on the logs. Examples are ash, shales of the maximum flooding surface (the time of the deepest water), and limestones and anhydrites that stay the same over long distances. Channel sandstones, reefs and coals of river floodplains are poor markers.
  • A structural cross-section in depths below sea level shows how the layers lie today. A section flattened on a marker removes later deformation and shows how the layers related to one another when the marker was deposited. If an interval thins towards the crest, the structure was growing while it was deposited.
  • We correlate deviated wells only in true vertical depths, and we compare thicknesses as TVT or TST. A hole drilled against the dip “squeezes” an interval in true vertical depth, and one drilled down the dip “stretches” it.
  • The line of correlation is drawn on a plot of “depth of an event in one well against its depth in the other well”. Its slope gives the ratio of the rates at which sediment accumulated. A horizontal step on the line is a hiatus in the section of the well on the vertical axis.
  • A missing or repeated interval, converging lines and a change in the pattern of the logs can have different causes: a fault, an unconformity, a pinch-out, a facies change or the geometry of the hole. We tell them apart by several signs. If it is doubtful whether a body continues between wells, we keep both versions and test them against data.
Key relationships
h=Δlcos⁡γ,hTVT=hcos⁡δ,ΔH=Δlcos⁡θh = \Delta l\cos\gamma,\qquad h_{\text{TVT}} = \frac{h}{\cos\delta},\qquad \Delta H = \Delta l\cos\theta

true stratigraphic thickness (TST), true vertical thickness (TVT) and thickness in true vertical depth of an interval of length Δl\Delta l along a deviated well

zB=zˉB+b (zA−zˉA),b=∑(zA−zˉA)(zB−zˉB)∑(zA−zˉA)2z_{\text{B}} = \bar z_{\text{B}} + b\,(z_{\text{A}} - \bar z_{\text{A}}),\qquad b = \frac{\sum (z_{\text{A}} - \bar z_{\text{A}})(z_{\text{B}} - \bar z_{\text{B}})}{\sum (z_{\text{A}} - \bar z_{\text{A}})^2}

the line of correlation fitted by least squares (the line closest to all the points)

P=1−d sin⁡αWP = 1 - \frac{d\,\sin\alpha}{W}

probability that a body of width WW penetrated by one well reaches another well at a distance dd

Common pitfalls
  • Correlating on the absolute values of a log instead of the pattern of the curves, and comparing the logs of deviated wells by measured depth (the length of hole drilled).
  • Taking the top or bottom of a log run, a casing shoe (the bottom end of a steel casing string), a splice between runs or jumps in the sonic log for geological boundaries.
  • Drawing correlation lines through an unconformity as through a continuous layer.
  • Taking the apparent thinning of a layer in a deviated well for a missing interval.
  • Taking the bases of fossil ranges from cuttings. Caving (pieces falling from the borehole wall) shifts them downwards.

After this lesson you can

  • list the data used for correlation and prepare them: true vertical depths below sea level, uniform scales, log normalisation (bringing the logs of different wells to a common scale) and the rejection of technical artefacts (false features caused by the tool or the hole)
  • choose marker beds (distinctive layers that can be recognised in many wells) and a datum (the surface used as a reference), build a structural cross-section and a cross-section flattened on a marker, and read from them how a structure grew
  • explain how lithostratigraphic correlation (matching rock types) differs from chronostratigraphic correlation (matching time), and what errors of connectivity “sand-to-sand” correlation leads to
  • convert intervals measured along deviated (slanted) wells into true vertical thickness (TVT) and true stratigraphic thickness (TST)
  • construct the line of correlation of two sections from biostratigraphic events (the first or last appearance of a fossil species), and use it to find hiatuses (gaps in the record) and the expected depths of layers
  • recognise correlation pitfalls (faults, unconformities, facies changes and pinch-outs) and choose between the “layer-cake” and “jigsaw” models

What a correlation line asserts

The section of each well is a one-dimensional record: a sequence of rocks and log curves against depth. Correlation ties these records into a three-dimensional picture. On a correlation panel we place the well sections side by side and join equal boundaries with lines. In Russian-language practice, people distinguish general correlation from detailed correlation. General correlation matches formations, horizons and large members through the whole section. Detailed correlation works layer by layer and bed by bed in the productive interval. The result of detailed correlation is an agreed table of picks and a panel, from which we build maps and models.

The lines on a panel are not observations. They are assertions. When you join the top of a sandstone in two wells, you assert that the same sandstone lies between them: it does not pinch out (thin out to nothing), it is not cut by a fault, and it is not eroded. When you join two ash beds, you assert that they fell at the same time. These two kinds of assertion are fundamentally different, and a later section deals with them. A new well, seismic data or a response to injection can refute any of them, so we revise the correlation throughout the life of a field.

What correlation gives
Result
Agreed picks of the top and base of layers
Where it is used
structure and thickness maps (Structure maps and contouring)
Result
Division of the productive interval into layers and thin beds
Where it is used
development units, layers of the geological model
Result
Continuity of reservoirs and shale partings (thin shale layers that separate reservoir beds) between wells
Where it is used
expected connection between injection and production wells (Why water is injected: pressure maintenance and oil displacement)
Result
Missing and repeated intervals
Where it is used
faults and unconformities (Faults)
Result
Position of the bit relative to the target layer
Where it is used
steering the well, taking core, running casing

Data for correlation and how to prepare them

Data for correlation
Source
Gamma ray log (GR, Gamma ray log and natural gamma ray spectrometry), spontaneous potential (SP, Spontaneous potential (SP) log)
What it gives
lithology and shale content; the pattern of the curves repeats in neighbouring wells
Limitations
absolute values depend on the tool, the hole and the drilling fluid (mud)
Source
Resistivity logs (Resistivity logs: electrical and induction tools)
What it gives
high-resistivity markers — limestones, anhydrites, coals
Limitations
depends on the fluid saturation: an oil–water contact looks like a change of rock type
Source
Sonic and density logs (Sonic log, Density log and photoelectric factor (Pe))
What it gives
dense and soft beds, coals; tie to seismic data
Limitations
cycle skipping (false jumps in the sonic log), washouts (places where the hole is washed wider than the bit)
Source
Dipmeter and borehole image logs
What it gives
dip angles, which show unconformities, faults and drag near faults (layers bent beside a fault)
Limitations
need conversion and checking against core
Source
Core, cuttings, gas logging
What it gives
rock, textures, fossils, signs of the depositional environment, shows (traces of oil or gas)
Limitations
cuttings are mixed over the interval and arrive late; there is little core
Source
Biostratigraphy
What it gives
levels of equal age: the first and last appearances of species
Limitations
from cuttings only the tops of ranges are reliable
Source
Seismic data
What it gives
tracing horizons between wells, unconformities, faults
Limitations
resolution is tens of metres; the pinch-out point of thin layers is not visible
Source
Reservoir pressures and contacts
What it gives
whether layers communicate: a common pressure gradient, a common contact
Limitations
pressures change as the field is developed

The more independent kinds of data agree with the correlation, the more reliable it is. Pressures from the wireline formation tester (a tool that measures formation pressure through the borehole wall) are covered in the topic Pressure measurement with the wireline formation tester.

Before correlation we bring the data to a common form. Picks are taken in measured depth, but we compare sections in true vertical depths below sea level (TVDSS; depth references). For deviated wells we calculate these from the directional survey (the measured path of the well), and we compare thicknesses as true vertical thicknesses or true stratigraphic thicknesses (see the section on deviated wells). We plot all logs at one vertical scale and with the same curve scales.

The absolute values of logs in different wells cannot be compared directly. Tools of different years and contractors, hole diameters and drilling fluids (muds) shift the scale (when GR misleads). So we correlate on the pattern of the curves. For quantitative comparisons we normalise the curves. In each well we take the same reference beds: a persistent shale member and a clean sandstone or limestone. Then we convert the scale linearly so that the readings in these beds coincide with those in the reference well.

Normalisation of the gamma ray log to a reference well
GRnorm=GRclref+IGR(GRshref−GRclref),IGR=GR−GRclGRsh−GRclGR_{\text{norm}} = GR_{cl}^{\text{ref}} + I_{GR}\left(GR_{sh}^{\text{ref}} - GR_{cl}^{\text{ref}}\right),\qquad I_{GR} = \frac{GR - GR_{cl}}{GR_{sh} - GR_{cl}}
where:
  • GRGRGR reading in the well being normalised, API
  • GRnormGR_{\text{norm}}the same reading on the scale of the reference well, API
  • GRcl, GRshGR_{cl},\ GR_{sh}readings against the reference clean bed and the reference shale member in the well being normalised, API
  • GRclref, GRshrefGR_{cl}^{\text{ref}},\ GR_{sh}^{\text{ref}}readings against the same beds in the reference well, API
  • IGRI_{GR}gamma ray index (what we find from the GR), —

Units: API units; the index is dimensionless

Holds when
  • The reference beds are present in all wells and have the same composition. We check this against core and regional data.
  • The differences between the scales are linear (a shift and a stretch). We correct for non-linear distortions, such as washouts and barite (a heavy mineral added to the drilling fluid), before normalisation.

Normalisation removes the differences between tools and well conditions, but it also removes any real changes in the reference beds. So we choose reference beds that are persistent and verified. We use the normalised curves for calculating shale volume and for automatic correlation, not in place of geological sense.

  • The top and the bottom of a log run are the limits of the measurement, not of the rocks. We do not join them with neighbouring wells.
  • A casing shoe, the passage of the record from open hole (the part of the well with no casing) into casing, and the splice of two runs give steps on the logs that look like geological boundaries.
  • Jumps in the sonic log and washouts seen on the caliper log (a log of the hole diameter) distort the curves in individual intervals (Borehole environment: invaded zone and mudcake).
  • The depths of different runs and of core are offset from one another by metres. Before correlating, we match them on the GR.

Marker beds, lithostratigraphic and chronostratigraphic correlation

A marker bed is a thin, persistent layer or surface that we can recognise in every well and that was deposited at practically the same time over the whole area. Marker beds are the framework of correlation.

  • Good markers: volcanic ash beds (a narrow peak on the GR), shales of maximum flooding surfaces (a GR maximum, often because of uranium), limestones and anhydrites that stay the same over long distances (a GR minimum, high resistivity and density), and coals of a delta plain that accumulated at the same time over a large area.
  • Poor markers: channel and delta sandstones, reef buildups, coals of river floodplains and fan deposits. These are bodies of limited extent, and their position in the section changes from well to well. Shale members that look alike are poor markers too, because they are easily confused.

Lithostratigraphic correlation joins like rocks: the top of a sandstone with the top of a sandstone. Its lines are the boundaries of rock bodies, and they are diachronous: when the shoreline shifts, the same sand member is younger in one place and older in another (lithostratigraphy). Chronostratigraphic correlation joins surfaces of equal age: marker beds, biostratigraphic events, flooding surfaces (surfaces made when the sea suddenly deepened) and sequence boundaries (sequence stratigraphy). Time lines never cross. Lithological boundaries may cross them at any angle.

Lithostratigraphic and chronostratigraphic correlation of the same wellsLithostratigraphic: “sand-to-sand”1234Chronostratigraphic: along flooding surfaces1234← landbasin →sand — yellow, shales — greyflooding surfaces — time lines

Two correlation panels of the same four wells, land on the left, basin on the right. In each well the sand intervals are yellow and the shales grey; in the left-hand wells a thick sand lies lower in the section, in the right-hand wells higher. On the upper panel the thickest sand interval of each well is joined to the thickest one in the neighbouring well, giving a single band of sand rising from land to basin. On the lower panel four gentle dashed lines, slightly inclined towards the basin, are drawn — flooding surfaces; beneath each of the upper three lies a sand wedge thinning in both directions, and each successively younger wedge is shifted to the right. The band of the upper panel crosses these time lines.

SchematicThe same four wells in a prograding succession (a succession built outwards as the shoreline moves towards the sea), correlated in two ways (schematic). Top: “sand-to-sand”, which gives one continuous layer. Bottom: along flooding surfaces (time lines), which gives three sand layers shifted towards the basin and separated by shales.

The difference is not just academic. In a prograding succession the shoreline moves steadily seawards. Each new sand layer is separated from the previous one by a thin shale deposited at the next flooding. “Sand-to-sand” correlation joins them into one connected layer. Correlation along flooding surfaces shows several layers that overlap only partly. Flooding shales persist over large areas and often act as barriers to crossflow. Water injected into one layer may never reach a producing well where the “sand” lies at the same level but belongs to another layer.

Biostratigraphic events in wells are the levels at which species appear and disappear (biostratigraphy). When we drill downwards, the first cuttings sample in which a species is found gives the top of its range, and this is reliable. Fragments that have caved from the borehole wall higher up the hole push the lower limit of the range downwards. So we take the lower limit from core and sidewall core (core cut from the borehole wall) (fossils in cuttings and core). Species that lived only in a particular environment appear and disappear together with the rock, not with time, so they are no use for correlation in time.

Order of correlation
  1. 1
    Prepare the data

    Use true vertical depths below sea level and uniform scales, normalise the curves, and reject technical artefacts.

  2. 2
    Build the framework

    Start with the most reliable and largest boundaries: regional shale members, ash beds, flooding surfaces and formation boundaries. Correlate them in all wells.

  3. 3
    Match the layers

    Within the framework, match layers and thin beds, from large to small, with the top and base of each. Show a pinch-out by converging lines.

  4. 4
    Check for realism

    Lines between neighbouring wells should be roughly parallel or diverge gradually, and they must not cross. Every change in thickness and attitude must be explained by a sequence of events.

  5. 5
    Check against independent data

    Seismic data, pressures and contacts, the behaviour of wells in production.

The datum: structural and flattened cross-sections

We can arrange the same wells on a panel in two ways. On a structural cross-section each well stands at its own depths below sea level. The section shows how the layers lie today: crests, troughs and faults. On a flattened cross-section (also called a stratigraphic cross-section; the operation is called flattening, or hanging) we shift all the wells vertically so that the chosen marker lies horizontal. This marker is the stratigraphic datum. If the marker was deposited horizontally, such a section removes all later deformation and shows the layers as they were when the marker was deposited.

Horizontally, we place the wells either at true distances or at equal intervals. With true distances we get a cross-section, in which shapes and slopes are visible. With equal intervals we get a correlation panel, on which it is easier to compare the sections themselves. A network of panels that intersect at wells and are drawn in perspective is a fence diagram.

  • A good datum is close to the interval of interest and above it. The farther the marker, the more later events enter the panel.
  • A good datum is present in all wells. It is not eroded and not cut out by a fault between them.
  • A good datum was deposited horizontally. Examples are flooding shales, ash beds and persistent coals and limestones. The top of a sand body, a reef or inclined progradational layers, and an unconformity surface with relief, are not good datums.
Markers M1 and M3 in the wells of the panel
Well
A
Distance from A, m
0
M1, m
2214.6
M3, m
2299.1
Well
B
Distance from A, m
1100
M1, m
2187.3
M3, m
2262.5
Well
C
Distance from A, m
2000
M1, m
2179.8
M3, m
2251.7
Well
D
Distance from A, m
3300
M1, m
2206.1
M3, m
2287.8

Depths below sea level. M1 is the shale of the maximum flooding surface (GR peak), M3 is an ash bed (narrow GR peak). The data are assumed for illustration.

Worked exampleGrowth of a structure from a section flattened on a marker
Given (assumed data)
Wells A, B, C, D
on one line across the structure; the depths of markers M1 and M3 are in the table above
Layers
almost horizontal (dip 1–2°), so thicknesses in vertical wells equal the true vertical thicknesses and are almost equal to the true ones
Solution
  1. Vertical thickness of the interval M1–M3: A — 2299.1−2214.6=84.52299.1 - 2214.6 = 84.5 m, B — 75.2 m, C — 71.9 m, D — 81.7 m.
  2. Flatten on M1. On the panel, M3 then stands in each well at a depth equal to this thickness. At C it is 84.5−71.9=12.684.5 - 71.9 = 12.6 m higher than at A. So the interval over the present-day crest is 15% thinner.
  3. The present-day difference in depth between A and C is 34.8 m for M1 and 47.4 m for M3.
  4. The difference 47.4−34.8=12.647.4 - 34.8 = 12.6 m is the growth of the structure in the time between the deposition of M3 and M1. It is 27% of the present-day amplitude at M3.
Answer: The structure was growing while the interval M1–M3 accumulated. Towards the crest the interval is thinner by 12.6 m, and about 27% of the present-day amplitude at M3 had already formed before M1 was deposited.

What it means. A section flattened on M1 shows the relief of M3 at the time M1 was deposited. This holds if M1 was deposited horizontally and if we can neglect differences in compaction. Sandstones compact less than shales, so if the wells have a different share of sand, part of the difference in thickness is not growth but compaction. The early growth of a structure matters for the timing of trap formation relative to migration (the movement of oil and gas into the trap) (timing of formation). It matters for the reservoir too: over a growing crest the layers are thinner and may differ in composition. The data in this example are assumed for illustration.

Wells on a structural cross-section and on a section flattened on marker M1Structural cross-section: depths below sea level2180220022202240226022802300ABCDM1M2M3IIIdistances to scale: A–D 3.3 kmFlattened on M1: equal spacing020406080ABCDM1M2M3IIIm below M1curves — GR, 0–150 API; sands in yellow

Two panels with four columns of gamma ray curves. In the upper one the columns stand at different heights: wells B and C are higher, A and D lower — an anticline; the lines of markers M1, M2 and M3 and of the boundaries of sand layers I and II are bent in an arc, the lower lines more strongly than the upper. In the lower panel the columns are aligned so that the M1 line is horizontal; the M3 line at wells B and C rises higher than at A and D — the interval over the crest is thinner. Sand II is shaded yellow between A, B and C and converges to a point between C and D, where shales replace it.

CalculatedWells A–D on a structural cross-section (top: depths below sea level, distances to scale) and on a section flattened on marker M1 (bottom: equal spacing). The curves are the GR, and yellow marks sand layers I and II. Between C and D, shales replace sand II, so the lines of its top and base converge. The lower panel shows that the interval M1–M3 is thinner over the crest.

Deviated and horizontal wells: TVD, TVT and TST

The log of a deviated well is recorded along the hole. For correlation we first convert it to true vertical depth (TVD). Otherwise all the intervals come out stretched and the markers come out deeper than they are. But even in true vertical depths a deviated well sees thicknesses differently from a vertical one. As it crosses a layer it moves sideways, and the top of the layer is already at a different depth there.

We describe the thickness of a layer by two quantities (true thickness and vertical thickness). The true stratigraphic thickness hh (TST) is measured normal (at right angles) to the bedding. The true vertical thickness hTVTh_\text{TVT} (TVT) is measured vertically at the given point. In geological texts they are often called simply the true thickness and the vertical thickness. In a vertical well the difference in depth between the top and the base equals the TVT. In a deviated well, the difference in true vertical depth ΔH\Delta H (simply called the TVD thickness) equals the TVT only in the special case when the hole runs along the strike of the layer.

Thicknesses of an interval in a deviated well
h=Δlcos⁡γ,hTVT=hcos⁡δ,ΔH=Δlcos⁡θh = \Delta l\cos\gamma,\qquad h_{\text{TVT}} = \frac{h}{\cos\delta},\qquad \Delta H = \Delta l\cos\theta
where:
  • Δl\Delta llength of hole between the top and the base of the interval (difference in measured depth), m
  • γ\gammaangle between the hole and the normal to the bedding: cos⁡γ=∣sin⁡δ sin⁡θ cos⁡(a−φ)−cos⁡δ cos⁡θ∣\cos\gamma = |\sin\delta\,\sin\theta\,\cos(a - \varphi) - \cos\delta\,\cos\theta|, °
  • δ, a\delta,\ adip and dip azimuth of the layers, °
  • θ, φ\theta,\ \varphiinclination and azimuth of the hole over the interval, °
  • h, hTVTh,\ h_{\text{TVT}}true stratigraphic thickness (TST) and true vertical thickness (TVT), m
  • ΔH\Delta Hdifference in true vertical depth between the top and the base in this well, m

Units: lengths in the same units; angles in degrees

Holds when
  • The layer is planar and its top and base are parallel. The hole is straight over the interval; otherwise we divide the interval into parts.
  • The azimuths of the layers and of the hole are measured from the same north. The formula for cos⁡γ\cos\gamma comes from the topic vector calculation.

The ratio is ΔH/hTVT=cos⁡θcos⁡δ/cos⁡γ\Delta H/h_\text{TVT} = \cos\theta\cos\delta/\cos\gamma. A hole running against the dip, up-dip, “squeezes” the interval in true vertical depth, and one running down the dip “stretches” it. The equality ΔH=hTVT\Delta H = h_\text{TVT} holds only for a vertical well and for a hole running along the strike of the layer. This is why we compare deviated wells with vertical ones using logs converted into true vertical thickness (TVT), and we compare thicknesses in areas of differing dip using true thickness (TST).

Worked exampleA layer “thinned” in a deviated well
Given (assumed data)
Layer (from image logs)
dip azimuth 150°, dip 12°
Vertical well 31
top and base of the layer at 2318.4 and 2351.1 m; vertical thickness 32.7 m
Nearby deviated well 32
over the interval inclination 55°, azimuth 315°; top and base at 2841.6 and 2885.7 m measured depth
Solution
  1. Length of hole in the layer Δl=2885.7−2841.6=44.1\Delta l = 2885.7 - 2841.6 = 44.1 m; in true vertical depth ΔH=44.1⋅cos⁡55°=25.3\Delta H = 44.1\cdot\cos 55° = 25.3 m — 7.4 m less than in well 31.
  2. Angle between the hole and the normal to the layer: cos⁡γ=∣p−q∣\cos\gamma = |p - q|, where p=sin⁡12° sin⁡55° cos⁡(150°−315°)=−0.16451p = \sin 12°\,\sin 55°\,\cos(150° - 315°) = -0.16451 and q=cos⁡12° cos⁡55°=0.56104q = \cos 12°\,\cos 55° = 0.56104; cos⁡γ=0.72555\cos\gamma = 0.72555, γ=43.5°\gamma = 43.5°.
  3. True stratigraphic thickness h=44.1⋅0.72555=32.0h = 44.1\cdot 0.72555 = 32.0 m; true vertical thickness hTVT=32.0/cos⁡12°=32.7h_\text{TVT} = 32.0/\cos 12° = 32.7 m — the same as in well 31.
Answer: The thickness of the layer is the same in both wells: 32.7 m vertical, 32.0 m true. The “thinning” to 25.3 m is a consequence of the geometry of the hole.

What it means. Without conversion, correlation would show that about 7.4 m of section is missing in well 32, as though the hole had crossed a normal fault, or that the layer is thinning. Here the hole runs almost against the dip — the azimuths differ by 165°; the steeper the dip and the greater the inclination, the stronger the “squeezing”. A hole running down the dip, by contrast, “stretches” the interval. The data in this example are assumed for illustration.

Thicknesses of a layer in a vertical and a deviated welllayerverticalhTVTdeviated, against the dipΔlΔHhschematic; dip to the right, ΔH < hTVT

Vertical section: a sandstone layer dipping down to the right. On the left a vertical well crosses the layer; the vertical segment between the top and the base is labelled h TVT. Further right a deviated well runs up the dip — from right to left going down — and crosses the layer obliquely; the length of hole inside the layer is labelled delta l, and its projection onto the vertical is labelled delta H, noticeably shorter than h TVT. Between the top and the base, perpendicular to the layers, the true thickness h is shown.

SchematicThe same layer in a vertical and a deviated well (schematic, section along the dip). The vertical well cuts the true vertical thickness hTVTh_\text{TVT}. The deviated well runs against the dip. It crosses the layer over a length of hole Δl\Delta l, but it descends vertically by only ΔH<hTVT\Delta H < h_\text{TVT}. The true thickness hh is measured at right angles to the layers.

A horizontal hole runs almost along the layer and may cross the same boundaries repeatedly. In true vertical depths this looks like a repetition of the section, even though there is no fault (a fault in a well). We describe the position of such a hole by its distance from a chosen boundary, measured along the normal. This distance is the stratigraphic depth. We compare the measured log with one calculated from a model of the layer. This is how geosteering (steering the well by the rock it meets) is done (Horizontal wells).

Graphic correlation: the line of correlation

Graphic correlation is a quantitative way to match two sections by events that date the layers. We plot the depths of the same events in the two wells: the reference well on the horizontal axis, and the well being matched on the vertical axis. The events are tops and bases of species ranges, ash beds, GR peaks of flooding surfaces, and reversals of magnetic polarity. If both wells accumulated sediment without breaks, the points fall near a straight line, the line of correlation.

  • The slope of the line is the ratio of accumulation rates. A slope of 1.3 means that in the same time 1.3 times as much rock was deposited in the well being matched (or the rock there was compacted less).
  • A horizontal step means that the depth in the well being matched does not change while the depth in the reference well does. Part of the section of the reference well is absent in the matched well: a hiatus with erosion, or a sharply reduced interval. A vertical step means the same for the reference well.
  • A kink (a sharp bend in the line) is a change in the ratio of accumulation rates.
  • Points off the line are events whose position is distorted. Examples are range bases from cuttings that caving has shifted downwards, and species that depended on the environment or arrived in the area later.
Least-squares line of correlation
zB=zˉB+b (zA−zˉA),b=∑k(zA,k−zˉA)(zB,k−zˉB)∑k(zA,k−zˉA)2z_{\text{B}} = \bar z_{\text{B}} + b\,(z_{\text{A}} - \bar z_{\text{A}}),\qquad b = \frac{\sum_k (z_{\text{A},k} - \bar z_{\text{A}})(z_{\text{B},k} - \bar z_{\text{B}})}{\sum_k (z_{\text{A},k} - \bar z_{\text{A}})^2}
where:
  • zA,k, zB,kz_{\text{A},k},\ z_{\text{B},k}depths of the kk-th event in the reference well A and in the well being matched, B, m
  • zˉA, zˉB\bar z_{\text{A}},\ \bar z_{\text{B}}mean depths of the selected events, m
  • bbslope of the line — the ratio of accumulation rates in B and in A, —

Units: metres; the same depth reference in both wells (true vertical)

Holds when
  • The events are of the same age in both wells. The fit uses reliable events from one stretch without breaks.
  • The ratio of accumulation rates is constant over the stretch. Where it changes or there is a hiatus, we build the line in parts.
  • The errors in the depths of events in B are larger than in A (for example, cuttings against core), so we fit zBz_{\text{B}} against zAz_{\text{A}}.

With the line we can convert any depth in one well into a depth in the other. We can find where to expect in well B the top of a layer that is known in A, or which interval of A corresponds to a hiatus in B. This works for log markers too: GR peaks, limestones and ash beds are events of the same kind.

Events in reference well A (core) and in well B (cuttings and log)
Event
species a
Type
top of range
A, m
2104.0
B, m
2230
Event
species b
Type
top of range
A, m
2121.5
B, m
2250
Event
ash
Type
log marker
A, m
2138.0
B, m
2263.8
Event
species c
Type
top of range
A, m
2152.0
B, m
2280
Event
species d
Type
top of range
A, m
2168.5
B, m
2290
Event
species e
Type
top of range
A, m
2183.0
B, m
2290
Event
species f
Type
top of range
A, m
2201.0
B, m
2305
Event
MFS
Type
log marker
A, m
2214.0
B, m
2318.0
Event
species g
Type
top of range
A, m
2229.5
B, m
2340
Event
species f
Type
base of range
A, m
2240.0
B, m
2375
Event
species h
Type
top of range
A, m
2252.0
B, m
2370

Tops and bases of ranges in well B come from cuttings sampled every 5 m. We record each one in the first sample in which the species is found (top) or in the sample where it is found for the last time (base). The ash bed and the GR maximum of the maximum flooding surface (MFS) come from the log. The data are assumed for illustration.

Worked exampleLine of correlation, hiatus and prediction of the top of a layer
Given (assumed data)
Events
11 events from the table above
Hiatus in well B
from the log and cuttings — at a depth of 2288 m
Top of the pay layer in well A
2205.6 m
Solution
  1. Take the lower group of reliable events (tops of species f, g, h and MFS): zˉA=2224.125\bar z_{\text{A}} = 2224.125 m and zˉB=2333.25\bar z_{\text{B}} = 2333.25 m. The sums are ∑(zA−zˉA)(zB−zˉB)=1868.4\sum(z_{\text{A}} - \bar z_{\text{A}})(z_{\text{B}} - \bar z_{\text{B}}) = 1868.4 and ∑(zA−zˉA)2=1443.2\sum(z_{\text{A}} - \bar z_{\text{A}})^2 = 1443.2 m². The slope is b=1.295b = 1.295.
  2. The base of species f in well B (2375 m) lies 21 m below the line, because fragments caved from the borehole wall have “extended” the species downwards. We leave this event out of the fit.
  3. The upper group (tops a, b, c and the ash bed) gives a slope of 1.020. The tops of species d and e in well B were found in a single sample, at 2290 m, just below the hiatus. So the points lie on a horizontal step.
  4. The hiatus depth of 2288 m in well B corresponds in well A to a depth of 2160 m on the upper line and 2189 m on the lower one. So about 29 m of the section of well A is missing in well B.
  5. Prediction: the top of the layer is expected in well B at 2333.25+1.295⋅(2205.6−2224.125)=2309.32333.25 + 1.295\cdot(2205.6 - 2224.125) = 2309.3 m.
Answer: Below the hiatus the section of well B is 1.295 times as thick as that of well A. At the hiatus about 29 m of section is missing in well B. We expect the top of the pay layer in well B at about 2309 m.

What it means. Tops of species from cuttings are recorded in a sample 0–5 m below the true top. So the slope of a line fitted to such points is uncertain by a few per cent. This is why events tied accurately by log, the ash bed and MFS, are so valuable. The hiatus at 2288 m in well B would cut across direct “layer-to-layer” correlation. There is no point in looking in well B for those layers of well A that lie in the missing interval. We check the predicted top with the log while drilling. The data in this example are assumed for illustration.

Graphic correlation of two wells2100214021802220226022202260230023402380Depth of event in well A (core), mDepth in well B, mhiatus in well Bslope 1.02slope 1.29abashcdefMFSgf — cavedhspecies topspecies baselog marker

Plot: the horizontal axis is the depth of the event in well A from 2095 to 2265 metres, the vertical axis is the depth in well B from 2220 to 2385 metres, depth increasing downwards. Four points of the upper group lie on a straight line with a slope of about 1. Then two points at the same depth of 2290 metres form a horizontal step. Below, four points of the lower group lie on a steeper straight line with a slope of about 1.3. One point — the base of species f — lies about 21 metres below the lower line and is labelled as shifted by caving.

CalculatedGraphic correlation of wells A and B from the data of the example. Downward triangles are tops of species ranges, the upward triangle is a base, and diamonds are events on the log. The figure shows two lines of correlation and the horizontal step between them, which is the hiatus in well B at 2288 m. The base of species f, shifted by caving, lies far from the line.

Correlation pitfalls: faults, unconformities, facies changes, pinch-outs

What is seen on the panel and what it may be
What is seen
Part of the section is missing
Possible causes
normal fault; unconformity with erosion; pinch-out; geometry of a deviated hole
How to tell
at a normal fault, the missing interval is nearly constant along the segment of the fault; at an unconformity, it grows towards the crest and cuts out different layers; a pinch-out is a gradual thinning; for a hole, check the TVT
What is seen
Part of the section is repeated
Possible causes
reverse fault or thrust; overturned limb (a fold limb tilted past the vertical); horizontal hole
How to tell
in an overturned limb the sequence is reversed; for a hole, use the directional survey
What is seen
An interval is sharply thicker
Possible causes
growth of a fault during deposition; channel incision (a river channel cut into the layers below); repetition
How to tell
growth shows as thickening in the hanging wall in one interval; incision shows as a sandstone with an eroded base
What is seen
The pattern of the curves in a member has changed
Possible causes
facies change; a different fluid; hole conditions
How to tell
a facies change keeps the markers above and below; a different fluid changes only the resistivity; washouts show on the caliper log

How to estimate the throw of a normal fault from a missing interval, and how repetition differs from the other three cases, is covered in the topic a fault in a well.

An unconformity. This is a surface where part of the rock record is missing. Beneath an erosion surface the correlation lines run up to it and end: in neighbouring wells different layers lie below the unconformity. We must not join layers across an unconformity as if they were continuous. We recognise an unconformity by the truncation of markers (they are cut off at the surface), by a sharp change of dip on the dipmeter, by the absence of biostratigraphic zones, and by weathering crusts and conglomerates at the boundary (unconformities).

Facies change and pinch-out. In a facies change the layer is still there, but the rock in it is different: a sandstone passes sideways into shale. We correlate the surfaces that bound the layer, not the sandstone, and the panel shows where the rock changes. In a pinch-out the layer thins to zero between the wells, and on the panel the lines of its top and base are brought together at a point. The wells cannot say where exactly this happens between them. Seismic data sometimes give a hint, but thin layers and the pinch-out point itself are usually below its resolution. The vertical sequences of rocks, from which we recognise the environment and the shape of bodies, and the characteristic shapes of GR curves are covered in the topic Facies, facies successions and sedimentary logs.

The last check is for realism. Correlation lines between neighbouring wells are roughly parallel or diverge gradually. They do not cross. Every change in thickness and attitude is explained by an understandable sequence of events: deposition, growth of the structure, erosion and faulting. A panel on which lines converge and diverge without any system is almost certainly wrong.

Choosing between models: “layer-cake” and “jigsaw”

How we correlate sand bodies depends on how large they are compared with the distance between wells. If the bodies extend for kilometres, we can trace the layers from well to well, and we build the section as a “layer-cake”. Examples are shoreface and shelf layers (sands deposited near the shore), turbidite lobes (deposits of underwater sediment flows) and platform carbonates. If the bodies are comparable to the well spacing (channels and bars of river and delta systems), they abut one another with almost no gaps. But some contacts let fluid through, and others are separated by thin partings of low permeability. This is a “jigsaw”. If the bodies are much smaller than the well spacing, we cannot match them one by one at all (deterministically). Such sections are described statistically (Reservoir architecture and performance).

Probability that a body reaches the neighbouring well
P=1−d sin⁡αW(P≥0)P = 1 - \frac{d\,\sin\alpha}{W}\quad (P \ge 0)
where:
  • PPprobability that a body penetrated by one well is also penetrated by the other, —
  • dddistance between the wells, m
  • WWwidth of the body (channel, channel belt, sand layer) across its axis, m
  • α\alphaangle between the axis of the body and the line joining the wells, °

Units: any consistent units of length

Holds when
  • The body is straight and of constant width, and its position relative to the wells is random.
  • The formula becomes crude at small α\alpha, when the body is elongated almost along the line of the wells. Then the sinuosity (how much the body winds) and the length of the body decide, not its width.
  • When dsin⁡α>Wd\sin\alpha > W, the probability is zero, because the body is narrower than the distance between the wells measured across its axis.

The width of the body along the line of the wells is W/sin⁡αW/\sin\alpha. If one well is in the body, the second well, a distance dd away, falls in the same body with probability 1−dsin⁡α/W1 - d\sin\alpha/W. The ratio of the size of the bodies to the well spacing decides whether we have a layer-cake or a jigsaw.

Worked exampleWhether to join the sandstones of two wells
Given (assumed data)
Distance between wells 1 and 2
600 m
Option 1: a belt of river channels
width about 800 m (from analogues), axis at 60° to the line of the wells
Option 2: a shoreface sand layer
extent across the shoreline about 6 km, the line of the wells across the shore
Solution
  1. Width of the belt along the line of the wells: 800/sin⁡60°=924800/\sin 60° = 924 m.
  2. Probability that a belt penetrated by well 1 is also penetrated by well 2: P=1−600/924=0.35P = 1 - 600/924 = 0.35. If the belt ran across the line of the wells (α=90°\alpha = 90°): 1−600/800=0.251 - 600/800 = 0.25.
  3. For a well 1200 m away on the same line, 1200>9241200 > 924, so the probability is zero.
  4. Shoreface layer: P=1−600/6000=0.90P = 1 - 600/6000 = 0.90.
Answer: The sandstone of the shoreface layer is almost certainly the same in both wells (0.90). The channel sandstone is the same only with a probability of about 0.35, so these are more likely different bodies.

What it means. Sandstones at the same depth that look alike on the logs are not necessarily one layer. Everything depends on the depositional environment, which determines the size of the bodies. In a “jigsaw”, we make the correlation of channel sandstones as one of several versions and test it with pressures, contacts and the response to injection. In the geological model such bodies are placed statistically. The width of the channel belt in the example is an estimate from analogues, and its uncertainty passes directly into the probability of connection. The data in this example are assumed for illustration.

Probability that a sand body continues to the neighbouring well0300600900120015000.000.250.500.751.00Distance between wells, mProbability of continuation0.350.90400 m924 m6 km

Plot: the horizontal axis is the distance between wells from 0 to 1500 metres, the vertical axis is the probability from 0 to 1. Three straight lines fall from one: for a width of 400 metres to zero already at 400 metres, for 924 metres to zero at 924 metres, for 6 kilometres almost horizontally, to 0.75 at 1500 metres. Points are marked at 600 metres: 0.35 on the middle line and 0.9 on the upper one.

CalculatedThe probability that a sand body continues to the neighbouring well, plotted against the well spacing. The curves are for body widths along the line of the wells of 400 m, 924 m (the channel belt of the example) and 6 km (the shoreface layer). The points show the example at a spacing of 600 m.
  • The depositional environment suggests the size and shape of the bodies. We learn it from core, the vertical sequences of rocks and the shapes of the curves (Facies, facies successions and sedimentary logs).
  • The net-to-gross ratio (the share of reservoir rock in the section): at a high ratio even channel sandstones merge with one another and form a connected body.
  • Pressures and contacts: a common pressure gradient and a common oil–water contact show that the bodies communicate. Different gradients and contacts show separate bodies (Pressure measurement with the wireline formation tester).
  • The behaviour of wells: the response of producers to injection, and interference tests (a pressure change in one well is watched in another).
  • Seismic data: the shapes of bodies and their boundaries, if the bodies are thicker than its resolution.

If data are lacking, we build several versions of the correlation into the geological model and examine how the hydrocarbons in place and the production forecast depend on the choice (Geological and technical uncertainty of a project). New wells and production history then weed out the unsuccessful versions.

Terms and symbols

Well correlation
Matching of sections: which layers and surfaces of one well correspond to which in another. General correlation works by formations and horizons, and detailed correlation works by layers and thin beds.
Marker bed (marker horizon)
A thin, persistent layer or surface that can be recognised in all wells and that was deposited at the same time. It is the framework of correlation.
Lithostratigraphic and chronostratigraphic correlation
Lithostratigraphic correlation joins like rocks (the boundaries are diachronous). Chronostratigraphic correlation joins surfaces of equal age (time lines do not cross).
Structural and flattened cross-sections
A structural cross-section places the wells at depths below sea level and shows the present-day attitude. A flattened cross-section shifts the wells so that a marker is horizontal and shows the relations at the time the marker was deposited.
Stratigraphic datum
The marker on which the panel is flattened. It must have been deposited horizontally and be present in all wells.
TVT and TST
The true vertical thickness hTVTh_\text{TVT}, and the true thickness normal to the layers, h=hTVTcos⁡δh = h_\text{TVT}\cos\delta.
Line of correlation
A straight line on a plot of the depths of the same events in two wells. The slope is the ratio of accumulation rates, and a horizontal step is a hiatus.
Pinch-out
Thinning of a layer to zero between wells. On the panel, the lines of the top and the base converge.
“Layer-cake” and “jigsaw”
A layer-cake is a section made of persistent layers. A jigsaw is a section made of bodies that are comparable to the well spacing and connected only partly.

Common mistakes and misconceptions

Wrong

“Sandstones at the same depth in two wells are one layer.”

In fact

Depth proves nothing: in a prograding succession and in river systems, different bodies lie at the same level. Markers above and below the layer, pressures, contacts and the response to injection prove a connection.

Wrong

“In a deviated well the layer is 7 m thinner — so there is a fault.”

In fact

First convert the interval to true vertical thickness or true thickness: a hole running against the dip “squeezes” the layer in true vertical depth. In the example of the lesson the “loss” of 7 m is pure geometry.

Wrong

“A sharp boundary on the resistivity log is a good marker.”

In fact

It may be an oil–water contact. A contact is horizontal in depths below sea level, cuts across the layers, and cannot be seen on the GR.

Wrong

“Flattening on a marker shows the ancient relief.”

In fact

Only if the marker was deposited horizontally and the differences in compaction are small. Shales that drape a reef, and progradational layers, were inclined from the start.

Wrong

“The base of a species range from cuttings can be used for correlation in the same way as the top.”

In fact

Fragments caved from the borehole wall shift the base downwards, by 21 m in the example of the lesson. From cuttings the tops are reliable. We take bases from core.

In depth: advanced detail

The composite standard: how graphic correlation builds a scale

In no single well does a species occur over the whole interval of its existence. In some places it was absent from the given environment, and in others the section is eroded. Graphic correlation lets us restore the full intervals step by step. We match the reference well with another well. By the line of correlation, the top of a species in the second well may turn out to be higher than in the reference well. Then we “extend” the interval in the reference well upwards to this level. If the base is lower, we extend the interval downwards. This gives a composite standard, a conventional scale in which every species has the fullest known boundaries. We then match all the wells of the area with it in turn, and each new well refines it.

It is important that when we extend intervals, we move the tops only up and the bases only down. The absence of a species from a sample may mean that it was simply not preserved there, while its presence in core cannot be false. For cuttings, with their caving, the rule works only for tops. At each step we draw the line of correlation through the most reliable events (ash beds, polarity reversals and log markers), and we check the other events against it.

Where the ratio of thicknesses in a deviated well comes from

Let a hole of length Δl\Delta l cross a layer. Its projection onto the normal to the bedding (the part of its length that points along that direction) is the true thickness: h=Δlcos⁡γh = \Delta l\cos\gamma. The vertical thickness is the distance between the top and the base measured vertically at any point. It is related to the true thickness through the dip angle: hTVT=h/cos⁡δh_\text{TVT} = h/\cos\delta. The projection of the hole onto the vertical is ΔH=Δlcos⁡θ\Delta H = \Delta l\cos\theta. Hence ΔH/hTVT=cos⁡θcos⁡δ/cos⁡γ\Delta H/h_\text{TVT} = \cos\theta\cos\delta/\cos\gamma.

For a hole in the vertical plane of the dip, the formula simplifies. Against the dip, γ=∣θ−δ∣\gamma = |\theta - \delta|. For θ=55°\theta = 55° and δ=12°\delta = 12° the ratio is cos⁡55°cos⁡12°/cos⁡43°≈0.77\cos 55°\cos 12°/\cos 43° \approx 0.77, so the interval is “squeezed” by almost a quarter. Down the dip, γ=θ+δ\gamma = \theta + \delta, and the ratio is cos⁡55°cos⁡12°/cos⁡67°≈1.44\cos 55°\cos 12°/\cos 67° \approx 1.44, so the interval is “stretched” to almost one and a half times its thickness. Along the strike, cos⁡γ=cos⁡θcos⁡δ\cos\gamma = \cos\theta\cos\delta and ΔH=hTVT\Delta H = h_\text{TVT}: a hole running along the strike sees the vertical thickness without distortion.

Automatic correlation of logs

Correlation programs compare the normalised curves of two wells numerically. The simplest way is a sliding correlation coefficient. The program shifts a stretch of the curve of one well along the curve of the other well and looks for the shift with the greatest similarity. This is good for finding a single marker, but it cannot cope with the fact that intervals in neighbouring wells are stretched differently. Dynamic time warping methods let the program stretch and compress one curve relative to the other while keeping the order of the layers. They find the best match for the whole section.

No algorithm knows geology. It will readily join the top of one channel to the top of another, and it will miss a fault if the pattern of the curves is similar. So we use automatic correlation as a draft and to check the manual one, while the geologist places the markers, hiatuses, faults and pinch-outs.

Connections to other topics

Check yourself

  1. 1

    A section of a delta succession has to be flattened on a marker. Which surface should be chosen as the datum?

    1. A.The top of a channel sandstone penetrated by all the wells
    2. B.An angular unconformity surface
    3. C.A GR peak in a persistent member of marine shales above the productive interval
    4. D.The top of a reef buildup
    Show the answer

    Correct answer: C. A GR peak in a persistent member of marine shales above the productive interval

    Flooding shales were deposited almost horizontally and at the same time over the whole area. This is a time line. The top of a channel sandstone is diachronous and uneven, a reef originally stood above the sea floor, and an unconformity surface has erosional relief and cuts across different layers. Flattening on any of these creates a false “ancient relief”.

  2. 2

    The layers dip at 15° towards azimuth 90°. A well with an inclination of 45° runs towards azimuth 90° (down the dip) and crosses the layer over 50 m of hole. What difference in true vertical depth will its log show, and what is the vertical thickness of the layer?

    Show the answer

    For a hole down the dip, γ=45°+15°=60°\gamma = 45° + 15° = 60°, h=50cos⁡60°=25.0h = 50\cos 60° = 25.0 m, and hTVT=25.0/cos⁡15°≈25.9h_\text{TVT} = 25.0/\cos 15° \approx 25.9 m. In true vertical depth ΔH=50cos⁡45°≈35.4\Delta H = 50\cos 45° \approx 35.4 m. The interval is “stretched” by 37% compared with a vertical well. Without conversion the layer would seem thickened or repeated.

  3. 3

    On a graphic correlation plot (depths in the reference well on the horizontal axis, and depths in the well being matched on the vertical axis) the line forms a horizontal step. What does this mean?

    1. A.A hiatus in the reference well
    2. B.The accumulation rates in the two wells are the same
    3. C.An error in identifying the fossils
    4. D.Part of the section of the reference well is absent in the matched well — a hiatus or a sharply reduced interval
    Show the answer

    Correct answer: D. Part of the section of the reference well is absent in the matched well — a hiatus or a sharply reduced interval

    On the step the depth in the matched well is constant while the depth in the reference well changes. So the interval of the reference well that corresponds to the step is not recorded in the matched well. Equal rates would give a slope of 1. A hiatus in the reference well would give a vertical step.

  4. 4

    We correlate on resistivity logs in four wells on the crest of an anticline. A sharp boundary lies at the same depth below sea level in all the wells, although the GR markers rise towards the crest by 30 m. What is this boundary, and how can we check it?

    Show the answer

    Almost certainly it is not a layer but an oil–water (or gas–oil) contact. It is horizontal in depths below sea level and cuts across layers that rise towards the crest. At this depth nothing changes on the GR. On the SP, at most, the deflection above the contact is slightly smaller. We check it with pressures (a change of gradient), testing and core. We must not join it as a lithological boundary.

  5. 5

    The distance between wells is 400 m, the width of the channel bodies is about 500 m, and their axis is perpendicular to the line of the wells. What is the probability that the channel sandstone of one well continues to the other?

    1. A.0.2
    2. B.0.5
    3. C.0.8
    4. D.1.25
    Show the answer

    Correct answer: A. 0.2

    P=1−400/500=0.2P = 1 - 400/500 = 0.2. Most probably the wells have different bodies. We can join them only as one of the versions, and we must test it with pressures and the response to injection.

  6. 6

    An injection well works in a layer that, by “sand-to-sand” correlation, is connected with a producer 1.5 km away. There has been no response in two years. Which alternative correlations does this explain, and what data would test them?

    Show the answer

    They may be different bodies: prograding layers separated by flooding shales (chronostratigraphic correlation), or channels that do not reach one another. A sealing fault is also possible. Checks: markers above and below the layer and biostratigraphy (are the sandstones of the same age?), formation-tester pressures (is it one hydrodynamic system?), contacts, seismic data and interference tests.

Updated: 15 October 2026