2.5D multifocusing imaging of crooked-line seismic surveys

Geophysics ◽  
2021 ◽  
pp. 1-67
Author(s):  
Hossein Jodeiri Akbari Fam ◽  
Mostafa Naghizadeh ◽  
Oz Yilmaz

Two-dimensional seismic surveys often are conducted along crooked line traverses due to the inaccessibility of rugged terrains, logistical and environmental restrictions, and budget limitations. The crookedness of line traverses, irregular topography, and complex subsurface geology with steeply dipping and curved interfaces could adversely affect the signal-to-noise ratio of the data. The crooked-line geometry violates the assumption of a straight-line survey that is a basic principle behind the 2D multifocusing (MF) method and leads to crossline spread of midpoints. Additionally, the crooked-line geometry can give rise to potential pitfalls and artifacts, thus, leads to difficulties in imaging and velocity-depth model estimation. We develop a novel multifocusing algorithm for crooked-line seismic data and revise the traveltime equation accordingly to achieve better signal alignment before stacking. Specifically, we present a 2.5D multifocusing reflection traveltime equation, which explicitly takes into account the midpoint dispersion and cross-dip effects. The new formulation corrects for normal, inline, and crossline dip moveouts simultaneously, which is significantly more accurate than removing these effects sequentially. Applying NMO, DMO, and CDMO separately tends to result in significant errors, especially for large offsets. The 2.5D multifocusing method can perform automatically with a coherence-based global optimization search on data. We investigated the accuracy of the new formulation by testing it on different synthetic models and a real seismic data set. Applying the proposed approach to the real data led to a high-resolution seismic image with a significant quality improvement compared to the conventional method. Numerical tests show that the new formula can accurately focus the primary reflections at their correct location, remove anomalous dip-dependent velocities, and extract true dips from seismic data for structural interpretation. The proposed method efficiently projects and extracts valuable 3D structural information when applied to crooked-line seismic surveys.

Geophysics ◽  
2006 ◽  
Vol 71 (5) ◽  
pp. U67-U76 ◽  
Author(s):  
Robert J. Ferguson

The possibility of improving regularization/datuming of seismic data is investigated by treating wavefield extrapolation as an inversion problem. Weighted, damped least squares is then used to produce the regularized/datumed wavefield. Regularization/datuming is extremely costly because of computing the Hessian, so an efficient approximation is introduced. Approximation is achieved by computing a limited number of diagonals in the operators involved. Real and synthetic data examples demonstrate the utility of this approach. For synthetic data, regularization/datuming is demonstrated for large extrapolation distances using a highly irregular recording array. Without approximation, regularization/datuming returns a regularized wavefield with reduced operator artifacts when compared to a nonregularizing method such as generalized phase shift plus interpolation (PSPI). Approximate regularization/datuming returns a regularized wavefield for approximately two orders of magnitude less in cost; but it is dip limited, though in a controllable way, compared to the full method. The Foothills structural data set, a freely available data set from the Rocky Mountains of Canada, demonstrates application to real data. The data have highly irregular sampling along the shot coordinate, and they suffer from significant near-surface effects. Approximate regularization/datuming returns common receiver data that are superior in appearance compared to conventional datuming.


Geophysics ◽  
2016 ◽  
Vol 81 (4) ◽  
pp. U25-U38 ◽  
Author(s):  
Nuno V. da Silva ◽  
Andrew Ratcliffe ◽  
Vetle Vinje ◽  
Graham Conroy

Parameterization lies at the center of anisotropic full-waveform inversion (FWI) with multiparameter updates. This is because FWI aims to update the long and short wavelengths of the perturbations. Thus, it is important that the parameterization accommodates this. Recently, there has been an intensive effort to determine the optimal parameterization, centering the fundamental discussion mainly on the analysis of radiation patterns for each one of these parameterizations, and aiming to determine which is best suited for multiparameter inversion. We have developed a new parameterization in the scope of FWI, based on the concept of kinematically equivalent media, as originally proposed in other areas of seismic data analysis. Our analysis is also based on radiation patterns, as well as the relation between the perturbation of this set of parameters and perturbation in traveltime. The radiation pattern reveals that this parameterization combines some of the characteristics of parameterizations with one velocity and two Thomsen’s parameters and parameterizations using two velocities and one Thomsen’s parameter. The study of perturbation of traveltime with perturbation of model parameters shows that the new parameterization is less ambiguous when relating these quantities in comparison with other more commonly used parameterizations. We have concluded that our new parameterization is well-suited for inverting diving waves, which are of paramount importance to carry out practical FWI successfully. We have demonstrated that the new parameterization produces good inversion results with synthetic and real data examples. In the latter case of the real data example from the Central North Sea, the inverted models show good agreement with the geologic structures, leading to an improvement of the seismic image and flatness of the common image gathers.


Geophysics ◽  
2020 ◽  
Vol 85 (2) ◽  
pp. V223-V232 ◽  
Author(s):  
Zhicheng Geng ◽  
Xinming Wu ◽  
Sergey Fomel ◽  
Yangkang Chen

