Imaging Noisy Seismic Data using a One Dimensional Inverse Scattering Algorithm

2015 ◽  
Vol 12 (4) ◽  
Author(s):  
Bogdan Nita ◽  
Christopher Smith

We test the capability of an inverse scattering algorithm for imaging noisy seismic data. The algorithm does not require a velocity model or any other a priori information about the medium under investigation. We use three different geometries which capture different types of one-dimensional media with variable velocity. We show that the algorithm can precisely locate the interfaces and discover the correct velocity changes at those interfaces under moderate noise condition. When the signal to noise ratio is too small, the data is de-noised using a threshold filter and then imaged with excellent results. KEYWORDS: Seismic Imaging, Inversion, Amplitude Correction, Scattering Theory, Noise, Threshold Filter. 2000 MATHEMATICS SUBJECT CLASSIFICATION 86A22, 35J05, 35R30.

2017 ◽  
Vol 5 (3) ◽  
pp. SJ81-SJ90 ◽  
Author(s):  
Kainan Wang ◽  
Jesse Lomask ◽  
Felix Segovia

Well-log-to-seismic tying is a key step in many interpretation workflows for oil and gas exploration. Synthetic seismic traces from the wells are often manually tied to seismic data; this process can be very time consuming and, in some cases, inaccurate. Automatic methods, such as dynamic time warping (DTW), can match synthetic traces to seismic data. Although these methods are extremely fast, they tend to create interval velocities that are not geologically realistic. We have described the modification of DTW to create a blocked dynamic warping (BDW) method. BDW generates an automatic, optimal well tie that honors geologically consistent velocity constraints. Consequently, it results in updated velocities that are more realistic than other methods. BDW constrains the updated velocity to be constant or linearly variable inside each geologic layer. With an optimal correlation between synthetic seismograms and surface seismic data, this algorithm returns an automatically updated time-depth curve and an updated interval velocity model that still retains the original geologic velocity boundaries. In other words, the algorithm finds the optimal solution for tying the synthetic to the seismic data while restricting the interval velocity changes to coincide with the initial input blocking. We have determined the application of the BDW technique on a synthetic data example and field data set.


2011 ◽  
Vol 403-408 ◽  
pp. 2337-2340
Author(s):  
Shu Cong Liu ◽  
Yan Xing Song ◽  
Jing Song Yang

Seismic illumination analysis was an effective means of recognizing and studying the energy distributions in the underground geological structure in seismic data acquisition. Effective seismic illumination analysis to a priori targeted-geological model to identify the energy distribution of seismic waves, can apply to seismic analysis and amplitude compensation analysis. To increase the signal to noise ratio and resolution of seismic data when vibrator seismic exploration, it was necessary to strengthen the energy of a certain direction to get the High-Precision imaging and the best illumination of the target areas.Simulation research were done on single source directional illumination seismic technology, with seismic illumination analysis, and the impact of source number, spacing change on directional illumination seismic technology were also analyzed. Simulation results showed that the directional seismic technology could improved SNR of seismic data, and could be used for seismic signal processing.


Geophysics ◽  
1981 ◽  
Vol 46 (8) ◽  
pp. 1116-1120 ◽  
Author(s):  
A. B. Weglein ◽  
W. E. Boyse ◽  
J. E. Anderson

We present a formalism for obtaining the subsurface velocity configuration directly from reflection seismic data. Our approach is to apply the results obtained for inverse problems in quantum scattering theory to the reflection seismic problem. In particular, we extend the results of Moses (1956) for inverse quantum scattering and Razavy (1975) for the one‐dimensional (1-D) identification of the acoustic wave equation to the problem of identifying the velocity in the three‐dimensional (3-D) acoustic wave equation from boundary value measurements. No a priori knowledge of the subsurface velocity is assumed and all refraction, diffraction, and multiple reflection phenomena are taken into account. In addition, we explain how the idea of slant stack in processing seismic data is an important part of the proposed 3-D inverse scattering formalism.


