Adaptive linear prediction filtering for random noise attenuation Mauricio D. Sacchi* and Mostafa Naghizadeh, University of Alberta
|
|
- Julian Cole
- 6 years ago
- Views:
Transcription
1 Adaptive linear prediction filtering for random noise attenuation Mauricio D. Sacchi* and Mostafa Naghizadeh, University of Alberta SUMMARY We propose an algorithm to compute time and space variant prediction filters for signal-to-noise ratio enhancement. Prediction filtering for seismic signal enhancement is, in general, implemented via filters that are estimated from the inversion of a system of equations in the t x or f x domain. In addition, prediction error filters are applied in small windows where the data can be modeled via a finite number of plane waves. Our algorithm, on the other hand, does not require the inversion of matrices. Furthermore, it does not require spatio-temporal windowing; the algorithm is implemented via a recursive scheme where the filter is continuously adapted to predict the signal. We postulate the prediction problem as a local smoothing problem and use a quadratic constraint to avoid solutions that model the noise. The algorithm uses a t x recursive implementation where the prediction filter for a given observation point is estimated via a simple rule. It turns out that the proposed algorithm is equivalent to the LMS (Least Mean Squares) filter often used for adaptive filtering. It is important to mention, however, that our derivation follows the framework that it is often used to solve underdetermined linear inverse problems. The latter involves the minimization of a cost function that includes a quadratic constraint to guarantee a stable solution. Synthetic and real data examples are used to test the algorithm. In particular, a field data test shows that adaptive t x filtering could offer an efficient and versatile alternative to classical f x deconvolution filtering. INTRODUCTION Prediction filters play an important role in seismic data processing with applications ranging from seismic deconvolution, signal-to-noise-ratio enhancement (Canales, 1984; Gulunay, 1986; Abma, 1995) and trace interpolation (Spitz, 1991). Prediction filters are often estimated from small spatio-temporal windows where waveforms can be approximated by events with constant dip. The latter is required for the optimal performance of the prediction filter. We propose to avoid windowing via a recursive algorithm where one prediction filter is estimated for each data sample. We show that the aforementioned problem is underdetermined and, as it is wellknown, admits an infinite number of solutions. A unique and stable solution if found by formulating the problem in terms of a regularization constraint. The filter required to smooth a given data point is constrained to be similar to the filter used to smooth an adjacent data point. Our formulation follows the classical approach used to solve an underdetermined problem. The final algorithm is equivalent to the LMS (Least Mean Squares) algorithm described in Widrow and Stearns (1985) and Hornbostel (1989). ADAPTIVE PREDICTION FILTERS 1-D case We start to formulate our problem by introducing an auto-regressive (AR) or linear prediction operator modeling operator of order p. To avoid notational clutter we first consider a 1D problem with p = 3 x k = a 1 x k 1 + a 2 x k 2 + a 3 x k 3 + ε k. (1) It is important to point out that one could have also proposed an analysis in terms of prediction error filters x k a 1 x k 1 a 2 x k 2 a 3 x k 3 = ε k, (2) where by calling f 0 = 1, and f k = a k, k = 1,2,3, equation (1) becomes x k f 0 + f 1 x k 1 + f 2 x k 2 + f 3 x k 3 = ε k. (3) We now assume that the coefficients of the autoregressive model are varying in time or space. In other words, each output sample k is predicted via coefficients that depend on k. The latter can be mathematically expressed as follows or, in matrix form x k = a k 1 x k 1 + a k 2 x k 2 + a k 3 x k 3 + ε k. (4) a k 1 x k = (x k 1 x k 2 x k 3 ) a k 2 a k 3 = M k a k. It is clear that for each sample k we have p = 3 unknown coefficients. This is an underdetermined problem that admits an infinite number of solutions. Uniqueness of the solution is guaranteed by introducing a quadratic regularization term. The local prediction filter for sample k is found by minimizing the following cost function (5) J k = x k M k a k µ 2 a k a (6) The first term in the right hand side of equation (6), the misfit function, measures how well the filter predicts the signal. The second term, the regularization term, measures the closeness of the local prediction filter to an a priori filter a 0. The regularization parameter µ determines the importance of the regularization term relative to the data misfit. Taking the derivatives with respect to the unknown filter coefficients and, after setting them to zero, leads to the following system of equations a k = (I (M H k M k + µ 2 I) 1 M H k M K)a 0 +(M H M + µ 2 I) 1 M H k x k. Where (.) H is used to indicate the conjugate transpose. Notice that our analysis considers the general case where the signal x k is complex. We now adopt the following identity to simplify equation (7) (7) (M H k M k + µ 2 I) 1 M H k = MH k (MMH + µ 2 I) 1 (8) 3347
2 Adaptive filtering and after considering that M k = (x n 1 x n 2 x n 3 ), δ 2 k = M km H k = 3 n=1 x k n 2. (9) After a few algebraic manipulations we arrive to the following expression â k = (I MH k M k δk 2 )a +µ2 0 + δk 2 x +µ2 k. (10) It is clear that the locally smoothed signal (equation (5)) is now given by ˆx k = M k â k = (1 δ k 2 δk 2 )M ka 0 + +µ2 MH k δ k 2 δk 2 x +µ2 k (11) The trade-off parameter µ controls the degree of adaptability to changes in the signal. In particular, we examine the influence of µ in the following two situations: 1 Minimum norm solution with perfect data fitting µ 0, ˆx k = x k. (12) The filter is predicting the original data and therefore, signal and noise are simultaneously modeled. This is often referred in inverse problems as over-fitting. 2 No adaptability µ, ˆx k = M k a 0 = x k 1 a 1,0 + x k 2 a 2,0 + x k 3 a 3,0. (13) The prediction filter cannot adapt to changes in x k and, therefore, smoothing becomes independent of k. Last equation is equivalent to conventional non-adaptive f x or t x prediction filtering (Gulunay, 1986; Abma, 1995) where the filter a 0 is estimated from the data via the least-squares method. An optimal value of µ needs to be selected in order to achieve noise suppression and adaptability to signal changes (e.g. changes in dip in the t x domain). We will return to this point in the section devoted to tests. Recursive algorithm We will assume that the signal varies slowly in time or space and, consequently, we propose The latter leads to the following algorithm Initialize a 0 and for k = 1,2,3... a 0 = a k 1. (14) a k = a k 1 η k M H k (M ka k 1 x k ) ˆx k = M k a k (15) where η k = 1 δ 2 k +µ2. 2D (t x) algorithm The 1-D algorithm can be easily extended to process 2-D and 3-D data. For that purpose, we first replace the 1-D inner product M H k a k = j x k j a k j by its multidimensional counterpart which, for a 2D prediction filter, is given by i x k i,l j a k,l i, j. j With this rule in mind, it is easy to generalize (15) to 2-D and 3-D algorithms. For instance, the 2-D t x adaptive algorithm can be implemented via the following expressions δ 2 k,l = p n=1 q m=1 x k n,l m 2 η k,l = 1 δ 2 k,l +µ2 a k,l i, j = a k 1,l 1 i, j η k,l x k i,l j ( p i=1 q j=1 x k i,l j a k,l i, j x k,l) ˆx k,l = p i=1 q j=1 x k i,l ja k,l i, j where we have considered prediction filters of size p q. Tests Synthetic example Figure 1 portrays a synthetic data composed of 5 events with parabolic moveout. The data were contaminated with Gaussian band-limited noise with signal-to-noise ratio SNR = 1. Figure 1a portrays the noisy data. Figure 1b portrays the filtered data obtained via the algorithm described in this paper. We have adopted a parameter µ = 20 and 2D (t x) adaptive prediction filters of size Figure 1c displays the noise section (data minus filtered data). The tradeoff parameter is obtained by inspecting the noise section. The optimal µ yields a noise section with little amount of coherent energy. In this example the adaptive filter was run in space-time starting from the first observation (t = 0,x = 0) and moving toward the last observation (t = t max,x = x max ) with initial prediction filter a 0,0 k,l = 0. The algorithm is rerun starting from the last observation and moving the recursion toward the first observation. The initial filter in this case is also given by a tmax,xmax k,l = 0. By this procedure, two smoothed signals were estimated. The average signal was adopted as the final estimate of the filtered data. We have also tried to initialize the algorithm with a least-squares t x filter computed in a small window (Abma and Claerbout, 1995). The aforementioned initialization did not provide important improvements. Field data example We continue with our tests by filtering a real data set from the Western Canadian Basin. The data consist of a seismic section with an appreciable deterioration in SNR for the deep part of the image. Figure 2 shows the mean normalized spectra of the data after and before filtering. Figure 3a shows the section prior to filtering. Figure 3b provides the filtered section via f x deconvolution. The f x filter was computed using the least-squares solution to the forward/backward prediction filtering problem. The method avoids finite aperture artifacts by considering non-toeplitz covariance matrices ( Marple, 1987). The length of the f x prediction filter is L = 4, the 3348
3 Adaptive filtering 0 (a) (b) (c) time (s) Figure 1: a) Synthetic seismic gather contaminated with Gaussian band-limited noise (SNR = 1). b) Output of the t x adaptive prediction filtering algorithm. c) Difference between the data and the filtered data. covariance matrix was inverted with a tradeoff parameter that represents 0.1% of its trace. Figure 3c present the filtered data after applying the t x adaptive prediction filter with an operator of size (time space). The tradeoff parameter for this example is set to µ = 15. The tradeoff parameter was found by examining the residual section. In other words, the algorithm is run with different tradeoff parameters; the optimal tradeoff is the one that leads to a residual section with a minimum amount of coherent energy. Power Input f x decon t x adaptive Frequency (Hz) strategies to cope with spatial changes in dip. As a final remark, it is interesting to point out that the adaptive method can be easily used to design 2D algorithms that operate in t x or in f x y for f slices. It can also be extended to the 3D case to cope with noise attenuation in t x y. In summary, adaptive filtering offers a versatile design that can be easily accommodated to noise attenuation in different domains. CONCLUSIONS We have presented an inversion-based approach to the design of adaptive filters. Our approach involves considering adaptive filtering as an undetermined linear inverse problem where a set of filter coefficients are used to predict one observation point at the time. The problem, clearly, involves estimating more than one parameter from a single observation. A quadratic regularization constraint was used to retrieve a filter that can locally predict the signal and, at the same time, adapt to temporal and spatial changes in the data. When the a priori filter (initial filter) is the filter compute from a nearby observation point, the method leads to a recursive algorithm that is equivalent to the LMS filter. Synthetic and real data tests were used to validate the adaptive algorithm as a tool for noise suppression. In particular, our field data tests show that adaptive filtering is a viable alternative to classical f x filtering for noise suppression. ACKNOWLEDGMENTS This research has been supported by the sponsors of the Signal Analysis and Imaging Group, at the University of Alberta and the National Sciences and Engineering Research Council of Canada via a Discovery Grant to MDS. Discussions with Sam Kaplan and Tad Ulrych were instrumental in the preparation of this paper. Figure 2: Power spectra of the data and filtered data in Figure 3. The yellow boxes in Figure3a,b, and c were redisplayed in Figures 4 a, b and c, respectively. In addition, Figures 4d and e show the estimators of the noise. The field data example clearly shows that both method have lead to similar results. The main advantage of the adaptive procedure roots in its simplicity. In addition, one does not have to worry about inverting large matrices (Abma and Clearbout, 1995) or to specify windowing 3349
4 Adaptive filtering Figure 3: a) Seismic section from a survey in the Western Canadian Basin. b) Seismic section after f x deconvolution. c) Seismic section after noise attenuation via adaptive t x prediction filtering. The boxes are redisplayed in Figure 4. (a) (b) (c) (d) (e) Figure 4: a) Original data. b) Data after signal-to-noise ratio enhancement via f x deconvolution. c) Data after signal-to-noise ratio enhancement via adaptive t x prediction filtering. d) Noise removed from the data via f x deconvolution. e) Noise removed from the data via adaptive t x prediction filtering. The data correspond to every second trace in the boxes in Figure
5 EDITED REFERENCES Note: This reference list is a copy-edited version of the reference list submitted by the author. Reference lists for the 2009 SEG Technical Program Expanded Abstracts have been copy edited so that references provided with the online metadata for each paper will achieve a high degree of linking to cited sources that appear on the Web. REFERENCES Abma, R., and J. Claerbout, 1995, Lateral prediction for noise attenuation by t-x and F-X techniques: Geophysics, 60, no. 6, Canales, L. L., 1984, Random noise reduction: 54th Annual International Meeting, SEG, Expanded Abstracts, Gulunay, N., 1986, Fx decon and complex Wiener prediction filter: 56th Annual International Meeting, SEG, Expanded Abstracts, Hornbostel S.,1989, Spatial prediction filtering in the t-x and f-x domains: Geophysics, 63, Marple, Lawrence S., 1997, Digital Spectral Analysis: With Applications, Prentice Hall. Spitz, S., 1991, Seismic trace interpolation in the F-X domain: Geophysics, 56, no. 6, Widrow, B., and S. D. Stearns, 1985 Adaptive Signal Processing, Prentice Hall. 3351
arxiv: v1 [physics.geo-ph] 23 Dec 2017
Statics Preserving Sparse Radon Transform Nasser Kazemi, Department of Physics, University of Alberta, kazemino@ualberta.ca Summary arxiv:7.087v [physics.geo-ph] 3 Dec 07 This paper develops a Statics
More informationStatics preserving projection filtering Yann Traonmilin*and Necati Gulunay, CGGVeritas
Yann Traonmilin*and Necati Gulunay, CGGVeritas Summary Projection filtering has been used for many years in seismic processing as a tool to extract a meaningful signal out of noisy data. We show that its
More informationInterpolation of nonstationary seismic records using a fast non-redundant S-transform
Interpolation of nonstationary seismic records using a fast non-redundant S-transform Mostafa Naghizadeh and Kris A. Innanen CREWES University of Calgary SEG annual meeting San Antonio 22 September 2011
More informationAdaptive multiple subtraction using regularized nonstationary regression
GEOPHSICS, VOL. 74, NO. 1 JANUAR-FEBRUAR 29 ; P. V25 V33, 17 FIGS. 1.119/1.343447 Adaptive multiple subtraction using regularized nonstationary regression Sergey Fomel 1 ABSTRACT Stationary regression
More informationMinimum entropy deconvolution with frequency-domain constraints
GEOPHYSICS, VOL. 59, NO. 6 (JUNE 1994); P. 938-945, 9 FIGS., 1 TABLE. Minimum entropy deconvolution with frequency-domain constraints Mauricio D. Sacchi*, Danilo R. Velis*, and Alberto H. Cominguez ABSTRACT
More informationElastic full waveform inversion for near surface imaging in CMP domain Zhiyang Liu*, Jie Zhang, University of Science and Technology of China (USTC)
Elastic full waveform inversion for near surface imaging in CMP domain Zhiyang Liu*, Jie Zhang, University of Science and Technology of China (USTC) Summary We develop an elastic full waveform inversion
More informationRegularizing seismic inverse problems by model reparameterization using plane-wave construction
GEOPHYSICS, VOL. 71, NO. 5 SEPTEMBER-OCTOBER 2006 ; P. A43 A47, 6 FIGS. 10.1190/1.2335609 Regularizing seismic inverse problems by model reparameterization using plane-wave construction Sergey Fomel 1
More informationStanford Exploration Project, Report 111, June 9, 2002, pages INTRODUCTION FROM THE FILTERING TO THE SUBTRACTION OF NOISE
Stanford Exploration Project, Report 111, June 9, 2002, pages 199 204 Short Note Theoretical aspects of noise attenuation Antoine Guitton 1 INTRODUCTION In Guitton (2001) I presented an efficient algorithm
More information( ) ( ), 1, 1. Downloaded 10/14/14 to Redistribution subject to SEG license or copyright; see Terms of Use at
Chao Zhang*, Jie Zhang, University of Science and Technology of China (USTC); Zhikun Sun, GeoTomo LLC Summary Refraction traveltimes have long been applied for deriving long-wavelength statics solutions.
More informationOn non-stationary convolution and inverse convolution
Stanford Exploration Project, Report 102, October 25, 1999, pages 1 137 On non-stationary convolution and inverse convolution James Rickett 1 keywords: helix, linear filtering, non-stationary deconvolution
More informationRobust one-step (deconvolution + integration) seismic inversion in the frequency domain Ivan Priezzhev* and Aaron Scollard, Schlumberger
Robust one-step (deconvolution + integration) seismic inversion in the frequency domain Ivan Priezzhev and Aaron Scollard, Schlumberger Summary Seismic inversion requires two main operations relative to
More informationPressure Normal Derivative Extraction for Arbitrarly Shaped Surfaces Endrias G. Asgedom,Okwudili Chuks Orji, Walter Söllner, PGS
Pressure Normal Derivative Extraction for Arbitrarly Shaped Surfaces Endrias G. Asgedom,Okwudili Chuks Orji, Walter Söllner, PGS Downloaded 9// to... Redistribution subject to SEG license or copyright;
More informationAdvanced Digital Signal Processing -Introduction
Advanced Digital Signal Processing -Introduction LECTURE-2 1 AP9211- ADVANCED DIGITAL SIGNAL PROCESSING UNIT I DISCRETE RANDOM SIGNAL PROCESSING Discrete Random Processes- Ensemble Averages, Stationary
More informationCompensating visco-acoustic effects in anisotropic resverse-time migration Sang Suh, Kwangjin Yoon, James Cai, and Bin Wang, TGS
Compensating visco-acoustic effects in anisotropic resverse-time migration Sang Suh, Kwangjin Yoon, James Cai, and Bin Wang, TGS SUMMARY Anelastic properties of the earth cause frequency dependent energy
More informationResidual Statics using CSP gathers
Residual Statics using CSP gathers Xinxiang Li and John C. Bancroft ABSTRACT All the conventional methods for residual statics analysis require normal moveout (NMO) correction applied on the seismic data.
More informationSeismic attributes of time-vs. depth-migrated data using self-adaptive window
Seismic attributes of time-vs. depth-migrated data using self-adaptive window Tengfei Lin*, Bo Zhang, The University of Oklahoma; Zhifa Zhan, Zhonghong Wan, BGP., CNPC; Fangyu Li, Huailai Zhou and Kurt
More informationTrace balancing with PEF plane annihilators
Stanford Exploration Project, Report 84, May 9, 2001, pages 1 333 Short Note Trace balancing with PEF plane annihilators Sean Crawley 1 INTRODUCTION Land seismic data often shows large variations in energy
More informationLayer thickness estimation from the frequency spectrum of seismic reflection data
from the frequency spectrum of seismic reflection data Arnold Oyem* and John Castagna, University of Houston Summary We compare the spectra of Short Time Window Fourier Transform (STFT) and Constrained
More informationADAPTIVE FILTER THEORY
ADAPTIVE FILTER THEORY Fourth Edition Simon Haykin Communications Research Laboratory McMaster University Hamilton, Ontario, Canada Front ice Hall PRENTICE HALL Upper Saddle River, New Jersey 07458 Preface
More informationDownloaded 05/01/17 to Redistribution subject to SEG license or copyright; see Terms of Use at
Mapping Imbricate Structures in the Thrust Belt of Southeast Turkey by Large-Offset Seismic Survey Oz Yilmaz*, Anatolian Geophysical, Istanbul, Turkey; and Serdar Uygun, Ali Ölmez, and Emel Çalı, Turkish
More informationMultiple realizations using standard inversion techniques a
Multiple realizations using standard inversion techniques a a Published in SEP report, 105, 67-78, (2000) Robert G Clapp 1 INTRODUCTION When solving a missing data problem, geophysicists and geostatisticians
More informationThe i-stats: An Image-Based Effective-Medium Modeling of Near-Surface Anomalies Oz Yilmaz*, GeoTomo LLC, Houston, TX
The i-stats: An Image-Based Effective-Medium Modeling of Near-Surface Anomalies Oz Yilmaz*, GeoTomo LLC, Houston, TX Summary Near-surface modeling for statics corrections is an integral part of a land
More informationSummary. Introduction
Detailed velocity model building in a carbonate karst zone and improving sub-karst images in the Gulf of Mexico Jun Cai*, Hao Xun, Li Li, Yang He, Zhiming Li, Shuqian Dong, Manhong Guo and Bin Wang, TGS
More informationStanford Exploration Project, Report 105, September 5, 2000, pages 41 53
Stanford Exploration Project, Report 105, September 5, 2000, pages 41 53 40 Stanford Exploration Project, Report 105, September 5, 2000, pages 41 53 Short Note Multiple realizations using standard inversion
More informationEvidence of an axial magma chamber beneath the ultraslow spreading Southwest Indian Ridge
GSA Data Repository 176 1 5 6 7 9 1 11 1 SUPPLEMENTARY MATERIAL FOR: Evidence of an axial magma chamber beneath the ultraslow spreading Southwest Indian Ridge Hanchao Jian 1,, Satish C. Singh *, Yongshun
More informationFrequency-Domain Rank Reduction in Seismic Processing An Overview
Frequency-Domain Rank Reduction in Seismic Processing An Overview Stewart Trickett Absolute Imaging Inc. Summary Over the last fifteen years, a new family of matrix rank reduction methods has been developed
More informationHigher-order statistics: improved event resolution?
Higher-order statistics: improved event resolution? David C. Henley and Arnim B. Haase Higher-order statistics ABSTRACT In any area where near-surface earth structure causes one or more small time delays
More informationAnisotropic tomography for TTI and VTI media Yang He* and Jun Cai, TGS
media Yang He* and Jun Cai, TGS Summary A simultaneous anisotropic tomographic inversion algorithm is developed. Check shot constraints and appropriate algorithm preconditioning play an important role
More informationSeismic data interpolation and denoising using SVD-free low-rank matrix factorization
Seismic data interpolation and denoising using SVD-free low-rank matrix factorization R. Kumar, A.Y. Aravkin,, H. Mansour,, B. Recht and F.J. Herrmann Dept. of Earth and Ocean sciences, University of British
More informationCalculation of the sun s acoustic impulse response by multidimensional
Calculation of the sun s acoustic impulse response by multidimensional spectral factorization J. E. Rickett and J. F. Claerbout Geophysics Department, Stanford University, Stanford, CA 94305, USA Abstract.
More informationDirect nonlinear traveltime inversion in layered VTI media Paul J. Fowler*, Alexander Jackson, Joseph Gaffney, and David Boreham, WesternGeco
Paul J. Fowler*, Alexander Jackson, Joseph Gaffney, and David Boreham, WesternGeco Summary We present a scheme for direct nonlinear inversion of picked moveout traveltimes in block-layered isotropic or
More information3D VTI traveltime tomography for near-surface imaging Lina Zhang*, Jie Zhang, Wei Zhang, University of Science and Technology of China (USTC)
Downloaded 01/03/14 to 16.01.198.34. Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/ 3D VTI traveltime tomography for near-surface imaging Lina Zhang*, Jie
More informationTOM 2.6. SEG/Houston 2005 Annual Meeting 2581
Oz Yilmaz* and Jie Zhang, GeoTomo LLC, Houston, Texas; and Yan Shixin, PetroChina, Beijing, China Summary PetroChina conducted a multichannel large-offset 2-D seismic survey in the Yumen Oil Field, Northwest
More informationImproving subsurface imaging in geological complex area: Structure PReserving INTerpolation in 6D (SPRINT6D)
Improving subsurface imaging in geological complex area: Structure PReserving INTerpolation in 6D (SPRINT6D) Dan Negut* and Mark Ng, Divestco Inc. Copyright 2008, ACGGP This paper was selected for presentation
More informationComparison between least-squares reverse time migration and full-waveform inversion
Comparison between least-squares reverse time migration and full-waveform inversion Lei Yang, Daniel O. Trad and Wenyong Pan Summary The inverse problem in exploration geophysics usually consists of two
More informationV003 How Reliable Is Statistical Wavelet Estimation?
V003 How Reliable Is Statistical Wavelet Estimation? J.A. Edgar* (BG Group) & M. van der Baan (University of Alberta) SUMMARY Well logs are often used for the estimation of seismic wavelets. The phase
More informationImplicit 3-D depth migration by wavefield extrapolation with helical boundary conditions
Stanford Exploration Project, Report 97, July 8, 1998, pages 1 13 Implicit 3-D depth migration by wavefield extrapolation with helical boundary conditions James Rickett, Jon Claerbout, and Sergey Fomel
More informationEstimating a pseudounitary operator for velocity-stack inversion
Stanford Exploration Project, Report 82, May 11, 2001, pages 1 77 Estimating a pseudounitary operator for velocity-stack inversion David E. Lumley 1 ABSTRACT I estimate a pseudounitary operator for enhancing
More informationRough sea estimation for phase-shift de-ghosting Sergio Grion*, Rob Telling and Seb Holland, Dolphin Geophysical
Rough sea estimation for phase-shift de-ghosting Sergio Grion*, Rob Telling and Seb Holland, Dolphin Geophysical Summary This paper discusses rough-sea de-ghosting for variabledepth streamer data. The
More informationAdaptive Filtering. Squares. Alexander D. Poularikas. Fundamentals of. Least Mean. with MATLABR. University of Alabama, Huntsville, AL.
Adaptive Filtering Fundamentals of Least Mean Squares with MATLABR Alexander D. Poularikas University of Alabama, Huntsville, AL CRC Press Taylor & Francis Croup Boca Raton London New York CRC Press is
More informationSeismic attributes for fault/fracture characterization
Satinder Chopra* and Kurt J. Marfurt + *Arcis Corporation, Calgary; + University of Houston, Houston Summary Seismic attributes have proliferated in the last three decades at a rapid rate and have helped
More informationImproved Radon transforms for filtering of coherent noise
Radon transforms for coherent-noise filtering Improved Radon transforms for filtering of coherent noise Shauna K. Oppert and R. James Brown ABSTRACT Radon transforms rely on the ability to predict the
More informationElastic least-squares migration with two-way wave equation forward and adjoint operators
Elastic least-squares migration with two-way wave equation forward and adjoint operators Ke Chen and Mauricio D. Sacchi, Department of Physics, University of Alberta Summary Time domain elastic least-squares
More information3D INTERPOLATION USING HANKEL TENSOR COMPLETION BY ORTHOGONAL MATCHING PURSUIT A. Adamo, P. Mazzucchelli Aresys, Milano, Italy
3D INTERPOLATION USING HANKEL TENSOR COMPLETION BY ORTHOGONAL MATCHING PURSUIT A. Adamo, P. Mazzucchelli Aresys, Milano, Italy Introduction. Seismic data are often sparsely or irregularly sampled along
More informationImproving Resolution with Spectral Balancing- A Case study
P-299 Improving Resolution with Spectral Balancing- A Case study M Fatima, Lavendra Kumar, RK Bhattacharjee, PH Rao, DP Sinha Western Offshore Basin, ONGC, Panvel, Mumbai Summary: The resolution limit
More informationEffect of Noise in Blending and Deblending Guus Berkhout* and Gerrit Blacquière, Delft University of Technology
Effect of Noise in Blending and De Guus Berkhout* and Gerrit Blacquière, Delft University of Technology Downloaded //3 to 3.8.3.78. Redistribution subject to SEG license or copyright; see Terms of Use
More informationSummary. We present the results of the near-surface characterization for a 3D survey in thrust belt area in Sharjah, United Arab Emirates.
Near-surface characterization, challenges, and solutions for high-density, high-productivity, Alexander Zarkhidze*, Claudio Strobbia, Abdallah Ibrahim, WesternGeco; Luis Viertel Herrera, Abdulla Al Qadi,
More informationAttenuation compensation in least-squares reverse time migration using the visco-acoustic wave equation
Attenuation compensation in least-squares reverse time migration using the visco-acoustic wave equation Gaurav Dutta, Kai Lu, Xin Wang and Gerard T. Schuster, King Abdullah University of Science and Technology
More informationAdaptive subtraction of free surface multiples through order-by-order prediction, matching filters and independent component analysis
Adaptive subtraction of free surface multiples through order-by-order prediction, matching filters and independent component analysis Sam T. Kaplan and Kristopher A. Innanen M-OSRP Annual Meeting; May
More informationSeismic processing of numerical EM data John W. Neese* and Leon Thomsen, University of Houston
Seismic processing of numerical EM data John W. Neese* and Leon Thomsen, University of Houston Summary The traditional methods for acquiring and processing CSEM data are very different from those for seismic
More information2010 SEG SEG Denver 2010 Annual Meeting
Anisotropic model building with wells and horizons: Gulf of Mexico case study comparing different approaches Andrey Bakulin*, Olga Zdraveva, Yangjun (Kevin) Liu, Kevin Lyons, WesternGeco/Schlumberger Summary
More informationPART 1. Review of DSP. f (t)e iωt dt. F(ω) = f (t) = 1 2π. F(ω)e iωt dω. f (t) F (ω) The Fourier Transform. Fourier Transform.
PART 1 Review of DSP Mauricio Sacchi University of Alberta, Edmonton, AB, Canada The Fourier Transform F() = f (t) = 1 2π f (t)e it dt F()e it d Fourier Transform Inverse Transform f (t) F () Part 1 Review
More informationStepwise Multiple Elimination using Linear Transforms: An alternative approach to SRME for stacking multiple-free near offsets in the Barents Sea
: An alternative approach to SRME for stacking multiple-free near offsets in the Barents Sea Anthony Hardwick*, Carl Lang and Bent Kjølhamar, TGS. Summary A novel, de-aliased method called Stepwise Multiple
More informationComplex-beam Migration and Land Depth Tianfei Zhu CGGVeritas, Calgary, Alberta, Canada
Page 1 of 10 Home Articles Interviews Print Editions Complex-beam Migration and Land Depth Tianfei Zhu CGGVeritas, Calgary, Alberta, Canada DECEMBER 2012 FOCUS ARTICLE Summary Gaussian-beam depth migration
More informationThe Deconvolution of Multicomponent Trace Vectors
The Deconvolution of Multicomponent Trace Vectors Xinxiang Li, Peter Cary and Rodney Couzens Sensor Geophysical Ltd., Calgary, Canada xinxiang_li@sensorgeo.com Summary Deconvolution of the horizontal components
More informationTOM 1.7. Sparse Norm Reflection Tomography for Handling Velocity Ambiguities
SEG/Houston 2005 Annual Meeting 2554 Yonadav Sudman, Paradigm and Dan Kosloff, Tel-Aviv University and Paradigm Summary Reflection seismology with the normal range of offsets encountered in seismic surveys
More informationARMA SPECTRAL ESTIMATION BY AN ADAPTIVE IIR FILTER. by JIANDE CHEN, JOOS VANDEWALLE, and BART DE MOOR4
560 R. BRU AND J. VITbRIA REFERENCES Y. H. Au-Yeung and Y. T. Poon, 3X3 orthostochastic matrices and the convexity of generalized numerical ranges, Linear Algebra Appl. 27:69-79 (1979). N. Bebiano, Some
More informationSEG/New Orleans 2006 Annual Meeting
Carmen C. Dumitrescu, Sensor Geophysical Ltd., and Fred Mayer*, Devon Canada Corporation Summary This paper provides a case study of a 3D seismic survey in the Leland area of the Deep Basin of Alberta,
More information3-D seismic continuity attribute for mapping discontinuities in a target zone: optimum parameters and evaluation
JOURNAL OF THE BALKAN GEOPHYSICAL SOCIETY, Vol. 2, No 4, November 1999, p. 112-119, 7 figs. 3-D seismic continuity attribute for mapping discontinuities in a target zone: optimum parameters and evaluation
More informationSeismic wavepropagation concepts applied to the interpretation of marine controlled-source electromagnetics
Seismic wavepropagation concepts applied to the interpretation of marine controlled-source electromagnetics Rune Mittet, EMGS SUMMARY are very useful for the interpretation of seismic data. Moreover, these
More informationDownloaded 07/03/14 to Redistribution subject to SEG license or copyright; see Terms of Use at
Downloaded 07/03/14 to 129.237.143.21. Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/ Rayleigh-wave diffractions due to a void in the layered half space
More informationGeophysical Data Analysis: Discrete Inverse Theory
Geophysical Data Analysis: Discrete Inverse Theory MATLAB Edition William Menke Lamont-Doherty Earth Observatory and Department of Earth and Environmental Sciences Columbia University. ' - Palisades, New
More informationAchieving depth resolution with gradient array survey data through transient electromagnetic inversion
Achieving depth resolution with gradient array survey data through transient electromagnetic inversion Downloaded /1/17 to 128.189.118.. Redistribution subject to SEG license or copyright; see Terms of
More informationRecursive Least Squares for an Entropy Regularized MSE Cost Function
Recursive Least Squares for an Entropy Regularized MSE Cost Function Deniz Erdogmus, Yadunandana N. Rao, Jose C. Principe Oscar Fontenla-Romero, Amparo Alonso-Betanzos Electrical Eng. Dept., University
More informationShort Note. Plane wave prediction in 3-D. Sergey Fomel 1 INTRODUCTION
Stanford Exploration Project, Report SERGEY, November 9, 2000, pages 291?? Short Note Plane wave prediction in 3-D Sergey Fomel 1 INTRODUCTION The theory of plane-wave prediction in three dimensions is
More informationEffect of velocity uncertainty on amplitude information
Stanford Exploration Project, Report 111, June 9, 2002, pages 253 267 Short Note Effect of velocity uncertainty on amplitude information Robert G. Clapp 1 INTRODUCTION Risk assessment is a key component
More informationSource estimation for frequency-domain FWI with robust penalties
Source estimation for frequency-domain FWI with robust penalties Aleksandr Y. Aravkin, Tristan van Leeuwen, Henri Calandra, and Felix J. Herrmann Dept. of Earth and Ocean sciences University of British
More informationMagnetic Case Study: Raglan Mine Laura Davis May 24, 2006
Magnetic Case Study: Raglan Mine Laura Davis May 24, 2006 Research Objectives The objective of this study was to test the tools available in EMIGMA (PetRos Eikon) for their utility in analyzing magnetic
More informationSummary. 2D potential field migration of gravity fields and their gradients. For a 2D gravity field, we can define the complex intensity:
Michael S. Zhdanov*, The University of Utah and TechnoImaging, Xiaojun Liu, The University of Utah, and Glenn A. Wilson, TechnoImaging Summary We introduce 2D potential field migration for the interpretation
More informationVollständige Inversion seismischer Wellenfelder - Erderkundung im oberflächennahen Bereich
Seminar des FA Ultraschallprüfung Vortrag 1 More info about this article: http://www.ndt.net/?id=20944 Vollständige Inversion seismischer Wellenfelder - Erderkundung im oberflächennahen Bereich Thomas
More informationDownloaded 07/03/14 to Redistribution subject to SEG license or copyright; see Terms of Use at
Applications of the JARS Method to study levee sites in southern Texas and southern New Mexico Julia Ivanov*, Richard D. Miller, Jianghai Xia, Kansas Geological Survey, Lawrence, KS Joseph B. Dunbar, Engineer
More informationEstimating vertical and horizontal resistivity of the overburden and the reservoir for the Alvheim Boa field. Folke Engelmark* and Johan Mattsson, PGS
Estimating vertical and horizontal resistivity of the overburden and the reservoir for the Alvheim Boa field. Folke Engelmark* and Johan Mattsson, PGS Summary Towed streamer EM data was acquired in October
More informationNOISE ROBUST RELATIVE TRANSFER FUNCTION ESTIMATION. M. Schwab, P. Noll, and T. Sikora. Technical University Berlin, Germany Communication System Group
NOISE ROBUST RELATIVE TRANSFER FUNCTION ESTIMATION M. Schwab, P. Noll, and T. Sikora Technical University Berlin, Germany Communication System Group Einsteinufer 17, 1557 Berlin (Germany) {schwab noll
More informationStructured tensor missing-trace interpolation in the Hierarchical Tucker format Curt Da Silva and Felix J. Herrmann Sept. 26, 2013
Structured tensor missing-trace interpolation in the Hierarchical Tucker format Curt Da Silva and Felix J. Herrmann Sept. 6, 13 SLIM University of British Columbia Motivation 3D seismic experiments - 5D
More informationAn Adaptive Sensor Array Using an Affine Combination of Two Filters
An Adaptive Sensor Array Using an Affine Combination of Two Filters Tõnu Trump Tallinn University of Technology Department of Radio and Telecommunication Engineering Ehitajate tee 5, 19086 Tallinn Estonia
More informationInterval anisotropic parameters estimation in a least squares sense Case histories from West Africa
P-263 Summary Interval anisotropic parameters estimation in a least squares sense Patrizia Cibin*, Maurizio Ferla Eni E&P Division (Milano, Italy), Emmanuel Spadavecchia - Politecnico di Milano (Milano,
More informationSeismic modeling evaluation of fault illumination in the Woodford Shale Sumit Verma*, Onur Mutlu, Kurt J. Marfurt, The University of Oklahoma
Seismic modeling evaluation of fault illumination in the Woodford Shale Sumit Verma*, Onur Mutlu, Kurt J. Marfurt, The University of Oklahoma Summary The Woodford Shale is one of the more important resource
More informationChurning seismic attributes with principal component analysis
Satinder Chopra + * and Kurt J. Marfurt + Arcis Seismic Solutions, Calgary; The University of Oklahoma, Norman Summary Seismic attributes are an invaluable aid in the interpretation of seismic data. Different
More information23855 Rock Physics Constraints on Seismic Inversion
23855 Rock Physics Constraints on Seismic Inversion M. Sams* (Ikon Science Ltd) & D. Saussus (Ikon Science) SUMMARY Seismic data are bandlimited, offset limited and noisy. Consequently interpretation of
More informationMain Menu. Douglas Oldenburg University of British Columbia Vancouver, BC, Canada
Using ERA low frequency E-field profiling and UBC 3D frequency-domain inversion to delineate and discover a mineralized zone in Porcupine district, Ontario, Canada. Vlad Kaminski* University of British
More informationLow-rank Promoting Transformations and Tensor Interpolation - Applications to Seismic Data Denoising
Low-rank Promoting Transformations and Tensor Interpolation - Applications to Seismic Data Denoising Curt Da Silva and Felix J. Herrmann 2 Dept. of Mathematics 2 Dept. of Earth and Ocean Sciences, University
More informationEE 367 / CS 448I Computational Imaging and Display Notes: Image Deconvolution (lecture 6)
EE 367 / CS 448I Computational Imaging and Display Notes: Image Deconvolution (lecture 6) Gordon Wetzstein gordon.wetzstein@stanford.edu This document serves as a supplement to the material discussed in
More informationA New Approach to Tune the Vold-Kalman Estimator for Order Tracking
A New Approach to Tune the Vold-Kalman Estimator for Order Tracking Amadou Assoumane, Julien Roussel, Edgard Sekko and Cécile Capdessus Abstract In the purpose to diagnose rotating machines using vibration
More informationCOMPLEX TRACE ANALYSIS OF SEISMIC SIGNAL BY HILBERT TRANSFORM
COMPLEX TRACE ANALYSIS OF SEISMIC SIGNAL BY HILBERT TRANSFORM Sunjay, Exploration Geophysics,BHU, Varanasi-221005,INDIA Sunjay.sunjay@gmail.com ABSTRACT Non-Stationary statistical Geophysical Seismic Signal
More informationWe A10 12 Common Reflection Angle Migration Revealing the Complex Deformation Structure beneath Forearc Basin in the Nankai Trough
We A10 12 Common Reflection Angle Migration Revealing the Complex Deformation Structure beneath Forearc Basin in the Nankai Trough K. Shiraishi* (JAMSTEC), M. Robb (Emerson Paradigm), K. Hosgood (Emerson
More informationADAPTIVE FILTER THEORY
ADAPTIVE FILTER THEORY Fifth Edition Simon Haykin Communications Research Laboratory McMaster University Hamilton, Ontario, Canada International Edition contributions by Telagarapu Prabhakar Department
More informationThin-Bed Reflectivity An Aid to Seismic Interpretation
Thin-Bed Reflectivity An Aid to Seismic Interpretation Satinder Chopra* Arcis Corporation, Calgary, AB schopra@arcis.com John Castagna University of Houston, Houston, TX, United States and Yong Xu Arcis
More informationIn the derivation of Optimal Interpolation, we found the optimal weight matrix W that minimizes the total analysis error variance.
hree-dimensional variational assimilation (3D-Var) In the derivation of Optimal Interpolation, we found the optimal weight matrix W that minimizes the total analysis error variance. Lorenc (1986) showed
More informationAnalysis of multicomponent walkaway vertical seismic profile data
Analysis of multicomponent walkaway vertical seismic profile data Bona Wu, Don C. Lawton, and Kevin W. Hall ABSTRACT A multicomponent walkaway VSP data processed for PP and PS imaging as well to study
More informationUsing Curvature to Map Faults, Fractures
Using Curvature to Map Faults, Fractures by SATINDER CHOPRA and KURT J. MARFURT Editor s note: Chopra is with Arcis Corp., Calgary, Canada; Marfurt is with the University of Oklahoma. Both are AAPG members.
More informationLeast Squares with Examples in Signal Processing 1. 2 Overdetermined equations. 1 Notation. The sum of squares of x is denoted by x 2 2, i.e.
Least Squares with Eamples in Signal Processing Ivan Selesnick March 7, 3 NYU-Poly These notes address (approimate) solutions to linear equations by least squares We deal with the easy case wherein the
More information2012 SEG SEG Las Vegas 2012 Annual Meeting Page 1
Wei Huang *, Kun Jiao, Denes Vigh, Jerry Kapoor, David Watts, Hongyan Li, David Derharoutian, Xin Cheng WesternGeco Summary Since the 1990s, subsalt imaging in the Gulf of Mexico (GOM) has been a major
More informationDeconvolution imaging condition for reverse-time migration
Stanford Exploration Project, Report 112, November 11, 2002, pages 83 96 Deconvolution imaging condition for reverse-time migration Alejandro A. Valenciano and Biondo Biondi 1 ABSTRACT The reverse-time
More information3. ESTIMATION OF SIGNALS USING A LEAST SQUARES TECHNIQUE
3. ESTIMATION OF SIGNALS USING A LEAST SQUARES TECHNIQUE 3.0 INTRODUCTION The purpose of this chapter is to introduce estimators shortly. More elaborated courses on System Identification, which are given
More informationLithology-based model transform. Summary
cdf( ) pdf( ) 3D lithological inversion of geophysical data Leif H. Cox*, Glenn A. Wilson, TechnoImaging, Michael S. Zhdanov, University of Utah and TechnoImaging, and Julian Vrbancich, Defense Science
More informationCONTENTS NOTATIONAL CONVENTIONS GLOSSARY OF KEY SYMBOLS 1 INTRODUCTION 1
DIGITAL SPECTRAL ANALYSIS WITH APPLICATIONS S.LAWRENCE MARPLE, JR. SUMMARY This new book provides a broad perspective of spectral estimation techniques and their implementation. It concerned with spectral
More informationReprocessing strategy for shallower prospects from the available 3D data set Case history of Cambay Basin
P-378 Reprocessing strategy for shallower prospects from the available 3D data set Case history of Cambay Basin M.Singh*, Vikash Chandra, USD Pandey Summary The occurrence of non commercial hydrocarbons
More informationSummary. Introduction
: Red Lodge, Montana Jingqiu Huang*, University of Houston, and Joe Wong, University of Calgary. Summary We report on a suite of geophysical surveys conducted on glacial sediments near Red Lodge, Montana.
More informationNonstationary filters, pseudodifferential operators, and their inverses
Nonstationary filters, pseudodifferential operators, and their inverses Gary F. Margrave and Robert J. Ferguson ABSTRACT An inversion scheme for nonstationary filters is presented and explored. Nonstationary
More informationarxiv: v3 [physics.geo-ph] 9 Jan 2013
Short-time homomorphic wavelet estimation arxiv:1209.0196v3 [physics.geo-ph] 9 Jan 2013 Roberto Henry Herrera and Mirko van der Baan Department of Physics, University of Alberta, Edmonton T6G 2E1, CA E-mail:
More information