The seislet transform uses the wavelet-lifting scheme and local slopes to analyze the seismic data. In its definition, the designing of prediction operators specifically for seismic images and data is an important issue. We have developed a new formulation of the seislet transform based on the relative time (RT) attribute. This method uses the RT volume to construct multiscale prediction operators. With the new prediction operators, the seislet transform gets accelerated because distant traces get predicted directly. We apply our method to synthetic and real data to demonstrate that the new approach reduces computational cost and obtains excellent sparse representation on test data sets.


2019 ◽  
Vol 7 (3) ◽  
pp. T701-T711
Author(s):  
Jianhu Gao ◽  
Bingyang Liu ◽  
Shengjun Li ◽  
Hongqiu Wang

Hydrocarbon detection is always one of the most critical sections in geophysical exploration, which plays an important role in subsequent hydrocarbon production. However, due to the low signal-to-noise ratio and weak reflection amplitude of deep seismic data, some conventional methods do not always provide favorable hydrocarbon prediction results. The interesting dolomite reservoirs in Central Sichuan are buried over an average depth of 4500 m, and the dolomite rocks have a low porosity below approximately 4%, which is measured by well-logging data. Furthermore, the dominant system of pores and fractures as well as strong heterogeneity along the lateral and vertical directions lead to some difficulties in describing the reservoir distribution. Spectral decomposition (SD) has become successful in illuminating subsurface features and can also be used to identify potential hydrocarbon reservoirs by detecting low-frequency shadows. However, the current applications for hydrocarbon detection always suffer from low resolution for thin reservoirs, probably due to the influence of the window function and without a prior constraint. To address this issue, we developed sparse inverse SD (SISD) based on the wavelet transform, which involves a sparse constraint of time-frequency spectra. We focus on investigating the applications of sparse spectral attributes derived from SISD to deep marine dolomite hydrocarbon detection from a 3D real seismic data set with an area of approximately [Formula: see text]. We predict and evaluate gas-bearing zones in two target reservoir segments by analyzing and comparing the spectral amplitude responses of relatively high- and low-frequency components. The predicted results indicate that most favorable gas-bearing areas are located near the northeast fault zone in the upper reservoir segment and at the relatively high structural positions in the lower reservoir segment, which are in good agreement with the gas-testing results of three wells in the study area.


Geophysics ◽  
2019 ◽  
Vol 84 (6) ◽  
pp. R989-R1001 ◽  
Author(s):  
Oleg Ovcharenko ◽  
Vladimir Kazei ◽  
Mahesh Kalita ◽  
Daniel Peter ◽  
Tariq Alkhalifah

Low-frequency seismic data are crucial for convergence of full-waveform inversion (FWI) to reliable subsurface properties. However, it is challenging to acquire field data with an appropriate signal-to-noise ratio in the low-frequency part of the spectrum. We have extrapolated low-frequency data from the respective higher frequency components of the seismic wavefield by using deep learning. Through wavenumber analysis, we find that extrapolation per shot gather has broader applicability than per-trace extrapolation. We numerically simulate marine seismic surveys for random subsurface models and train a deep convolutional neural network to derive a mapping between high and low frequencies. The trained network is then tested on sections from the BP and SEAM Phase I benchmark models. Our results indicate that we are able to recover 0.25 Hz data from the 2 to 4.5 Hz frequencies. We also determine that the extrapolated data are accurate enough for FWI application.


Geophysics ◽  
1989 ◽  
Vol 54 (2) ◽  
pp. 181-190 ◽  
Author(s):  
Jakob B. U. Haldorsen ◽  
Paul A. Farmer

Occasionally, seismic data contain transient noise that can range from being a nuisance to becoming intolerable when several seismic vessels try simultaneously to collect data in an area. The traditional approach to solving this problem has been to allocate time slots to the different acquisition crews; the procedure, although effective, is very expensive. In this paper a statistical method called “trimmed mean stack” is evaluated as a tool for reducing the detrimental effects of noise from interfering seismic crews. Synthetic data, as well as field data, are used to illustrate the efficacy of the technique. Although a conventional stack gives a marginally better signal‐to‐noise ratio (S/N) for data without interference noise, typical usage of the trimmed mean stack gives a reduced S/N equivalent to a fold reduction of about 1 or 2 percent. On the other hand, for a data set containing high‐energy transient noise, trimming produces stacked sections without visible high‐amplitude contaminating energy. Equivalent sections produced with conventional processing techniques would be totally unacceptable. The application of a trimming procedure could mean a significant reduction in the costs of data acquisition by allowing several seismic crews to work simultaneously.


Geophysics ◽  
2010 ◽  
Vol 75 (4) ◽  
pp. V51-V60 ◽  
Author(s):  
Ramesh (Neelsh) Neelamani ◽  
Anatoly Baumstein ◽  
Warren S. Ross