2015 ◽  
Vol 3 (1) ◽  
pp. SB17-SB22 ◽  
Author(s):  
Richard C. Bain

Reliance on prestack time-migrated seismic data to define structural highs without incorporating all subsurface data and without taking into account the regional and local lateral depositional trends may result in dry holes or poorly positioned production wells due to local velocity changes, which are usually caused by some depositional or structural phenomenon. Tying check-shot control to depositional units may reveal those phenomena and permit assumptions to be made about velocities in areas beyond check-shot control points. We discovered a significant gas accumulation in an area surrounded by dry holes and marginal wells in the Vicksburg Formation in McAllen Ranch Field, Hidalgo County, Texas, by treating a seismic velocity anomaly as a geologic problem and by simple application of arithmetic and geometry to a 3D velocity model. Due to the effects of the anomaly, seismic data displayed in time gave no indication of the existence of a 325 ha (800 ac), 150 BCFG anticlinal structure. A subsurface model that accounted for the velocity anomaly was able to predict its extent and severity by readily identifiable thickness changes in the anomalous units. The resulting discovery yielded a sevenfold increase in field production within a two-year time span.


Geophysics ◽  
1982 ◽  
Vol 47 (10) ◽  
pp. 1422-1430 ◽  
Author(s):  
S. Raz

A new space‐time approach to inverting multidimensional seismic data is formulated. This work constitutes a natural, but assuredly nontrivial, extension of a previously reported one‐dimensional (1-D) Bremmer inversion procedure. In its most general format, the scheme applies to the three‐dimensional (3-D) inversion problem; however, appropriately reduced forms applicable to two‐dimensional (2-D) as well as 1-D configurations are also presented. The basic scheme and all its variants utilize exclusively zero‐offset data. Although interpretable in terms of a distorted wave Born model, the inversion procedure and subsequent algorithms differ substantially from their Born‐based predecessors. Here, the background is not only inhomogeneous but is not known a priori. The choice of a uniquely defined “phase‐corrected” background is cardinal to the proposed reconstruction scheme. The improvement in the inversion accuracy is considerable. Comparatively large discontinuities can be handled, and the phenomenon of error accumulation with depth is overcome.


Geophysics ◽  
2011 ◽  
Vol 76 (4) ◽  
pp. S151-S155 ◽  
Author(s):  
Mikhail Baykulov ◽  
Stefan Dümmong ◽  
Dirk Gajewski

A processing workflow was introduced for reflection seismic data that is based entirely on common-reflection-surface (CRS) stacking attributes. This workflow comprises the CRS stack, multiple attenuation, velocity model building, prestack data enhancement, trace interpolation, and data regularization. Like other methods, its limitation is the underlying hyperbolic assumption. The CRS workflow provides an alternative processing path in case conventional common midpoint (CMP) processing is unsatisfactory. Particularly for data with poor signal-to-noise ratio and low-fold acquisition, the CRS workflow is advantageous. The data regularization feature and the ability of prestack data enhancement provide quality control in velocity model building and improve prestack depth-migrated images.


Geophysics ◽  
2020 ◽  
pp. 1-45
Author(s):  
German Garabito ◽  
Paul L. Stoffa ◽  
Yuri S. F. Bezerra ◽  
João L. Caldeira

