Regularizing 3-D data sets with DMO

Geophysics ◽  
1996 ◽  
Vol 61 (4) ◽  
pp. 1103-1114 ◽  
Author(s):  
Anat Canning ◽  
Gerald H. F. Gardner
Keyword(s):  
Computer Program ◽  
Seismic Data ◽  
Physical Modeling ◽  
Data Sets ◽  
Angle Of Incidence ◽  
Data Set ◽  
Velocity Function ◽  
Inherent Problem ◽  
Original Acquisition ◽  
Good Agreement

The combination of DMO and [Formula: see text] is used here to change the original acquisition geometry of a 3-D seismic data set into a more convenient form. For example, irregular 3-D surveys can be projected onto a regular midpoint‐offset grid with zero source‐receiver azimuth and equal increments in offset. The algorithm presented here is based on a new, nonaliased 3-D DMO algorithm in (f, x) domain. It does not require any knowledge of the velocity function for constant or rms velocity variations. The computer program was designed to process and to output very large multifold 3-D data sets. A synthetic example of a point diffractor in 3-D space and a 3-D experiment in a physical modeling tank are used to demonstrate the procedure. In both cases, the results obtained after the data set is regularized are compared with a data set that was acquired initially with the desired configuration. These comparisons show very good agreement. Analysis of the procedure indicates that it may not reconstruct AVO correctly. This is an inherent problem that occurs because the reorganization procedure changes the angle of incidence.

Sparse reflectivity inversion for nonstationary seismic data with surface-related multiples: Numerical and field-data experiments

Geophysics ◽  
2017 ◽  
Vol 82 (3) ◽  
pp. R199-R217 ◽  
Author(s):  
Xintao Chai ◽  
Shangxu Wang ◽  
Genyang Tang
Keyword(s):  
Seismic Data ◽  
Resolution Enhancement ◽  
Synthetic Data ◽  
Data Sets ◽  
Data Set ◽  
Anelastic Attenuation ◽  
Seismic Resolution ◽  
Text Filtering ◽  
The Stability ◽  
Reflectivity Inversion

Seismic data are nonstationary due to subsurface anelastic attenuation and dispersion effects. These effects, also referred to as the earth’s [Formula: see text]-filtering effects, can diminish seismic resolution. We previously developed a method of nonstationary sparse reflectivity inversion (NSRI) for resolution enhancement, which avoids the intrinsic instability associated with inverse [Formula: see text] filtering and generates superior [Formula: see text] compensation results. Applying NSRI to data sets that contain multiples (addressing surface-related multiples only) requires a demultiple preprocessing step because NSRI cannot distinguish primaries from multiples and will treat them as interference convolved with incorrect [Formula: see text] values. However, multiples contain information about subsurface properties. To use information carried by multiples, with the feedback model and NSRI theory, we adapt NSRI to the context of nonstationary seismic data with surface-related multiples. Consequently, not only are the benefits of NSRI (e.g., circumventing the intrinsic instability associated with inverse [Formula: see text] filtering) extended, but also multiples are considered. Our method is limited to be a 1D implementation. Theoretical and numerical analyses verify that given a wavelet, the input [Formula: see text] values primarily affect the inverted reflectivities and exert little effect on the estimated multiples; i.e., multiple estimation need not consider [Formula: see text] filtering effects explicitly. However, there are benefits for NSRI considering multiples. The periodicity and amplitude of the multiples imply the position of the reflectivities and amplitude of the wavelet. Multiples assist in overcoming scaling and shifting ambiguities of conventional problems in which multiples are not considered. Experiments using a 1D algorithm on a synthetic data set, the publicly available Pluto 1.5 data set, and a marine data set support the aforementioned findings and reveal the stability, capabilities, and limitations of the proposed method.

Common-focus point-based target-oriented imaging approach for continuous seismic reservoir monitoring

Geophysics ◽  
2018 ◽  
Vol 83 (4) ◽  
pp. M41-M48 ◽  
Author(s):  
Hongwei Liu ◽  
Mustafa Naser Al-Ali
Keyword(s):  
Seismic Data ◽  
Rock Physics ◽  
Time Lapse ◽  
Velocity Estimation ◽  
Estimation Methods ◽  
Data Sets ◽  
Data Set ◽  
Velocity Models ◽  
Reservoir Monitoring ◽  
Focus Point

