Imaging internal multiples from subsalt VSP data — Examples of target-oriented interferometry

Geophysics ◽  
2008 ◽  
Vol 73 (4) ◽  
pp. S157-S168 ◽  
Author(s):  
Ivan Vasconcelos ◽  
Roel Snieder ◽  
Brian Hornby

Seismic interferometry has become a technology of growing interest for imaging borehole seismic data. We demonstrate that interferometry of internal multiples can be used to image targets above a borehole receiver array. By internal multiples, we refer to all types of waves that scatter multiple times inside the model. These include, for instance, interbed, intrasalt, and water-bottom multiples as well as conversions among them. We use an interferometry technique that is based on representation theorems for perturbed media and targets the reconstruction of specific primary reflections from multiply reflected waves. In this interferometry approach, we rely on shot-domain wavenumber separation to select the directions of waves arriving at a given receiver. Using a numerical walkaway (WAW) VSP experiment recorded by a subsalt borehole receiver array in the Sigsbee salt model, we use the interference of internal multiples to image the salt structure from below. In this numerical example, the interferometric image that uses internal multiples reconstructs the bottom- and top-of-salt reflectors above the receiver array as well as the subsalt sediment structure between the array and the salt. Because of the limited source summation in this interferometry example, the interferometric images show artifact reflectors within the salt body. We apply this method to a field walkaway VSP from the Gulf of Mexico. With the field data, we demonstrate that the choice of shot-domain wavenumbers in the target-oriented interferometry procedure controls the wavenumbers in the output pseudoshot gathers. Target-oriented interferometric imaging from the 20-receiver array recovers the top-of-salt reflector that is consistent with surface seismic images. We present our results with both correlation-based and deconvolution-based interferometry.

Geophysics ◽  
1997 ◽  
Vol 62 (3) ◽  
pp. 723-729 ◽  
Author(s):  
Colin M. Sayers

Wide‐aperture walkaway vertical seismic profile (VSP) data acquired through transversely isotropic horizontal layers can be used to determine the P phase‐slowness surface, local to a receiver array in a borehole. In the presence of dip, errors in the slowness surface may occur if the medium is assumed to be layered horizontally. If the acquisition plane is oriented parallel to the dip direction, the derived slowness is too large for sources offset from the well in the down‐dip direction and too small for sources offset from the well in the up‐dip direction. For acquisition parallel to the strike of the layers, the recovery of the P phase‐slowness in the vicinity of the receiver array is excellent. It is therefore preferable to orient the walkaway VSP in the strike direction to estimate the anisotropic parameters of the medium in the vicinity of a receiver array. However, this may not be possible. If the dip direction of all layers has the same azimuth, the variation of walkaway traveltimes with azimuth has a simple form. This allows data from a single walkaway VSP extending both sides of a well to be inverted for the local anisotropic P phase‐slowness surface at the receivers even in the presence of dip. If data are acquired at more than one azimuth, the dip direction can be determined.


Geophysics ◽  
2013 ◽  
Vol 78 (3) ◽  
pp. B103-B110 ◽  
Author(s):  
Charles Hurich ◽  
Sharon Deemer

Seismic images are inherently directionally biased by the source-receiver geometry. This directional bias is particularly problematic for seismic imaging in hard rock terrains where structural dips may have any orientation with respect to the surface. We tested a technique for partially mitigating directional bias by combining surface and borehole seismic data and evaluated the results of a first field test of the technique. In this technique, surface data acquired using standard 2D acquisition procedures were combined with borehole data derived from a walk-away vertical seismic profile (VSP). The VSP data were transformed into the borehole datum using seismic interferometry. The interferometry created virtual shot records comprising sources and receivers in the borehole. The virtual shot records were then processed, using standard common midpoint techniques, resulting in an image from the borehole datum. The combination of the surface and borehole data increased the range of illumination angles resulting in seismic images that included reflections from structures with a wider range of dips than is available to surface profiling alone. The field test demonstrated that the surface and borehole data provide complementary information, which is more than either data set alone can provide. The test also verified the robustness of the virtual source technique even when the original VSP data are highly contaminated by high-amplitude tube waves. These results demonstrated that the combined imaging approach has significant potential for application in the polydeformed hard rock domains often encountered in minerals exploration.


2007 ◽  
Vol 26 (6) ◽  
pp. 760-763 ◽  
Author(s):  
Brian E. Hornby ◽  
Jianhua Yu

