scholarly journals Deep seismic imaging of the eastern Nankai trough, Japan, from multifold ocean bottom seismometer data by combined travel time tomography and prestack depth migration

Author(s):  
J.-X. Dessa ◽  
S. Operto ◽  
S. Kodaira ◽  
A. Nakanishi ◽  
G. Pascal ◽  
...  
Solid Earth ◽  
2019 ◽  
Vol 10 (3) ◽  
pp. 765-784 ◽  
Author(s):  
Andrzej Górszczyk ◽  
Stéphane Operto ◽  
Laure Schenini ◽  
Yasuhiro Yamada

Abstract. Imaging via pre-stack depth migration (PSDM) of reflection towed-streamer multichannel seismic (MCS) data at the scale of the whole crust is inherently difficult. This is because the depth penetration of the seismic wavefield is controlled, firstly, by the acquisition design, such as streamer length and air-gun source configuration, and secondly by the complexity of the crustal structure. Indeed, the limited length of the streamer makes the estimation of velocities from deep targets challenging due to the velocity–depth ambiguity. This problem is even more pronounced when processing 2-D seismic data due to the lack of multi-azimuthal coverage. Therefore, in order to broaden our knowledge about the deep crust using seismic methods, we present the development of specific imaging workflows that integrate different seismic data. Here we propose the combination of velocity model building using (i) first-arrival tomography (FAT) and full-waveform inversion (FWI) of wide-angle, long-offset data collected by stationary ocean-bottom seismometers (OBSs) and (ii) PSDM of short-spread towed-streamer MCS data for reflectivity imaging, with the former velocity model as a background model. We present an application of such a workflow to seismic data collected by the Japan Agency for Marine-Earth Science and Technology (JAMSTEC) and the Institut Français de Recherche pour l'Exploitation de la Mer (IFREMER) in the eastern Nankai Trough (Tokai area) during the 2000–2001 Seize France Japan (SFJ) experiment. We show that the FWI model, although derived from OBS data, provides an acceptable background velocity field for the PSDM of the MCS data. From the initial PSDM, we refine the FWI background velocity model by minimizing the residual move-outs (RMOs) picked in the pre-stack-migrated volume through slope tomography (ST), from which we generate a better-focused migrated image. Such integration of different seismic datasets and leading-edge imaging techniques led to greatly improved imaging at different scales. That is, large to intermediate crustal units identified in the high-resolution FWI velocity model extensively complement the short-wavelength reflectivity inferred from the MCS data to better constrain the structural factors controlling the geodynamics of the Nankai Trough.


Geophysics ◽  
2004 ◽  
Vol 69 (4) ◽  
pp. 1053-1070 ◽  
Author(s):  
Einar Iversen

The isochron, the name given to a surface of equal two‐way time, has a profound position in seismic imaging. In this paper, I introduce a framework for construction of isochrons for a given velocity model. The basic idea is to let trajectories called isochron rays be associated with iso chrons in an way analogous to the association of conventional rays with wavefronts. In the context of prestack depth migration, an isochron ray based on conventional ray theory represents a simultaneous downward continuation from both source and receiver. The isochron ray is a generalization of the normal ray for poststack map migration. I have organized the process with systems of ordinary differential equations appearing on two levels. The upper level is model‐independent, and the lower level consists of conventional one‐way ray tracing. An advantage of the new method is that interpolation in a ray domain using isochron rays is able to treat triplications (multiarrivals) accurately, as opposed to interpolation in the depth domain based on one‐way traveltime tables. Another nice property is that the Beylkin determinant, an important correction factor in amplitude‐preserving seismic imaging, is closely related to the geometric spreading of isochron rays. For these reasons, the isochron ray has the potential to become a core part of future implementations of prestack depth migration. In addition, isochron rays can be applied in many contexts of forward and inverse seismic modeling, e.g., generation of Fresnel volumes, map migration of prestack traveltime events, and generation of a depth‐domain–based cost function for velocity model updating.