The ideal approach for continuous reservoir monitoring allows generation of fast and accurate images to cope with the massive data sets acquired for such a task. Conventionally, rigorous depth-oriented velocity-estimation methods are performed to produce sufficiently accurate velocity models. Unlike the traditional way, the target-oriented imaging technology based on the common-focus point (CFP) theory can be an alternative for continuous reservoir monitoring. The solution is based on a robust data-driven iterative operator updating strategy without deriving a detailed velocity model. The same focusing operator is applied on successive 3D seismic data sets for the first time to generate efficient and accurate 4D target-oriented seismic stacked images from time-lapse field seismic data sets acquired in a [Formula: see text] injection project in Saudi Arabia. Using the focusing operator, target-oriented prestack angle domain common-image gathers (ADCIGs) could be derived to perform amplitude-versus-angle analysis. To preserve the amplitude information in the ADCIGs, an amplitude-balancing factor is applied by embedding a synthetic data set using the real acquisition geometry to remove the geometry imprint artifact. Applying the CFP-based target-oriented imaging to time-lapse data sets revealed changes at the reservoir level in the poststack and prestack time-lapse signals, which is consistent with the [Formula: see text] injection history and rock physics.

Spectral analysis and inversion of experimental codas

Geophysics ◽  
1993 ◽  
Vol 58 (3) ◽  
pp. 408-418 ◽  
Author(s):  
L. R. Jannaud ◽  
P. M. Adler ◽  
C. G. Jacquin
Keyword(s):  
Spectral Analysis ◽  
Characteristic Length ◽  
Gaussian Model ◽  
Seismic Survey ◽  
Data Sets ◽  
Data Set ◽  
Anisotropic Elastic ◽  
And Inversion ◽  
Good Agreement

A method developed for the determination of the characteristic lengths of an heterogeneous medium from the spectral analysis of codas is based on an extension of Aki’s theory to anisotropic elastic media. An equivalent Gaussian model is obtained and seems to be in good agreement with the two experimental data sets that illustrate the method. The first set was obtained in a laboratory experiment with an isotropic marble sample. This sample is characterized by a submillimetric length scale that can be directly observed on a thin section. The spectral analysis of codas and their inversion yields an equivalent correlation length that is in good agreement with the observed one. The second data set is obtained in a crosshole experiment at the usual scale of a seismic survey. The codas are recorded, analysed, and inverted. The analysis yields a vertical characteristic length for the studied subsurface that compares well with the characteristic length measured by seismic and stratigraphic logs.

A robust workflow for performing joint impedance inversion and its applications

2020 ◽  
Vol 8 (1) ◽  
pp. T141-T149
Author(s):  
Ritesh Kumar Sharma ◽  
Satinder Chopra ◽  
Larry R. Lines
Keyword(s):  
Time Domain ◽  
Seismic Data ◽  
Sedimentary Basin ◽  
Vertical Component ◽  
Data Sets ◽  
Data Set ◽  
Matching Process ◽  
Multicomponent Seismic ◽  
Impedance Inversion ◽  
Straightforward Approach

Multicomponent seismic data offer several advantages for characterizing reservoirs with the use of the vertical component (PP) and mode-converted (PS) data. Joint impedance inversion inverts both of these data sets simultaneously; hence, it is considered superior to simultaneous impedance inversion. However, the success of joint impedance inversion depends on how accurately the PS data are mapped on the PP time domain. Normally, this is attempted by performing well-to-seismic ties for PP and PS data sets and matching different horizons picked on PP and PS data. Although it seems to be a straightforward approach, there are a few issues associated with it. One of them is the lower resolution of the PS data compared with the PP data that presents difficulties in the correlation of the equivalent reflection events on both the data sets. Even after a few consistent horizons get tracked, the horizon matching process introduces some artifacts on the PS data when mapped into PP time. We have evaluated such challenges using a data set from the Western Canadian Sedimentary Basin and then develop a novel workflow for addressing them. The importance of our workflow was determined by comparing data examples generated with and without its adoption.

scholarly journals Global validation of improved SCIAMACHY scientific ozone limb data using ozonesonde measurements

2015 ◽  
Vol 8 (5) ◽  
pp. 4817-4858
Author(s):  
J. Jia ◽  
A. Rozanov ◽  
A. Ladstätter-Weißenmayer ◽  
J. P. Burrows
Keyword(s):  
Global Scale ◽  
Data Sets ◽  
Vertical Profiles ◽  
Data Set ◽  
Latitude Range ◽  
Ozone Profile ◽  
The Tropics ◽  
Relative Differences ◽  
Significant Underestimation ◽  
Good Agreement

