scholarly journals Azimuthally Anisotropic 3D Velocity Continuation

2011 ◽  
Vol 2011 ◽  
pp. 1-8 ◽  
Author(s):  
William Burnett ◽  
Sergey Fomel

We extend time-domain velocity continuation to the zero-offset 3D azimuthally anisotropic case. Velocity continuation describes how a seismic image changes given a change in migration velocity. This description turns out to be of a wave propagation process, in which images change along a velocity axis. In the anisotropic case, the velocity model is multiparameter. Therefore, anisotropic image propagation is multidimensional. We use a three-parameter slowness model, which is related to azimuthal variations in velocity, as well as their principal directions. This information is useful for fracture and reservoir characterization from seismic data. We provide synthetic diffraction imaging examples to illustrate the concept and potential applications of azimuthal velocity continuation and to analyze the impulse response of the 3D velocity continuation operator.

2021 ◽  
Author(s):  
Farah Syazana Dzulkefli ◽  
Kefeng Xin ◽  
Ahmad Riza Ghazali ◽  
Guo Qiang ◽  
Tariq Alkhalifah

Abstract Salt is known for having a generally low density and higher velocity compared with the surrounding rock layers which causes the energy to scatter once the seismic wavefield hits the salt body and relatively less energy is transmitted through the salt to the deeper subsurface. As a result, most of imaging approaches are unable to image the base of the salt and the reservoir below the salt. Even the velocity model building such as FWI often fails to illuminate the deeper parts of salt area. In this paper, we show that Full Wavefield Redatuming (FWR) is used to retrieved and enhance the seismic data below the salt area, leading to a better seismic image quality and allowing us to focus on updating the velocity in target area below the salt. However, this redatuming approach requires a good overburden velocity model to retrieved good redatumed data. Thus, by using synthetic SEAM model, our objective is to study on the accuracy of the overburden velocity model required for imaging beneath complex overburden. The results show that the kinematic components of wave propagation are preserved through redatuming even with heavily smoothed overburden velocity model.


Geophysics ◽  
1996 ◽  
Vol 61 (6) ◽  
pp. 1846-1858 ◽  
Author(s):  
Claudio Bagaini ◽  
Umberto Spagnolini

Continuation to zero offset [better known as dip moveout (DMO)] is a standard tool for seismic data processing. In this paper, the concept of DMO is extended by introducing a set of operators: the continuation operators. These operators, which are implemented in integral form with a defined amplitude distribution, perform the mapping between common shot or common offset gathers for a given velocity model. The application of the shot continuation operator for dip‐independent velocity analysis allows a direct implementation in the acquisition domain by exploiting the comparison between real data and data continued in the shot domain. Shot and offset continuation allow the restoration of missing shot or missing offset by using a velocity model provided by common shot velocity analysis or another dip‐independent velocity analysis method.


2004 ◽  
Vol 44 (1) ◽  
pp. 123 ◽  
Author(s):  
G.P. Thomas ◽  
M.R. Lennane ◽  
F. Glass ◽  
T. Walker ◽  
M. Partington ◽  
...  

