This paper describes the implementation of a 3-D pre-stack depth migration and its application to deep seismic reflection data. Based on the Kirchhoff migration approach, the migrated image is generated by weighted integration along diffraction surfaces through the shot record sections. A finite-difference solution of the eikonal equation is used to compute the diffraction surfaces for a given 3-D velocity model. A test on synthetic data and an application to a 3-D deep seismic reflection data set (ISO89-3D), recorded in the vicinity of the German continental deep drill hole (KTB), are detailed. Because of strong lateral velocity contrasts in some parts of this area, the migration was performed in two steps. First, parts of the shot records were migrated using a constant velocity in order to obtain a general overview of the subsurface structure. Afterwards, the whole data set was migrated using a 3-D velocity model derived mainly from wide-angle measurements. Slices through the migrated volume covering an area of about 21 km × 21 km in the horizontal and 15 km in the vertical direction were compared with results from a 2-D survey (KTB8502) carried out within the same region. Over nearly the whole survey area, the SE1 reflector shows up as the continuation of the Franconian Lineament with varying strength of reflection. Its signature, which deviates partly from that of a simple plane reflector, can be followed through the whole upper crust down to about 9 km depth. There, the Erbendorf body appears as a complicated, strongly reflective bowl-shaped pattern with its uppermost top reflection shifted about 5 km to the south-southeast of the KTB drill hole. The results make it clear that the lateral extent and the true depth of these main geological structures (SE1, Erbendorf body) can be interpreted correctly only from the 3-D migration result. Finally, the parallel implementation procedure is described, both for constant- and variable-velocity models.
[1]
W. Rabbel,et al.
Origin and nature of crystal reflections: Results from integrated seismic measurements at the KTB superdeep drilling site
,
1997
.
[2]
Jean Brac,et al.
Can we image complex structures with first‐arrival traveltime?
,
1993
.
[3]
M. Kvasnička,et al.
3-D network ray tracing
,
1994
.
[4]
William W. Symes,et al.
Upwind finite-difference calculation of traveltimes
,
1991
.
[5]
Migration of the 3-D deep-seismic reflection survey at the KTB location, Oberpfalz, Germany
,
1997
.
[6]
J. Lauterjung,et al.
The German Continental Deep Drilling Program KTB: Overview and major results
,
1997
.
[7]
P. Podvin,et al.
Finite difference computation of traveltimes in very contrasted velocity models: a massively parallel approach and its associated tools
,
1991
.
[8]
J. Vidale.
Finite‐difference calculation of traveltimes in three dimensions
,
1990
.
[9]
T. Moser.
Shortest path calculation of seismic rays
,
1991
.
[10]
H. Gebrande,et al.
Pre-stack migration and true-amplitude processing of DEKORP near-normal incidence and wide-angle reflection measurements
,
1996
.
[11]
David Lumley,et al.
Imaging complex geologic structure with single‐arrival Kirchhoff prestack depth migration
,
1997
.