Abstract. In this manuscript, the latest SCIAMACHY limb ozone scientific vertical profiles, namely the current V2.9 and the upcoming V3.0, are extensively compared with ozone sonde data from the WOUDC database. The comparisons are made on a global scale from 2003 to 2011, involving 61 sonde stations. The retrieval processors used to generate V2.9 and V3.0 data sets are briefly introduced. The comparisons are discussed in terms of vertical profiles and stratospheric partial columns. Our results indicate that the V2.9 ozone profile data between 20–30 km is in good agreement with ground based measurements with less than 5% relative differences in the latitude range of 90° S–40° N (with exception of the tropical Pacific region where an overestimation of more than 10% is observed), which corresponds to less than 5 DU partial column differences. In the tropics the differences are within 3%. However, this data set shows a significant underestimation northwards of 40° N (up to ~15%). The newly developed V3.0 data set reduces this bias to below 10% while maintaining a good agreement southwards of 40° N with slightly increased relative differences of up to 5% in the tropics.

3D ray+Born migration/inversion—Part 2: Application to the SEG/EAGE overthrust experiment

Geophysics ◽  
2003 ◽  
Vol 68 (4) ◽  
pp. 1357-1370 ◽  
Author(s):  
Stéphane Operto ◽  
Gilles Lambaré ◽  
Pascal Podvin ◽  
Philippe Thierry
Keyword(s):  
Spatial Resolution ◽  
Seismic Data ◽  
Velocity Model ◽  
Nonuniform Distribution ◽  
Data Sets ◽  
Data Set ◽  
Quantitative Estimates ◽  
3D Processing ◽  
Type Data ◽  
New Formula

The SEG/EAGE overthrust model is a synthetic onshore velocity model that was used to generate several large synthetic seismic data sets using acoustic finite‐difference modeling. From this database, several realistic subdata sets were extracted and made available for testing 3D processing methods. For example, classic onshore‐type data‐acquisition geometries are available such as a swath acquisition, which is characterized by a nonuniform distribution of long offsets with azimuth and midpoints. In this paper, we present an application of 2.5D and 3D ray‐Born migration/inversion to several classical data sets from the SEG/EAGE overthrust experiment. The method is formulated as a linearized inversion of the scattered wavefield. The method allows quantitative estimates of short wavelength components of the velocity model. First, we apply a 3D migration/inversion formula formerly developed for marine acquisitions to the swath data set. The migrated sections exhibit significant amplitude artifacts and acquisition footprints, also revealed by the shape of the local spatial resolution filters. From the analysis of these spatial resolution filters, we propose a new formula significantly improving the migrated dip section. We also present 3D migrated results for the strike section and a small 3D target containing a channel. Finally, the applications demonstrate, that the ray+Born migration formula must be adapted to the acquisition geometry to obtain reliable estimates of the true amplitude of the model perturbations. This adaptation is relatively straightforward in the frame of the ray+Born formalism and can be guided by the analysis of the resolution operator.

scholarly journals TRANSFORMATION OF RANDOM BREAKING WAVES ON SURF BEAT

1986 ◽  
Vol 1 (20) ◽  
pp. 9 ◽  
Author(s):  
William R. Dally ◽  
Robert G. Dean
Keyword(s):  
Surf Zone ◽  
Wave Model ◽  
Breaking Waves ◽  
Wave Transformation ◽  
Random Wave ◽  
Data Sets ◽  
Data Set ◽  
Wave Decay ◽  
Mean Water Level ◽  
Good Agreement

Based on a previous study by the authors of regular breaking waves in the surf zone, a model for random wave transformation across the nearshore region is developed. The results of a laboratory investigation of the effect of a steady opposing current on the wave decay process are presented and a proposed governing equation verified. Surf beat effects on wave transformation are then included in the model by representing the long wave as a temporally and spatiallyvarying current and mean water level. The concept of an equivalent water depth, which contains the effect of the current, is introduced and then included in a stochastic form in the random wave model. Surf beat is found to noticeably increase the decay of the root mean square wave height, especially in the inner surf where the beat is strongest. Comparison of the models to two field data sets show very good agreement for Hotta and Mizuguchi (1980), but rather poor for Thornton and Guza (1983). Possible explanations for the unexpected behavior of the second data set, pertaining to filtering, are discussed. Finally, a possible explanation for the dependence of random wave decay on deepwater steepness, noted by Battjes and Stive (1985), is presented.

scholarly journals Comparison of 14C Ages Between LSC and AMS Measurements of Choukai Jindai Cedar Tree Rings at 2600 cal BP