The eastern Dampier Sub-basin on Australia’s northwestern margin has been subject to intensive exploration activity since the early 1960s. The commercial success rate for exploration drilling, however, has been a disappointing 8%, despite numerous indications of at least one active petroleum system. During 2002–2003, Woodside and its joint venture partners undertook an integrated review of the area, aimed at unlocking its remaining potential. Stratigraphy, hydrocarbon charge and 3D seismic data quality were addressed in parallel.The eastern Dampier Sub-basin stratigraphy was upgraded from the existing, conventional, second-order tectono-stratigraphic framework to a third-order, exploration-scale, genetic stratigraphic framework. The new framework has regional predictive capability in terms of reservoir (and seal) presence and facies, and has led to recognition of new plays and an enhanced understanding of known plays. One new play involves shoreface sands within the Calypso Formation. New light has been shed on the known Lower Cretaceous M.australis sands play (K30), by the creation of gross depositional environment maps at third-order sequence scale. The Upper Jurassic deepwater clastics play of the Lewis Trough has also been developed, by recognition of four prospective, sand-rich gravity-flow intervals in the early Oxfordian (J42 play).A 3D charge modelling study, underpinned by new geochemical analysis, has allowed delineation of areas of higher and lower risk in terms of hydrocarbon charge and phase (oil versus gas). Key source rocks for oil are identified in the early Oxfordian W.spectabilis biozone, although they are also a likely source for gas in the southwest of the area. The Bathonian-Callovian Upper Legendre Formation is a major source for gas, but could also have contributed minor oil in the northeast of the area. By a combination of geochemical fingerprinting and 3D forward modelling, most hydrocarbon occurrences in the area have been tied to these source intervals, complete with a consistent view of maturities and migration pathways.Some 1,500 km2 of the Panaeus multi-client 3D survey were reprocessed, with close attention to multiple removal, velocities and imaging. A step-change improvement in seismic quality was obtained, together with improved velocities for depth conversion.The prospect portfolio has been polarised and much enhanced through these studies, and the results of several existing wells have become better understood. Some new prospects were identified by apparent direct fluid indications, detected in one case by 3D volume AVO screening. Other new prospects are the result of a clearer seismic image, or of the revised velocity model for depth conversion. New plays are still being followed up, while the fresh light cast on existing plays (e.g. K30 and J42), in combination with improved seismic data, has led to development of several interesting opportunities.


2005 ◽  
Vol 45 (1) ◽  
pp. 421
Author(s):  
P. Bocca ◽  
L. Fava ◽  
E. Stolf

3D pre-stack depth migration (PSDM) reprocessing was conducted in 2003 on a portion of the Onnia 3D seismic cube, located in exploration permit AC/P-21, Timor Sea.The main objective of the reprocessing was to obtain the best seismic depth image and the most realistic structural reconstruction of the sub-surface to mitigate the risk factors associated with trap definition (trap retention and trap efficiency). This represents one of the main challenges for oil exploration in the area.The 3D PSDM methodology was chosen as the most appropriate imaging tool to define the correct sub-surface geometry and fault imaging through the use of an appropriate velocity field. An integrated approach to building the final velocity model was adopted, with a substantial contribution from the regional geological model.Several examples are given to demonstrate that the 3D PSDM reprocessing significantly improved the seismic image and thus the confidence in the interpretation, contributing strongly to the definition of the exploration targets.The interpretation of the new seismic data has resulted in a new structural picture in which higher confidence in seismic imaging has improved fault correlation. This has enabled better structural definition at the Middle Jurassic Plover Formation level that has reduced the complexity of the large Vesta Prospect, in the centre of the Swan Graben to the northwest of East Swan–1. Improved understanding of the fault reactivation mechanism and the structural elements of the trap (trap integrity) were eventually incorporated in the prospect risking.In the Swan Graben 3D PSDM has proved to be a very powerful instrument capable of producing significant impact on the exploration even in an area with a complex geological setting and a fairly poor seismic data quality.


Geophysics ◽  
1994 ◽  
Vol 59 (3) ◽  
pp. 351-361 ◽  
Author(s):  
M. Ali C. Tura ◽  
Robert J. Greaves ◽  
Wafik B. Beydoun

A crosswell seismic experiment at the San Emidio oil field in Bakersfield, California, is carried out to evaluate crosswell reflection/diffraction tomography and image the interwell region to locate a possible pinchout zone. In this experiment, the two wells used are 2500 ft (762 m) apart, and the zone to be imaged is 11 000 ft (3350 m) to 13 000 ft (3960 m) deep. With the considered distances, this experiment forms the first large scale reservoir characterization application of crosswell reflection/diffraction tomography. A subset of the intended data, formed of two common receiver gathers and one common shot gather, was collected at the San Emidio oil field. The crosswell data display a wide variety of wave modes including tube waves, singly and multiply reflected/diffracted waves, and refracted waves. The data are processed using frequency filters, median filters, and spatial muting filters to enhance the reflected/diffracted energy. A 2-D layered velocity model with gradients is built using zero‐offset VSPs and full‐waveform acoustic logs from the two wells. This model is used to generate synthetic finite‐difference data for the field data acquisition geometry. The synthetic data are processed and imaged using the elastic ray‐Born 𝓁2-migration/inversion (ERBMI) method. A smooth 2-D velocity model incorporating only gradients and a few layers is used as a background model for the imaging. Considering the limited data acquisition geometry, synthetic data images compare favorably with the initial velocity model. With the encouraging results obtained from synthetic data, the ERBMI method, with the smooth background velocity model is used next to image the processed field data. Images obtained from the crosswell data show a good match with the reflected field in the zero‐offset VSPs and with migrated surface seismic data. From the interpretation of these images, the potential of this crosswell seismic method for answering questions regarding reservoir continuity and existence of pinchout zones can be seen.


