Spatial aliasing in multichannel seismic data can be overcome by solving an inversion in which the model is the section that would be recorded in a well sampled zero-offset experiment, and the data are seismic data after normal moveout (NMO). The formulation of the (linear) relation between the data and the model is based on the wave equation and on Fourier analysis of aliasing. A processing se...

[1] The volcanic eruption on Miyakejima, Japan, in 2000 was marked by the largest earthquake swarm ever recorded in Japan, a seismicity migration accompanying a dike intrusion as the dike propagated from Miyakejima to the northwest, and formation of a caldera on Mount Oyama on Miyakejima. In this study, we propose a seismic source model which can be used to model both seismic and geodetic displ...

Marchenko redatuming allows one to use surface seismic reflection data to generate the seismic response at any point in the subsurface due to sources at the surface. Without requiring much information about the earth’s properties, the seismic response generated by Marchenko redatuming contains accurate estimates of not only the primaries, but also internal multiples. A target-oriented imaging m...

This paper presents least-squares wave equation AVP (Amplitude versus Ray Parameter) migration with non-quadratic regularization. We pose migration as an inverse problem and propose a cost function that makes use of a priori information about the AVP common image gather. In particular, we introduce two regularization goals: smoothness along the ray parameter direction and sparseness in the dept...

We introduce a novel finite-difference (FD) approach for seismic wave extrapolation in time. We derive the coefficients of the finite-difference operator from a lowrank approximation of the space-wavenumber, wave-propagator matrix. Applying the technique of lowrank finite-differences, we also improve the finite difference scheme of the two-way Fourier finite differences (FFD). We call the new o...

We introduce a novel technique for seismic wave extrapolation in time. The technique involves cascading a Fourier Transform operator and a finite difference operator to form a chain operator: Fourier Finite Differences (FFD). We derive the FFD operator from a pseudo-analytical solution of the acoustic wave equation. 2-D synthetic examples demonstrate that the FFD operator can have high accuracy...

Joint inversion of different geophysical datasets is an effective way to eliminate non-uniqueness in geophysical inversion problems. In this paper, we focus on a case study of joint inversion of seismic traveltime and gravity observations. The results are encouraging and we can have confidence that, comparing traveltime tomography alone, joint inversion of seismic traveltime and gravity data im...

Reliable seismic depth migrations require an accurate input velocity model. Inaccurate velocity estimates will distort point diffractors into smiles or frowns on a depth section. For both poststack and prestack migrated sections, high velocities cause deep smiles while low velocities cause shallow frowns on migrated gathers. However, for prestack images in the offset domain, high velocities cau...

Reverse time migration (RTM) is a commonly employed imaging technique in seismic applications (e.g., to image reservoirs of oil). Its standard implementation cannot account for multiple scattering/reverberation. For this reason it has not yet found application in nondestructive evaluation (NDE). This paper applies RTM imaging to NDE applications in bounded samples, where reverberation is always...

In order to produce an accurate depth migration of a seismic line, one needs an accurate velocity model. This model should reflect the lithology and lithostatic load that the rocks have experienced through burial and deformation. If we understand the geometry of these rocks, and their behaviour upon deformation, we can produce a more accurate velocity model. Performing palinspastic restoration ...