Nine modules of losses, ambiguities and error bars. This one adds them up. It counts how
many rocks are consistent with one measurement, how many the second measurement removes, which
question neither of them can answer, and it ends with a form of words that the numbers will
support.
This is the closing module and it assumes all eight before it.
The rocks that still fit — before and after the gradient
The measurement — and how well it was made
What was recorded
What is under the anomaly
This setting is not visible to the analysis below, which is given the amplitudes and
nothing else.
Work through the five steps in order, then try the exercises.
Method and references are in the last two tabs.
1Where the information went
Every step, and what each one cost
A rock has many independent properties. A seismic amplitude has one number, or two
where offset is available. The difference has to go somewhere, and the earlier modules located it.
The card on the left collects those losses, with the size of each one measured rather than
described.
The losses are of different kinds. Some are
collapses: several numbers becoming one, as when velocity and density become an impedance.
Some are saturations: a property that stops responding once it has changed a little, as gas
saturation does. Some are blurs: the wavelet in time, the noise in amplitude. And one is a
confusion: half the reflection belongs to a rock that was never the subject of the
question.
Only the blurs get better with better data. The collapses and the saturation are
properties of the physics and no amount of acquisition or processing touches them.
1 Every loss, and where it was measured
numbers from the earlier modules
2 Which of them better data can fix
and which are permanent
Losses that better data can reduce—
Losses that it cannot—
Properties going in—
Numbers coming out—
Measured with—
The one question none of them answers—
None of this is an argument against using amplitudes. It is an argument for establishing which
questions the data can answer before they are asked. The next three steps put that on a countable
footing.
2What one number leaves
Module 07's count, now with a real error bar on it
Module 07 counted the rocks consistent with a measured amplitude, using a tolerance
picked to represent the noise. Now that tolerance can be computed properly: module 08 gives the
uncertainty in a fitted intercept from the angle range, the fold and the noise, and the set of
admissible rocks is everything within two of those standard deviations.
The panel shows the survivors as a population — what porosities, what fluids, what
saturations. The count is large, and the important thing about it is not the number itself but which
of its dimensions are wide and which are narrow.
Changing the noise widens or narrows the whole set, and changing the mute does it
again. Neither of them changes the geology by so much as a percent of porosity.
1 The rocks that fit the intercept
porosity, and what is in the pores
2 How wide each unknown still is
longer is worse
Models searched—
Consistent with the intercept—
Tolerance on the intercept—
Porosity still admitted—
Saturation still admitted—
Fraction that contain any hydrocarbon—
The fraction containing hydrocarbon is a count of grid cells under equal weighting, not a
probability. It says a bright amplitude is much easier to build with hydrocarbons than without, on
this grid, which is why bright spots get drilled. It does not say the odds.
3What the second number removes
The whole case for AVO, as a number
Now require the gradient to match as well, within its own two standard deviations. The
surviving set shrinks. How much it shrinks is the practical case for doing AVO, stated as a
figure rather than as a conviction.
At moderate noise the gradient removes 55% of what the intercept left. That is a real and
useful return. At high noise it removes 25%, because the gradient's error bar
has grown so wide that it barely constrains anything, which is module 08's result in a form that
can be costed.
Then look at which dimensions shrank. Vp/Vs narrows sharply, which is exactly
what module 05 predicted: offset measures shear. Porosity narrows a little. And saturation does not
move at all. The second measurement constrains lithology and fluid type, and not how much is
there.
1 Before and after the gradient
each unknown, both times
2 What that removed
the survivors, on the crossplot
Consistent with the intercept—
…and with the gradient—
The second measurement removed—
Vp/Vs range, before and after—
Saturation range, before and after—
What adding offset removed here—
At the noisy setting the second column collapses toward the first. AVO is not a method that
either works or does not; its value is proportional to how well the gradient can be measured, and
module 08 set out what that depends on.
4Four ways to be confidently wrong
Each one is a rule of thumb that is usually true
Every trap below comes from a rule that works most of the time, which is why each has
survived in use: a rule that failed often would have been discarded. Each button loads the situation
and the analysis that goes with it.
1 The trap
what the rule says, and what is true
2 The two rocks, as gathers
what the data would have shown
The rule—
Where it fails—
What the rule expects—
What is actually there—
Measurable difference—
Where this was established—
The common thread is that each rule is a statement about a reservoir, while an
amplitude is a statement about a contrast — between two rocks, at a thickness, through a
wavelet, with noise on top. Every trap is that mismatch showing through.
5What the numbers support
Generated from the surviving set, not from a template
The statement below is assembled from the analysis on the previous steps. Each clause
is backed by a count or a range from the surviving models: nothing in it is asserted that the
numbers do not support, and nothing that they do support is left out.
Changing the noise, the mute or the fold changes the wording with them. What a
statement can claim depends on how well the measurement was made, so it is re-derived for each set
of conditions rather than carried over.
The shape of it is consistent. There is usually something firm to say about impedance
and about the presence of a hydrocarbon, usually nothing to say about saturation, and almost always
a caveat about the seal, because half of every amplitude belongs to a rock nobody characterized.
1 The statement
every clause traceable to a count
2 What would narrow it most
the value of the next measurement
Is a hydrocarbon established—
Is the amount established—
Is the porosity established—
Is the thickness established—
Most valuable next measurement—
The one-line version—
One thing this page leaves out: an anomaly that holds together across a mapped area is far
better evidence than any single location, because the noise is independent from trace to trace and
the geology is not. Everything here is the worst case of a single gather. Real interpretation is
better than that, though not in the dimensions the physics has closed off.
WWhy it matters
Using a method whose limits are known
The point of the whole set is not that AVO is unreliable. It is that AVO answers
some questions well, some badly and some not at all, and that which is which can be worked out in
advance rather than discovered after a well.
Ask the question the data can answer. "Is there a hydrocarbon here" is often answerable.
"How much" usually is not. Framing a prospect around the second one guarantees a disappointing
result no matter how good the processing is.
Compute the value of the gradient before paying for it. Step 3 puts a number on what the
second measurement removes. At high noise it removes a fifth of the solutions, which may not
justify the work.
Spend on what narrows the set. Step 5 ranks the candidates. Very often the cheapest
useful measurement is a constraint on the seal rather than anything about the reservoir.
Report a set, not a value. Every number in a quantitative interpretation is a
representative of a set of possibilities. Quoting the representative without the set is the single
most common way these methods mislead.
The practice that follows from all of this is to look for the other models once one
has been found that explains the data. The first fit reads as an answer, and it is one member of a
set whose size has not yet been checked. Everything in these ten modules is an elaboration of
that.
This site is itself a model, and inherits the limitation it describes. The rock
physics is one choice among several, the noise is simpler than real noise, nothing here is spatial,
and a real interpreter has structural conformance, well control and a hundred other things this page
does not. The numbers describe the parameter space that was searched, which is not the earth.
The result the whole set converges on: an amplitude constrains a contrast rather than a
reservoir, and the size of the set it leaves is computable before anything is drilled.
Exercises
Each answer is a number read off this page at the settings
described.
Price the gradient. In step 3, read how many models survive the intercept alone and how
many survive both numbers, at moderate noise and then at noisy.
AnswerAt moderate noise 441 models match the intercept and 198 match
both, so the gradient removes 55% of them. At high noise 1,156 match the intercept and 870
match both — the gradient removes only 25%. Same rocks,
same acquisition geometry, same physics. What changed is how well the gradient could be
measured, and module 08 explains precisely why. If a project cannot afford the processing to
get the noise down, the AVO may be worth much less than the proposal assumed.
Find the dimension that never narrows. Still in step 3, compare the Vp/Vs range and the
saturation range before and after the gradient is applied.
AnswerAt moderate noise Vp/Vs goes from 1.46–2.08 to 1.46–1.60, a
real narrowing, and it is why the hydrocarbon-bearing fraction goes from 98% to 100%: offset
measures shear, and the wet rocks are the ones with high Vp/Vs. Saturation goes from 0–100% to
10–100%, and the only thing that changed is that the wet rocks were removed. Among the models
that do contain a hydrocarbon, every saturation from 10% to 100% still fits exactly as well. That is the fizz problem from modules
02 and 03 arriving at its final form: no measurement on this page, at any noise level, can tell
a commercial column from a residual one.
Watch the mute change the conclusion. In step 5, set the widest angle to 15°,
read the statement, then take it to 45° and read it again.
AnswerAt a 15° mute the gradient removes 4% of the models the
intercept left — its error bar is so wide that it constrains almost nothing, and the statement
drops back to what an intercept alone can support. At 30° it removes 55%. At 45° it removes
75%. Nothing about the earth changed between the two readings. It is the most direct
demonstration in the site that a processing parameter decides what can be concluded.
Take each trap in turn. In step 4, work through all four and note which module first
established each one.
AnswerBright means gas fails because
a gas sand at 20% porosity gives an intercept of −0.004 and is invisible — module 03, exercise
2. Dim means wet fails for the same reason from the other side. A lot of gas fails because 10%
saturation produces two thirds of the full amplitude effect — module 02, exercise 3, and
module 07's contour. Class III means gas fails because the same gas sand is Class I, II or III
depending only on the shale above it — module 06, exercise 5. Four rules, four
counterexamples, all built from the same physics rather than from anecdotes.
Draft the statement first, then compare. Set the noise to moderate and the mute to
30°, write down what the data supports, then open step 5.
AnswerThe defensible version is
close to: the impedance is constrained to a band; the response is very hard to produce without
a hydrocarbon on this parameter grid; the fluid is more likely gas than oil on the gradient;
the saturation is entirely unconstrained between about 5% and 100%; and the porosity depends on
an assumption about the shale that has not been tested. Most people's first draft is firmer
than that on saturation and vaguer than that on the shale, which is exactly the pair of errors
this site was built to correct.
Key points
The losses are of different kinds and only some are fixable. Collapses like velocity and
density into impedance, and the saturation of the gas effect, are physics. The blurs — the
wavelet, the noise — are the only ones better data improves.
One number leaves a large set. With a properly computed tolerance, 441 rocks out of the
13,020 searched match the measured intercept at moderate noise, spanning porosities from 19% to
35% and every saturation from zero to full.
The second number removes about half of them, and only when it can be measured well. The gradient cuts 441 to 198 at moderate noise and 1,156 to
870 when noisy — 55% against 25%. Widen the mute instead and it cuts 75%; narrow it to 15° and it
cuts 4%.
The value of AVO is proportional to how well the gradient is determined.
What offset constrains is Vp/Vs, and therefore fluid type. The admitted Vp/Vs range narrows from 1.46–2.08 to 1.46–1.60,
and the hydrocarbon-bearing fraction goes from 98% to 100%.
Saturation is untouched by all of it. Among the models that contain a hydrocarbon, before and after
the gradient and at every noise level, the admitted saturation runs from 10% to 100%. Nothing in normal-incidence or
offset-dependent amplitude separates a commercial column from a residual one.
Every rule of thumb has a computable counterexample. Bright is not gas, dim is not wet,
a strong anomaly is not a lot of gas, and Class III is not a rock type — the same sand moves
through three classes on the shale slider alone.
MMethod, simplifications and references
What is computed, and where it departs from production software
Everything here is built from the earlier modules. The rock physics is module
01 and 02's, the reflection coefficients are exact Zoeppritz as in module 05, the intercept and gradient
are fitted as in module 06, and their uncertainties come from module 08's closed form. Each Method tab
lists its own assumptions and all of them apply here.
The tolerances are computed, not chosen. A model is kept if its intercept is
within two standard deviations of the measured one, and — in step 3 — its gradient likewise. Those
standard deviations follow from the angle range, the number of traces and the noise level on the
sliders, so moving any of those genuinely changes what is admissible rather than just relabelling
it.
The search grid is porosity 5–35% in 1% steps, clay 0–30% in 10% steps, fluid
brine, oil or gas, hydrocarbon saturation 10–100% in 10% steps, and shale velocity 2300–3100 m/s in
200 m/s steps, with a hydrocarbon at zero saturation counted once as brine rather than three times.
The fit for each model uses eight angles rather than sixteen, which changes the recovered gradient by
about 0.001 and makes the sweep fast enough to run while a slider moves.
Counts are not probabilities. This is the same caution as module 07's and it
matters more here, because the closing statement is generated from the counts. The survivor count is
the number of cells in a hand-chosen grid, weighted equally, that reproduce the measurement. Turning
it into a probability needs a prior that reflects the basin, the depositional setting and the trap,
and different priors give materially different answers from identical seismic. Read every count as
"this many distinct explanations exist", which is a statement about the measurement, and never as
"this is how likely each is", which is a statement about geology that seismic alone cannot make.
Where this departs from production practice.
Nothing is spatial. This is the largest omission on the page. Real interpretation gains
enormously from an anomaly conforming to structure and repeating over an area, and none of that is
available to a single-gather analysis.
No well control, no seismic inversion, no other attributes, all of which narrow the set
further and are the reason quantitative interpretation is worth doing at all.
The truth is known to the page and hidden from the analysis. That is a luxury: on real
data nobody can check the answer, which makes the discipline of counting alternatives more
important rather than less.
Thickness is fixed at 45 m and tuning is therefore not part of the search here, though
module 07 showed it widening the admissible set considerably.
The generated statement is a teaching device. It is assembled from the counts by a fixed
set of rules. It is not a report, and a real one would carry the geological reasoning that this
page has no access to.
References
Tarantola, A., 2005, Inverse Problem Theory and Methods for Model Parameter
Estimation: SIAM.
Cambois, G., 2000, Can P-wave AVO be quantitative?: The Leading Edge, 19, 1246–1251.
Avseth, P., Mukerji, T., and Mavko, G., 2005, Quantitative Seismic Interpretation:
Cambridge University Press.
Bosch, M., Mukerji, T., and González, E. F., 2010, Seismic inversion for reservoir properties
combining statistical rock physics and geostatistics: a review: Geophysics, 75,
75A165–75A176.
Simm, R., and Bacon, M., 2014, Seismic Amplitude: An Interpreter's Handbook:
Cambridge University Press.
Rutherford, S. R., and Williams, R. H., 1989, Amplitude-versus-offset variations in gas sands:
Geophysics, 54, 680–688.
That is the set
Ten modules, one idea
An amplitude constrains a contrast rather than a reservoir, and the size of the set it leaves is
computable before anything is drilled. Module 01, where the set starts, builds a rock and follows it
to a trace.