Radiocarbon ◽  
2010 ◽  
Vol 52 (3) ◽  
pp. 895-900 ◽  
Author(s):  
Yui Takahashi ◽  
Hirohisa Sakurai ◽  
Kayo Suzuki ◽  
Taiichi Sato ◽  
Shuichi Gunji ◽  
...  
Keyword(s):  
Mass Spectrometry ◽  
Tree Rings ◽  
Liquid Scintillation ◽  
Measurement Methods ◽  
Scintillation Counting ◽  
Age Difference ◽  
Data Sets ◽  
Data Set ◽  
The Difference ◽  
Good Agreement

Radiocarbon ages of Choukai Jindai cedar tree rings growing in the excess era of 14C concentrations during 2757–2437 cal BP were measured using 2 types of 14C measurement methods, i.e. liquid scintillation counting (LSC) and accelerator mass spectrometry (AMS). The difference between the 2 methods is 3.7 ± 5.2 14C yr on average for 61 single-year tree rings, indicating good agreement between the methods. The Choukai data sets show a small sharp bump with an average 14C age of 2497.1 ± 3.0 14C yr BP during 2650–2600 cal BP. Although the profile of the Choukai LSC data set compares well with that of IntCal04, having a 14C age difference of 4.6 ± 5.3 14C yr on average, the Choukai LSC 14C ages indicate variability against the smoothed profile of IntCal04.

Land FWI in the Delaware Basin, west Texas: A case study

2020 ◽  
Vol 39 (5) ◽  
pp. 324-331
Author(s):  
Gary Murphy ◽  
Vanessa Brown ◽  
Denes Vigh
Keyword(s):  
Seismic Data ◽  
Waveform Inversion ◽  
Surface Velocity ◽  
High Signal ◽  
Data Sets ◽  
Data Set ◽  
Low Frequencies ◽  
Near Surface ◽  
Delaware Basin ◽  
West Texas

As part of a wide-reaching full-waveform inversion (FWI) research program, FWI is applied to an onshore seismic data set collected in the Delaware Basin, west Texas. FWI is routinely applied on typical marine data sets with high signal-to-noise ratio (S/N), relatively good low-frequency content, and reasonably long offsets. Land seismic data sets, in comparison, present significant challenges for FWI due to low S/N, a dearth of low frequencies, and limited offsets. Recent advancements in FWI overcome limitations due to poor S/N and low frequencies making land FWI feasible to use to update the shallow velocities. The chosen area has contrasting and variable near-surface conditions providing an excellent test data set on which to demonstrate the workflow and its challenges. An acoustic FWI workflow is used to update the near-surface velocity model in order to improve the deeper image and simultaneously help highlight potential shallow drilling hazards.

Time-lapse acoustic impedance variations during CO2 injection in Weyburn oilfield, Canada

Geophysics ◽  
2019 ◽  
Vol 85 (1) ◽  
pp. M1-M13 ◽  
Author(s):  
Yichuan Wang ◽  
Igor B. Morozov
Keyword(s):  
Oil Recovery ◽  
Seismic Data ◽  
Acoustic Impedance ◽  
Time Lapse ◽  
Time Shift ◽  
Co2 Injection ◽  
Inversion Method ◽  
Data Sets ◽  
Storage Site ◽  
Data Set

For seismic monitoring injected fluids during enhanced oil recovery or geologic [Formula: see text] sequestration, it is useful to measure time-lapse (TL) variations of acoustic impedance (AI). AI gives direct connections to the mechanical and fluid-related properties of the reservoir or [Formula: see text] storage site; however, evaluation of its subtle TL variations is complicated by the low-frequency and scaling uncertainties of this attribute. We have developed three enhancements of TL AI analysis to resolve these issues. First, following waveform calibration (cross-equalization) of the monitor seismic data sets to the baseline one, the reflectivity difference was evaluated from the attributes measured during the calibration. Second, a robust approach to AI inversion was applied to the baseline data set, based on calibration of the records by using the well-log data and spatially variant stacking and interval velocities derived during seismic data processing. This inversion method is straightforward and does not require subjective selections of parameterization and regularization schemes. Unlike joint or statistical inverse approaches, this method does not require prior models and produces accurate fitting of the observed reflectivity. Third, the TL AI difference is obtained directly from the baseline AI and reflectivity difference but without the uncertainty-prone subtraction of AI volumes from different seismic vintages. The above approaches are applied to TL data sets from the Weyburn [Formula: see text] sequestration project in southern Saskatchewan, Canada. High-quality baseline and TL AI-difference volumes are obtained. TL variations within the reservoir zone are observed in the calibration time-shift, reflectivity-difference, and AI-difference images, which are interpreted as being related to the [Formula: see text] injection.

Sign in / Sign up

Export Citation Format

Share Document