SEP131 -- TABLE OF CONTENTS |

A complex velocity model produces shadow zones in an image due to focusing and defocusing of the seismic waves, and limited recording geometry. These shadow zones contain weak signal masked by artifacts. To recover the real signal, and reduce artifacts is necessary to go beyond migration. One option is to use a wave-equation target-oriented inversion scheme that explicitly computes the least squares inversion Hessian. The wave-equation target-oriented inversion has a big null space. It seeks to form an image where there is lack or very little data information. In this situation is where

Plane-wave migration in tilted coordinates [pdf 1.6M][src 226Mb]

Most existing one-way wave-equation migration algorithms have difficulty in imaging steep dips in a medium with strong lateral velocity variation. We propose a new one-way wave-equation-based migration, called ``plane-wave migration in tilted coordinates.'' The surface data are converted to plane-wave source data by slant-stacking processing, and each resulting plane-wave source dataset is migrated independently in a tilted coordinate system with an extrapolation direction determined by the source plane-wave direction at the surface. For most waves illuminating steeply dipping reflectors, the extrapolation direction is closer to their propagation direction in the tilted coordinates. Therefore, plane-wave migration in tilted coordinates can correctly image steeply dipping reflectors, even by applying one-way extrapolators. In a well-chosen tilted coordinate system, waves that overturn in conventional vertical Cartesian coordinates do not overturn in the new coordinate system. Using plane-wave migration in tilted coordinates, we can image overturned energy with much lower cost compared to reverse-time migration.

Downward continuation migration cannot provide reliable angle domain common image gathers (CIGs) for steeply dipping reflectors, because it cannot handle most waves that illuminate steep reflectors. Also there is a severe stretch in conventional horizontal subsurface offset at steep reflectors. Both reverse-time migration and plane-wave migration in tilted coordinates solve these two problems and provide robust angle domain CIGs for steeply dipping reflectors. A test on the BP velocity benchmark dataset shows that both migration methods generate robust angle domain CIGs that are comparable. When the migration velocity is not correct, the angle domain CIGs from both migration methods show useful moveout information for velocity estimation.

Phase unwrapping of angle-domain common image gathers [pdf 336K] [src 1.8Mb]

In this paper we adapt a phase unwrapping algorithm to estimate the depth shift in Angle-Domain Common Image Gathers (ADCIGs). We show how to set up a linear system of equations tailored to the seismic case and how to solve it by minimizing an measure via iterations of weighted least-squares problems. For this procedure a meaningful choice of initial weights is crucial. We propose to unwrap jointly several angle gathers and show that this can overcome sampling deficiencies in the angle domain, such as those that come from processing a limited number of subsurface offsets for angle-gather generation.

Angle-domain parameters computed via weighted slant-stack [pdf 392K] [src 4.9Mb]

Angle-domain common image gathers (ADCIGs), created from downward-continuation or reverse time migration, can provide useful lithological and velocity information (Prucha et al., 1999). In geologically complex areas, poor illumination causes undesirable kinematic effects and amplitude variations along the angle axis (Prucha et al., 2000; Valenciano, 2006). The subsurface-offset-to-angle transformation consists of a radial trace transform in the Fourier space with some regularization in the angle direction (Sava and Fomel, 2000) or slant-stack in the physical space plus an additional transformation from offset ray-parameter to reflection angle (Prucha et al., 1999). The regularization, to some extent, can diminish the amplitude variation caused by poor illumination. The more accurate solution to the illumination problem, however, is achieved by computing a regularized least-squares inverse image (Clapp, 2005) rather than the simply the adjoint (migration). The inverse image problem can be solved either by computing the Hessian implicitly (Clapp, 2005) or explicitly Valenciano (2006), preferentially, in the reflection-angle domain or, without any physically meaningful regularization direction, in the subsurface-offset domain. ...

We study the effectiveness of velocity-only time-domain waveform inversion for inverting synthetic data modeled with both velocity and density contrasts. We present a detailed review of the Born approximation for the constant-density acoustic wave equation and its application to the inversion of velocity models for seismic reflection data. We create synthetic models with both constant and variable density and compare the effectiveness of velocity-only waveform inversion in each case. Results from this simple test suggest that density contrasts can hamper the reconstruction of velocity perturbations.

Exact seismic velocities for TI media and extended Thomsen formulas for stronger anisotropies [pdf 448K][src 1.2Mb]

I explore a different type of approximation to the exact anisotropic wave velocities as a function of incidence angle in transversely isotropic (TI) media. This formulation extends the Thomsen weak anisotropy approach to stronger deviations from isotropy without significantly affecting the simplicity of the equations. One easily recognized improvement is that the extreme value of the quasi-SV-wave speed is located near the correct incidence angle , rather than always being at the position 45, which universally holds for Thomsen's approximation - although is actually never correct for

Prediction Error Filters (PEFs) capture the covariance of a dataset. In this paper I use PEFs to quantify and highlight difference between two volumes. A series of PEFs are estimated on one volume and then applied to a second. The resulting hypercube is an indicator of where, and how much, two volumes differ.

Accelerating seismic computations using customized number representations on FPGAs [pdf 260K] [src 9.7Mb]

Field Programmable Gate Arrays (FPGA) offer significant potential speedups over conventional microprocessors for some applications. For downward continued migration, complex math and Fast Fourier Transforms (FFT) are the dominant cost. Convolution is the dominant cost for reverse time migration. We implement these core algorithms on a FPGA and show speedups ranging from 5 to 15, including the transfer time to and from the processors. We consider methods to further speed up these migration algorithms.

SEP131 -- TABLE OF CONTENTS |

2007-09-18