2021 ◽  
Author(s):  
Alexander Bauer ◽  
Benjamin Schwarz ◽  
Dirk Gajewski

<p>Most established methods for the estimation of subsurface velocity models rely on the measurements of reflected or diving waves and therefore require data with sufficiently large source-receiver offsets. For seismic data that lacks these offsets, such as vintage data, low-fold academic data or near zero-offset P-Cable data, these methods fail. Building on recent studies, we apply a workflow that exploits the diffracted wavefield for depth-velocity-model building. This workflow consists of three principal steps: (1) revealing the diffracted wavefield by modeling and adaptively subtracting reflections from the raw data, (2) characterizing the diffractions with physically meaningful wavefront attributes, (3) estimating depth-velocity models with wavefront tomography. We propose a hybrid 2D/3D approach, in which we apply the well-established and automated 2D workflow to numerous inlines of a high-resolution 3D P-Cable dataset acquired near Ritter Island, a small volcanic island located north-east of New Guinea known for a catastrophic flank collapse in 1888. We use the obtained set of parallel 2D velocity models to interpolate a 3D velocity model for the whole data cube, thus overcoming possible issues such as varying data quality in inline and crossline direction and the high computational cost of 3D data analysis. Even though the 2D workflow may suffer from out-of-plane effects, we obtain a smooth 3D velocity model that is consistent with the data.</p>


2017 ◽  
Vol 5 (3) ◽  
pp. T279-T285 ◽  
Author(s):  
Parvaneh Karimi ◽  
Sergey Fomel ◽  
Rui Zhang

Integration of well-log data and seismic data to predict rock properties is an essential but challenging task in reservoir characterization. The standard methods commonly used to create subsurface model do not fully honor the importance of seismic reflectors and detailed structural information in guiding the spatial distribution of rock properties in the presence of complex structures, which can make these methods inaccurate. To overcome initial model accuracy limitations in structurally complex regimes, we have developed a method that uses the seismic image structures to accurately constrain the interpolation of well properties between well locations. A geologically consistent framework provides a more robust initial model that, when inverted with seismic data, delivers a highly detailed yet accurate subsurface model. An application to field data from the North Sea demonstrates the effectiveness of our method, which proves that incorporating the seismic structural framework when interpolating rock properties between wells culminates in the increased accuracy of the final inverted result compared with the standard inversion workflows.


Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE269-VE280 ◽  
Author(s):  
Priyank Jaiswal ◽  
Colin A. Zelt

Imaging 2D multichannel land seismic data can be accomplished effectively by a combination of traveltime inversion and prestack depth migration (PSDM), referred to as unified imaging. Unified imaging begins by inverting the direct-arrival times to estimate a velocity model that is used in static corrections and stacking velocity analysis. The interval velocity model (from stacking velocities) is used for PSDM. The stacked data and the PSDM image are interpreted for common horizons, and the corresponding wide-aperture reflections are identified in the shot gathers. Using the interval velocity model, the stack interpretations are inverted as zero-offset reflections to constrain the corresponding interfaces in depth; the interval velocity model remains stationary. We define a coefficient of congruence [Formula: see text] that measures the discrepancy between horizons from the PSDM image andtheir counterparts from the zero-offset inversion. A value of unity for [Formula: see text] implies that the interpreted and inverted horizons are consistent to within the interpretational uncertainties, and the unified imaging is said to have converged. For [Formula: see text] greater than unity, the interval velocity model and the horizon depths are updated by jointly inverting the direct arrivals with the zero-offset and wide-aperture reflections. The updated interval velocity model is used again for both PSDM and a zero-offset inversion. Interpretations of the new PSDM image are the updated horizon depths. The unified imaging is applied to seismic data from the Naga Thrust and Fold Belt in India. Wide-aperture and zero-offset data from three geologically significant horizons are used. Three runs of joint inversion and PSDM are required in a cyclic manner for [Formula: see text] to converge to unity. A joint interpretation of the final velocity model and depth image reveals the presence of a triangle zone that could be promising for exploration.