2010 ◽  
Vol 7 (1) ◽  
pp. 41-48 ◽  
Author(s):  
Bao-Li Wang ◽  
Guang-Ming Zhu ◽  
Jing-Huai Gao

Geophysics ◽  
2013 ◽  
Vol 78 (2) ◽  
pp. S93-S103 ◽  
Author(s):  
Jakob B. U. Haldorsen ◽  
W. Scott Leaney ◽  
Richard T. Coates ◽  
Steen A. Petersen ◽  
Helge Ivar Rutledal ◽  
...  

We evaluated a method for using 3C vertical seismic profile data to image acoustic interfaces located between the surface source and a downhole receiver array. The approach was based on simple concepts adapted from whole-earth seismology, in which observed compressional and shear wavefields are traced back to a common origin. However, unlike whole-earth and passive seismology, in which physical sources are imaged, we used the observed compressional and shear wavefields to image secondary sources (scatterers) situated between the surface source and the downhole receiver array. The algorithm consisted of the following steps: first, estimating the receiver compressional wavefield; second, using polarization to estimating the shear wavefield; third, deconvolving the shear wavefield using estimates of the source wavelet obtained from the direct compressional wave; fourth, the compressional and shear wavefields were back projected into the volume between the source and receivers; where, finally, an imaging condition was applied. When applied to rig-source VSP data acquired in an extended-reach horizontal well, this process was demonstrated to give images of formation features in the overburden, consistent with surface-seismic images obtained from the same area.


2019 ◽  
Vol 133 ◽  
pp. 01008
Author(s):  
Bartosz Gierlach ◽  
Tomasz Danek ◽  
Ayiaz Kaderali

Inversion of velocity parameters for the walkaway VSP data in a multilayered medium can be impeded by velocity gradients and anisotropy in some layers. A problem occurs if we compare velocities obtained from borehole seismic profiling which are equal to their vertical components with the velocities calculated with paths coming from far offsets where the horizontal component plays an important role, especially when the vertical gradient exists and the ray paths are curve-shaped. In this contribution we present the results of velocity model inversion for VSP data considering velocity gradient and elliptical anisotropy. The algorithm consists of two steps, optimization of velocity parameters and optimization of ray paths for the given model. Both procedures use the Nelder-Mead simplex method which finds local minima. Due to the character of optimization we performed also multistart analysis which can provide information about possible equivalences between parameters. Analysis was conducted for different parameterizations, in some cases allowing introduction of additional parameters: vertical gradient and elliptical anisotropy coefficient. The optimal model for a specific set of data is chosen with the help of Bayesian Information Criterion to balance complexity of model with quality of approximation of traveltimes.


Geothermics ◽  
2021 ◽  
Vol 94 ◽  
pp. 102095
Author(s):  
David Li ◽  
Lianjie Huang ◽  
Benxin Chi ◽  
Kai Gao ◽  
Clay Jones ◽  
...  

2013 ◽  
Vol 17 (6) ◽  
pp. 951-957
Author(s):  
Maxim I. Protasov ◽  
Maxim N. Dmitriev

Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. N29-N40
Author(s):  
Modeste Irakarama ◽  
Paul Cupillard ◽  
Guillaume Caumon ◽  
Paul Sava ◽  
Jonathan Edwards

Structural interpretation of seismic images can be highly subjective, especially in complex geologic settings. A single seismic image will often support multiple geologically valid interpretations. However, it is usually difficult to determine which of those interpretations are more likely than others. We have referred to this problem as structural model appraisal. We have developed the use of misfit functions to rank and appraise multiple interpretations of a given seismic image. Given a set of possible interpretations, we compute synthetic data for each structural interpretation, and then we compare these synthetic data against observed seismic data; this allows us to assign a data-misfit value to each structural interpretation. Our aim is to find data-misfit functions that enable a ranking of interpretations. To do so, we formalize the problem of appraising structural interpretations using seismic data and we derive a set of conditions to be satisfied by the data-misfit function for a successful appraisal. We investigate vertical seismic profiling (VSP) and surface seismic configurations. An application of the proposed method to a realistic synthetic model shows promising results for appraising structural interpretations using VSP data, provided that the target region is well-illuminated. However, we find appraising structural interpretations using surface seismic data to be more challenging, mainly due to the difficulty of computing phase-shift data misfits.


Sign in / Sign up

Export Citation Format

Share Document