Geophysics ◽  
2009 ◽  
Vol 74 (6) ◽  
pp. Q41-Q47 ◽  
Author(s):  
Ranjan Dash ◽  
George Spence ◽  
Roy Hyndman ◽  
Sergio Grion ◽  
Yi Wang ◽  
...  

The subseafloor structure offshore western Canada was imaged using first-order water-layer multiples from ocean-bottom seismometer (OBS) data and the results were compared to conventional imaging using primary reflections. This multiple-migration (mirror-imaging) method uses the downgoing pressure wavefield just above the seafloor, which is devoid of any primary reflections but consists of receiver-side ghosts of these primary reflections. The mirror-imaging method employs a primaries-only Kirchhoff prestack depth migration algorithm to image the receiver ghosts. The additional travel path of the multiples through the water layer is accounted for by a simple manipulation of the velocity model and processing datum: the receivers lie not on the seabed but on a sea surface twice as high as the true water column. Migration results show that the multiple-migrated image provides a much broader illumination of the subsurface than is possible for conventional imaging using the primaries, especially for the very shallow reflections and sparse OBS spacing. The resulting image from mirror imaging has illumination comparable to the vertical incidence surface streamer (single-channel) reflection data.


Geophysics ◽  
2001 ◽  
Vol 66 (3) ◽  
pp. 721-732 ◽  
Author(s):  
Lanlan Yan ◽  
Larry R. Lines

Seismic imaging of complex structures from the western Canadian Foothills can be achieved by applying the closely coupled processes of velocity analysis and depth migration. For the purposes of defining these structures in the Shaw Basing area of western Alberta, we performed a series of tests on both synthetic and real data to find optimum imaging procedures for handling large topographic relief, near‐surface velocity variations, and the complex structural geology of steeply dipping formations. To better understand the seismic processing problems, we constructed a typical foothills geological model that included thrust faults and duplex structures, computed the model responses, and then compared the performance of different migration algorithms, including the explicit finite difference (f-x) and Kirchhoff integral methods. When the correct velocity was used in the migration tests, the f-x method was the most effective in migration from topography. In cases where the velocity model was not assumed known, we determined a macrovelocity model by performing migration/velocity analysis by using smiles and frowns in common image gathers and by using depth‐focusing analysis. In applying depth imaging to the seismic survey from the Shaw Basing area, we found that imaging problems were caused partly by near‐surface velocity problems, which were not anticipated in the modeling study. Several comparisons of different migration approaches for these data indicated that prestack depth migration from topography provided the best imaging results when near‐surface velocity information was incorporated. Through iterative and interpretive migration/velocity analysis, we built a macrovelocity model for the final prestack depth migration.


Geophysics ◽  
2009 ◽  
Vol 74 (6) ◽  
pp. WCA35-WCA45 ◽  
Author(s):  
Chaoshun Hu ◽  
Paul L. Stoffa

Subsurface images based on low-fold seismic reflection data or data with geometry acquisition limitations, such as obtained from ocean-bottom seismography (OBS), are often corrupted by migration swing artifacts. Incorporating prestack instantaneous slowness information into the imaging condition can significantly reduce these migration swing artifacts and improve image quality, especially for areas with poor illumination. We combine the horizontal surface slowness information of observed seismic data with Gaussian-beam depth migration to implement a new slowness-driven Gaussian-beam prestack depth migration whereby Fresnel weighting is combined naturally with beam summation. The prestack instantaneous slowness information is extracted from the original OBS or shot gathers using local slant stacks and is combined with a local semblance analysis. During migration, we propagate the seismic energy downward, knowing its instantaneous slowness information. At each image location, the beam summation is localized in a resolution-dependent Fresnel zone; the instantaneous slowness information controls the beam summation. Synthetic and real data examples confirm that slowness-driven Gaussian-beam migration can suppress most noise from inadequate stacking and give a clearer migration result.


Sign in / Sign up

Export Citation Format

Share Document