The application of the reverse time migration (RTM) in land seismic data is still a great challenge due to its low quality, low signal-to-noise ratio, irregular spatial sampling, acquisition gaps, missing traces, etc. Therefore, prior to the application of this kind of depth migration, the input pre-stack data must be conveniently preconditioned, that is, it must be interpolated, regularized, and enhanced. There are several methods for seismic data preconditioning, but for 2D real land data, the regularization of pre-stack data based on common reflection surface (CRS) stack method provides high quality enhanced preconditioned data, which is suitable for pre-stack depth migration and velocity model building. This work demonstrates the potential of RTM combined with CRS-based pre-stack data regularization, applied to real land seismic data with low quality and irregularly sparse spatial sampled, from geologically complex areas with the presence of diabase sills and steep dip reflections. Usually, determining the wavelet of the seismic source from land data is a challenge, because of this, RTM migration is often applied using artificial sources (e.g. Ricker). In this work, from the power spectrum of the pre-stacked data, we determine the wavelet of the seismic source to apply the RTM to real land data. We present applications of the pre-stack data preconditioning based on CRS stack and of the RTM in 2D land data of Tacutu and Parnaiba Basins, Brazil. Comparisons with the standard Kirchhoff depth migration reveals that the RTM improves the quality and resolution of the migrated images.


2007 ◽  
Author(s):  
Sverre Brandsberg-Dahl ◽  
Brian E. Hornby ◽  
Xiang Xiao

Author(s):  
Ehsan Jamali Hondori ◽  
Chen Guo ◽  
Hitoshi Mikada ◽  
Jin-Oh Park

AbstractFull-waveform inversion (FWI) of limited-offset marine seismic data is a challenging task due to the lack of refracted energy and diving waves from the shallow sediments, which are fundamentally required to update the long-wavelength background velocity model in a tomographic fashion. When these events are absent, a reliable initial velocity model is necessary to ensure that the observed and simulated waveforms kinematically fit within an error of less than half a wavelength to protect the FWI iterative local optimization scheme from cycle skipping. We use a migration-based velocity analysis (MVA) method, including a combination of the layer-stripping approach and iterations of Kirchhoff prestack depth migration (KPSDM), to build an accurate initial velocity model for the FWI application on 2D seismic data with a maximum offset of 5.8 km. The data are acquired in the Japan Trench subduction zone, and we focus on the area where the shallow sediments overlying a highly reflective basement on top of the Cretaceous erosional unconformity are severely faulted and deformed. Despite the limited offsets available in the seismic data, our carefully designed workflow for data preconditioning, initial model building, and waveform inversion provides a velocity model that could improve the depth images down to almost 3.5 km. We present several quality control measures to assess the reliability of the resulting FWI model, including ray path illuminations, sensitivity kernels, reverse time migration (RTM) images, and KPSDM common image gathers. A direct comparison between the FWI and MVA velocity profiles reveals a sharp boundary at the Cretaceous basement interface, a feature that could not be observed in the MVA velocity model. The normal faults caused by the basal erosion of the upper plate in the study area reach the seafloor with evident subsidence of the shallow strata, implying that the faults are active.


2021 ◽  
Vol 13 (1) ◽  
pp. 168781402098732
Author(s):  
Ayisha Nayyar ◽  
Ummul Baneen ◽  
Syed Abbas Zilqurnain Naqvi ◽  
Muhammad Ahsan

Localizing small damages often requires sensors be mounted in the proximity of damage to obtain high Signal-to-Noise Ratio in system frequency response to input excitation. The proximity requirement limits the applicability of existing schemes for low-severity damage detection as an estimate of damage location may not be known  a priori. In this work it is shown that spatial locality is not a fundamental impediment; multiple small damages can still be detected with high accuracy provided that the frequency range beyond the first five natural frequencies is utilized in the Frequency response functions (FRF) curvature method. The proposed method presented in this paper applies sensitivity analysis to systematically unearth frequency ranges capable of elevating damage index peak at correct damage locations. It is a baseline-free method that employs a smoothing polynomial to emulate reference curvatures for the undamaged structure. Numerical simulation of steel-beam shows that small multiple damages of severity as low as 5% can be reliably detected by including frequency range covering 5–10th natural frequencies. The efficacy of the scheme is also experimentally validated for the same beam. It is also found that a simple noise filtration scheme such as a Gaussian moving average filter can adequately remove false peaks from the damage index profile.


Sign in / Sign up

Export Citation Format

Share Document