[First Hit]

Datapages, Inc.Print this page

AAPG GEO 2010 Middle East
Geoscience Conference & Exhibition
Innovative Geoscience Solutions – Meeting Hydrocarbon Demand in Changing Times
March 7-10, 2010 – Manama, Bahrain

Non-Hyperbolic Reflection Tomography for Better Imaging and Interpretation

Zhaobo Meng2; Chengbin Peng1

(1) Nexus Geosciences Inc., Sugar Land, TX.

(2) ConocoPhillips, Houston, TX.


In this abstract we will demonstrate that the application of non-hyperbolic reflection tomography can produce a 3D subsurface Previous HitvelocityNext Hit model that is smooth for prestack Previous HitdepthNext Hit migration and at the same time geological for accurate subsurface pore pressure prediction and lithology interpretation. Our non-hyperbolic reflection tomography workflow is very different from the conventional reflection tomography workflow: no automatic volume picking of residual moveout is used. Instead a set of prestack events are interpreted on prestack seismic image volumes. The prestack events carry detailed moveout information that is more accurate than a single parameter fit to the common image gathers. In addition the prestack events carry structure Previous HitdipNext Hit information that is needed for accurate 3D ray tracing offset by offset. Examples will be given for both compaction driven and lithology controlled geological environments.


Most production tomography workflow relies on automatic volume picking of residual moveout coefficients, a single parameter description of the moveout variation across offset. Modern seismic acquisition has large source receiver offset (8 kilometer or more). In many minibasin plays the maximum offset is equal to or larger than the basin width. The underlying assumption in most hyperbolic tomography workflow starts to break down: the moveout variation on common image point gathers can no longer be described by a single parameter, the curvature. Hence the need of non-hyperbolic residual moveout descriptions !

The quantity “non-hyperbolic residual moveout” is “interpretive” because measuring non-hyperbolic moveouts requires prestack seismic interpretation tools and procedures, which is not available in most desktop applications.

When seismic data is migrated with the correct Previous HitvelocityNext Hit model, typically the resulting common offset image gathers are flat, as opposed to when the migration Previous HitvelocityNext Hit is not accurate which would result in offset gathers that are not flat. The “non-flatness” depends on the migration model, and is sensitive to other factors such as the Previous HitdipNext Hit field. When the model is simple, the offset gathers analytically follow a hyperbolic curve (Meng & Bleistein, 2001). Otherwise the residual moveout curve would be “non-hyperbolic” — i.e. anything that is not hyperbolic. Differentiating hyperbolic from non-hyperbolic moveouts (Meng et al., 2004) requires a significant change in the tomography workflow, as hyperbolic moveouts can be characterized by one parameter - the curvature - and most reflection tomographic approaches in the literature are based on the single-parameter hyperbolic moveout.

We describe our non-hyperbolic tomography workflow Previous HitusingNext Hit one synthetic example and one real data example. The first example is a synthetic example that simulates a compaction driven environment such as the GOM. The second example is a real data example in a lithology driven environment.
It is worthwhile to point out that non-hyperbolic tomography will be more applicable to VTI and TTI migration and Previous HitvelocityNext Hit model updating, as the migration operators are non-hyperbolic in nature. Previous HitUsingNext Hit hyperbolic moveout assumption for VTI and TTI tomography is logically inconsistent.

Synthetic Example: GOM Compaction Driven

The synthetic data was provided as a tomographic benchmark test. Only prestack migrated image gathers and the migration Previous HitvelocityNext Hit model were made available during this blind test. The true Previous HitvelocityNext Hit model used to generate the unmigrated synthetic data was made available only after the completion of the benchmark test.

Figure 1 shows both the common image gathers and our hyperbolic picks (shallow) and non-hyperbolic picks (deep) made on the gathers. Substantial residual moveouts in the whole data volume are visible.

A total of 20+ prestack events were picked on the full-volume common image gathers. In Figure 1 each color corresponds to an individual prestack event.

After prestack event tracking, the prestack events were demigrated and tomographic inversion was conducted. Fig. 2 (right) shows the final Previous HitvelocityNext Hit model perturbation from tomographic inversion. Figure 2 (left) shows the true Previous HitvelocityNext Hit model perturbation that was made available after completion of the benchmark test. As one can see, our tomography has resolved both the low Previous HitvelocityNext Hit anomaly and the two high Previous HitvelocityNext Hit anomalies very well: both in the magnitude of Previous HitvelocityNext Hit correction and in the geometry of the anomalies.

Real Data Example: Lithology Driven Environment

The real data example covers a 1600 km2 area in a lithology-driven environment. The goal of this study was to improve a Previous HitvelocityNext Hit model starting from existing Kirchhoff image gathers. We have employed 3D prestack demigration and remigration (Peng and Sheng, 2009) in the iterative tomography workflow. A third party tomography result was available for comparison. After analyzing the third party’s Previous HitvelocityNext Hit model, some issues were noted, including: (1) the tomographic result is not usable for lithological intrepretation (Fig. 5, left) because it does not look geological; (2) the starting image gathers show large residual errors in the shallow carbonates, while the gathers appeared reasonably flat in the deep target area under a Mass Transport Complex (MTC), a Previous HitvelocityNext Hit compensating effect over a Previous HitdepthNext Hit range. These errors in the carbonates introduce significant uncertainty to the interpretation in the deep siliclastical target area (Fig. 4, left).