We propose a complex-valued curvelet transform-based (CCT-based) algorithm that adaptively subtracts from seismic data those noises for which an approximate template is available. The CCT decomposes a geophysical data set in terms of small reflection pieces, with each piece having a different characteristic frequency, location, and dip. One can precisely change the amplitude and shift the location of each seismic reflection piece in a template by controlling the amplitude and phase of the template's CCT coefficients. Based on these insights, our approach uses the phase and amplitude of the data's and template's CCT coefficients to correct misalignment and amplitude errors in the noise template, thereby matching the adapted template with the actual noise in the seismic data, reflection event-by-event. We also extend our approach to subtract noises that require several templates to be approximated. By itself, the method can only correct small misalignment errors ([Formula: see text] in [Formula: see text] data) in the template; it relies on conventional least-squares (LS) adaptation to correct large-scale misalignment errors, such as wavelet mismatches and bulk shifts. Synthetic and real-data results illustrate that the CCT-based approach improves upon the LS approach and a curvelet-based approach described by Herrmann and Verschuur.


Geophysics ◽  
1992 ◽  
Vol 57 (12) ◽  
pp. 1623-1632 ◽  
Author(s):  
Richard E. Duren ◽  
Stan V. Morris

Null steering refers to the removal (or zeroing) of interferences at specified dips by creating receiving patterns with nulls that are aligned on the interferences. This type of beamforming is more effective than forming a simple crossline array and can be applied to both multistreamer and swath data for reducing out‐of‐plane interferences (sideswipe, boat interference, etc.) that corrupt two‐dimensional (2-D) data (the desired signal). Many beamforming techniques lead to signal cancellation when the interferences are correlated with the desired signal. However, a beamforming technique that has been developed is effective in the presence of signal correlated interferences. The technique can be effectively extended to prestack and poststack seismic data. The number of interferences and their dips are identified by a visual examination of the plotted data. This information can be used to design filters that are applied to the total data set. The resulting 2-D data set is free from the crossline interferences with the inline 2-D data remaining unaltered. Model and real data comparisons between null steering and simple crossline array summation show that: (1) null steering significantly attenuates crossline interference, and (2) 2-D inline data, masked by sideswipe, can be revealed once sideswipe is attenuated by null steering. The real data examples show the identification and effective attenuation of interferences that could easily be interpreted as inline 2-D data: (1) an apparent steeply dipping event, and (2) an apparent flat “bright spot.”


Geophysics ◽  
2016 ◽  
Vol 81 (1) ◽  
pp. E57-E68 ◽  
Author(s):  
Martin Panzner ◽  
Jan Petter Morten ◽  
Wiktor Waldemar Weibull ◽  
Børge Arntsen

Subbasalt imaging has gained significant interest in the last two decades, driven by the urge to better understand the geologic structures beneath volcanic layers, which can be up to several kilometers thick. This understanding is crucial for the development and risking of hydrocarbon play models in these areas. However, imaging based on the reflection seismic data alone suffers from severe amplitude transmission losses and interbed multiples in the volcanic sequence, as well as from poor definition of the subbasalt velocity structure. We have considered a sequential imaging workflow, in which the resistivity model from joint controlled-source electromagnetic and magnetotelluric data inversion was used to update the velocity model and to improve the structural definition in the migrated seismic image. The quantitative link between resistivity and velocity was derived from well data. The workflow used standard procedures for seismic velocity analysis, electromagnetic data inversion, and well analysis, and thereby allowed detail control and input based on additional geophysical knowledge and experience in each domain. Using real data sets from the Faroe-Shetland Basin, we can demonstrate that the integration of seismic and electromagnetic data significantly improved the imaging of geologic structures covered by up to several-kilometer-thick extended volcanic sequences. The improved results might alter the interpretation compared with the imaging results from seismic data alone.


Geophysics ◽  
2014 ◽  
Vol 79 (6) ◽  
pp. B243-B252 ◽  
Author(s):  
Peter Bergmann ◽  
Artem Kashubin ◽  
Monika Ivandic ◽  
Stefan Lüth ◽  
Christopher Juhlin

A method for static correction of time-lapse differences in reflection arrival times of time-lapse prestack seismic data is presented. These arrival-time differences are typically caused by changes in the near-surface velocities between the acquisitions and had a detrimental impact on time-lapse seismic imaging. Trace-to-trace time shifts of the data sets from different vintages are determined by crosscorrelations. The time shifts are decomposed in a surface-consistent manner, which yields static corrections that tie the repeat data to the baseline data. Hence, this approach implies that new refraction static corrections for the repeat data sets are unnecessary. The approach is demonstrated on a 4D seismic data set from the Ketzin [Formula: see text] pilot storage site, Germany, and is compared with the result of an initial processing that was based on separate refraction static corrections. It is shown that the time-lapse difference static correction approach reduces 4D noise more effectively than separate refraction static corrections and is significantly less labor intensive.


Sign in / Sign up

Export Citation Format

Share Document