Near surface velocity investigation with STEM and its application on the static correction for seismic data

2009 ◽  
Author(s):  
Liangjun Yan ◽  
Zhuliu Su ◽  
Xingbing Xie
Geophysics ◽  
2020 ◽  
pp. 1-57
Author(s):  
Daniele Colombo ◽  
Ernesto Sandoval ◽  
Diego Rovetta ◽  
Apostolos Kontakis

Land seismic velocity modeling is a difficult task largely related to the description of the near surface complexities. Full waveform inversion is the method of choice for achieving high-resolution velocity mapping but its application to land seismic data faces difficulties related to complex physics, unknown and spatially varying source signatures, and low signal-to-noise ratio in the data. Large parameter variations occur in the near surface at various scales causing severe kinematic and dynamic distortions of the recorded wavefield. Some of the parameters can be incorporated in the inversion model while others, due to sub-resolution dimensions or unmodeled physics, need to be corrected through data preconditioning; a topic not well described for land data full waveform inversion applications. We have developed novel algorithms and workflows for surface-consistent data preconditioning utilizing the transmitted portion of the wavefield, signal-to-noise enhancement by generation of CMP-based virtual super shot gathers, and robust 1.5D Laplace-Fourier full waveform inversion. Our surface-consistent scheme solves residual kinematic corrections and amplitude anomalies via scalar compensation or deconvolution of the near surface response. Signal-to-noise enhancement is obtained through the statistical evaluation of volumetric prestack responses at the CMP position, or virtual super (shot) gathers. These are inverted via a novel 1.5D acoustic Laplace-Fourier full waveform inversion scheme using the Helmholtz wave equation and Hankel domain forward modeling. Inversion is performed with nonlinear conjugate gradients. The method is applied to a complex structure-controlled wadi area exhibiting faults, dissolution, collapse, and subsidence where the high resolution FWI velocity modeling helps clarifying the geological interpretation. The developed algorithms and automated workflows provide an effective solution for massive full waveform inversion of land seismic data that can be embedded in typical near surface velocity analysis procedures.


2020 ◽  
Vol 39 (5) ◽  
pp. 324-331
Author(s):  
Gary Murphy ◽  
Vanessa Brown ◽  
Denes Vigh

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.


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.


Geophysics ◽  
2016 ◽  
Vol 81 (4) ◽  
pp. S195-S206 ◽  
Author(s):  
Mrinal Sinha ◽  
Gerard T. Schuster

Imaging seismic data with an erroneous migration velocity can lead to defocused migration images. To mitigate this problem, we first choose a reference reflector whose topography is well-known from the well logs, for example. Reflections from this reference layer are correlated with the traces associated with reflections from deeper interfaces to get crosscorrelograms. Interferometric least-squares migration (ILSM) is then used to get the migration image that maximizes the crosscorrelation between the observed and the predicted crosscorrelograms. Deeper reference reflectors are used to image deeper parts of the subsurface with a greater accuracy. Results on synthetic and field data show that defocusing caused by velocity errors is largely suppressed by ILSM. We have also determined that ILSM can be used for 4D surveys in which environmental conditions and acquisition parameters are significantly different from one survey to the next. The limitations of ILSM are that it requires prior knowledge of a reference reflector in the subsurface and the velocity model below the reference reflector should be accurate.


2017 ◽  
Vol 5 (4) ◽  
pp. SR23-SR33 ◽  
Author(s):  
Xin Cheng ◽  
Kun Jiao ◽  
Dong Sun ◽  
Zhen Xu ◽  
Denes Vigh ◽  
...  