With access only to the original Previous HitvelocityNext Hit model and PreSDM gathers, we: (1) applied 3D tomography only to the base of carbonates (Tertiary) to correct the large residual errors; (2) demigrated the image gathers Previous HitusingNext Hit the original Previous HitvelocityNext Hit model; (3) re-migrated Previous HitusingNext Hit the new supra-Tertiary Previous HitvelocityNext Hit model; (4) performed another 3D tomography to the deep target area.

Figure 3 (left) shows the original multiple-stage tomographic Previous HitvelocityNext Hit model overlaid on the original stacked image. Figure 3 (right) shows our tomographic Previous HitvelocityNext Hit model overlaid on the same stacked image. Our tomography reveals the high resolution details that are consistent with the lithology. The flatness of the new image gathers (Fig. 4, right) validates our tomography result. This is in contrast to the original tomographically updated model (Fig. 4, left). The lower Previous HitvelocityNext Hit (in blue color) may be caused by over-pressurized fluids trapped under the carbonates (Fig. 3 right).
Figure 5 shows an interesting area that involves a normal fault in extension setting. The original model was built by a third party Previous HitusingNext Hit a tomographic workflow (Fig. 5, left). Fig. 5 (right) shows our tomographic results which reveal a higher Previous HitvelocityNext Hit on the up-throw side (older stratigraphy), and a lower Previous HitvelocityNext Hit on the down-throw side (younger sedimentary rocks) across the fault. The fault clearly defines a boundary in our Previous HitvelocityNext Hit model. The new model shows higher Previous HitvelocityNext Hit at the synclines and lower Previous HitvelocityNext Hit at the crests, the subtle Previous HitvelocityNext Hit difference caused possibly by a fluid pressure effect. Supported by both flat image gathers and geological consistency, this kind of Previous HitvelocityNext Hit/impedance inversion may provide invaluable information about fluid migration and caprock sealing effects.


We have developed a 3D tomographic Previous HitvelocityNext Hit updating workflow that uses the non-hyperbolic moveout residues on common image gathers. The enabler of this workflow is a prestack seismic interpretation system that permits prestack event tracking on large prestack seismic volumes. A variety of interpretive tools have also been built around the 3D tomography workflow to enable quick and efficient incorporation of geological and geophysical information/knowledge into the process, including QC and editing.

The tomographic toolkits and the workflow have been tested in GOM salt model building and non-GOM lithology-driven areas. In both cases, the 3D tomographic results proved to be geologically consistent. Our final Previous HitvelocityNext Hit models also provide higher resolution details in the target areas that can be used for other applications such as pore pressure prediction and fluid migration evaluations.


The authors sincerely thank BHP Billiton Petroleum USA Inc and its partners for permission to publish these results and the seismic data examples. The synthetic data is provided by ConocoPhillips.

680352_A.jpgFig. 1. Non-hyperbolic moveout picks are produced by interpreting horizons on the near offset volume and propagating the seed horizons automatically in offset dimension. They are used as inputs in tomographic updating.

680352_B.jpgFig. 2. Synthetic model perturbation (Left) and non-hyperbolic tomographic inversion (Right). Synthetic prestack data is migrated Previous HitusingNext Hit a Previous HitvelocityNext Hit model that is a function of V(z) only. Common imaging gathers and the V(z) model were available for the tomography benchmark test. The true Previous HitvelocityNext Hit model was disclosed to all participants after submission of the benchmark result. Our non-hyperbolic tomography yields Previous HitvelocityNext Hit perturbations that are very close to the truth.

680352_C.jpgFig. 3. In this non-GOM example, hard carbonate (limestone) covers a lower Previous HitvelocityNext Hit mass transport complex (MTC). The left picture shows the original multiple stage tomographic Previous HitvelocityNext Hit model overlaid on the original stacked image; the right picture shows the higher resolution of the new model overlaid on the same stacked image.

680352_D.jpgFig. 4. The left picture shows the gathers before the interpretive tomography (Previous HitusingNext Hit the original tomographically updated model) and right shows the gathers after the interpretive tomography with some muting applied

680352_E.jpgFig. 5. The left picture shows the original Previous HitvelocityNext Hit overlaid on the stacked image, while the right picture shows the new Previous HitvelocityNext Hit overlaid on the same stacked image. Our tomography generates a model that shows lower Previous HitvelocityNext Hit at the down throw and higher Previous HitvelocityNext Hit at the up-throw side of a normal fault due to expansion. On the left, the new model reveals higher Previous HitvelocityNext Hit at the synclines and lower Previous HitvelocityTop at the crests, the velocities on the right are likely caused by fluid pressure effect.