2010 ◽  
Vol 50 (2) ◽  
pp. 723
Author(s):  
Sergey Birdus ◽  
Erika Angerer ◽  
Iftikhar Abassi

Processing of multi and wide-azimuth seismic data faces some new challenges, and one of them is depth-velocity modelling and imaging with azimuthal velocity anisotropy. Analysis of multi-azimuth data very often reveals noticeable fluctuations in moveout between different acquisition directions. They can be caused by several factors: real azimuthal interval velocity anisotropy associated with quasi-vertical fractures or present day stress field within the sediments; short-wavelength velocity heterogeneities in the overburden; TTI (or VTI) anisotropy in the overburden; or, random distortions due to noise, multiples, irregularities in the acquisition geometry, etcetera. In order to build a velocity model for multi-azimuth pre-stack depth migration (MAZ PSDM) taking into account observed azimuthal anisotropy, we need to recognise, separate and estimate all the effects listed above during iterative depth-velocity modelling. Analysis of seismic data from a full azimuth 3D seismic land survey revealed the presence of strong spatially variable azimuthal velocity anisotropy that had to be taken into consideration. Using real data examples we discuss major steps in depth processing workflow that took such anisotropy into account: residual moveout estimation in azimuth sectors; separation of different effects causing apparent azimuthal anisotropy (see A–D above); iterative depth-velocity modelling with azimuthal anisotropy; and, subsequent MAZ anisotropic PSDM. The presented workflow solved problems with azimuthal anisotropy in our multi-azimuth dataset. Some of the lessons learned during this MAZ project are relevant to every standard narrow azimuth seismic survey recorded in complex geological settings.


2020 ◽  
Author(s):  
Alexander Bauer ◽  
Benjamin Schwarz ◽  
Richard Delf ◽  
Dirk Gajewski

<p>In the recent years, the diffracted wavefield has gained increasing attention in the field of applied seismics. While classical seismic imaging and inversion schemes mainly focus on high-amplitude reflected measurements, the faint and often masked diffracted wavefield is neglected or even treated as noise. In order to be able to extract depth-velocity models from seismic reflection data, sufficiently large source-receiver offsets are needed. However, the acquisition of such multi-channel seismic data is expensive and often only feasible for the hydrocarbon industry, while academia has to cope with low-fold or zero-offset data. The diffracted wavefield is the key for extracting depth velocities from such data, as the moveout of diffractions – in contrast to reflections – can be measured in the zero-offset domain. Recently, we have demonstrated on multi-channel, single-channel and passive seismic data that by means of wavefront tomography depth-velocity models can be retrieved solely based on diffractions or passive seismic events along with the localizations of these scatterers. The input for wavefront tomography are so-called wavefront attributes, which can be extracted from the data in an unsupervised fashion by means of coherence analysis. In order to obtain the required diffraction-only data, we use a recently proposed scheme that adaptively subtracts the high-amplitude reflected wavefield from the raw data. Due to their most common acquisition geometry, most ground-penetrating-radar (GPR) data inherently lack offsets. In addition, GPR data generally contain a rich diffracted wavefield, which in turn contains information about sought-after structures, as diffractions are caused by small-scale heterogeneities such as faults, tips or edges. In this work, we show an application of the suggested workflow – coherence analysis, diffraction separation and diffraction wavefront tomography – to GPR data acquired at a glacier, resulting in a depth-velocity model and the localizations of the scatterers, both obtained in a fully unsupervised fashion. While the resulting  velocity model may be used for depth migration of the raw data, the localizations of the scatterers may in addition provide important information on the inner structure of the glacier in order to, for instance, localize water intrusions or fractures.</p>


Sign in / Sign up

Export Citation Format

Share Document