Over the past decade, acoustic full-waveform inversion (FWI) has become one of the standard methods in the industry to construct high-resolution velocity fields from the seismic data acquired. While most of the successful applications are for marine acquisition data with rich low-frequency diving or postcritical waves at large offsets, the application of acoustic FWI on land data remains a challenging topic. Land acoustic FWI application faces many severe difficulties, such as the presence of strong elastic effects, large near-surface velocity contrast, and heterogeneous, topography variations, etc. In addition, it is well-known that low-frequency transmitted seismic energy is crucial for the success of FWI to overcome sensitivity to starting velocity fields; unfortunately, those are the parts of the data that suffer the most from a low signal-to-noise ratio (S/N) in land acquisition. We have developed an acoustic FWI application on a land data set from North Kuwait, and demonstrated our solutions to mitigate some of the challenges posed by land data. More specifically, we have developed a semblance-based high-resolution Radon (HR-Radon) inversion approach to enhance the S/N of the low-frequency part of the FWI input data and to ultimately improve the convergence of the land FWI workflow. To mitigate the impact of elastic effects, we included only the diving and postcritical early arrivals in the waveform inversion. Our results show that, with the aid of HR-Radon preconditioning and a carefully designed workflow, acoustic FWI has the ability to derive a reliable high-resolution near-surface model that could not be otherwise recovered through traditional tomographic methods.


Geophysics ◽  
2018 ◽  
Vol 83 (6) ◽  
pp. Q49-Q69
Author(s):  
Jixiang Xu ◽  
Shitai Dong ◽  
Huajuan Cui ◽  
Yan Zhang ◽  
Ying Hu ◽  
...  

Near-surface scattered waves (NSWs) are the main noise in seismic data in areas with a complex near surface and can be divided into surface-to-surface scattered waves and body-to-surface scattered waves. We have developed a method for NSW enhancement that uses modified source-receiver interferometry. The method consists of two parts. First, deconvolutional intersource interferometry is used to cancel the common raypath of seismic waves from a near-surface scatterer to the common receiver and the receiver function. Second, convolutional interreceiver interferometry is used to compensate the common raypath of seismic waves from the common source to the near-surface scatterer and the source function. For an isotropic point scatterer near the earth’s surface in modified source-receiver interferometry, a body-to-surface scattered wave can be reconstructed by constructive interference not only among three body-to-surface scattered waves but also among a body-to-surface scattered wave and two surface-to-surface scattered waves; a surface-to-surface scattered wave can be reconstructed by constructive interference not only among three surface-to-surface scattered waves but also among a surface-to-surface scattered wave and two body-to-surface scattered waves. According to stationary phase analysis based on the superposition principle, we have developed a so-called dual-wheel driving configuration of modified source-receiver interferometry for enhancing NSWs in the data of conventional seismic exploration. The main advantages of the scheme are that (1) it can be used to enhance NSWs without the need for any a priori knowledge of topography and near-surface velocity, (2) it can be used to reconstruct NSWs from real sources to real receivers, including 3D near-surface side-scattered waves, and (3) it can be applied to conventional seismic data with finite-frequency bandwidth, spatially limited and sparse arrays, different source and receiver functions, and static correction. Numerically simulated data and field seismic data are used to demonstrate the feasibility and effectiveness of the scheme.


Geophysics ◽  
1997 ◽  
Vol 62 (5) ◽  
pp. 1558-1569 ◽  
Author(s):  
Dimitri Bevc

Wave‐equation datuming overcomes some of the problems that seismic data recorded on rugged surface topography present in routine image processing. The main problems are that (1) standard, optimized migration and processing algorithms assume data are recorded on a flat surface, and that (2) the static correction applied routinely to compensate for topography is inaccurate for waves that do not propagate vertically. Wave‐based processes such as stacking, dip‐moveout correction, normal‐moveout correction, velocity analysis, and migration after static shift can be severely affected by the nonhyperbolic character of the reflections. To alleviate these problems, I apply wave‐equation datuming early in the processing flow to upward continue the data to a flat datum, above the highest topography. This is what I refer to as “flooding the topography.” This approach does not require detailed a priori knowledge of the near‐surface velocity, and it streamlines subsequent processing because the data are regridded onto a regularly sampled datum. Wave‐equation datuming unravels the distortions caused by rugged topography, and unlike the static shift method, it does not adversely effect subsequent wave‐based processing. The image obtained after wave‐equation datuming exhibits better reflector continuity and more accurately represents the true structural image than the image obtained after static shift.


Sign in / Sign up

Export Citation Format

Share Document