Commonly, seismic data processing procedures, such as stacking and prestack migration, require the ability to detect bad traces/shots and restore or replace them by interpolation, particularly when the seismic observations are noisy or there are malfunctioned components in the recording system. However, currently available trace/shot interpolation methods in the spatial or Fourier domain must deal with requirements such as evenly sampled traces/shots, infinite bandwidth of the signals, and linear seismic events. In this paper, we present a novel method, termed the E-S (eigenspace seismic) method, using principal component analysis (PCA) of the seismic signal to address the issue of reliable detection or interpolation of bad traces/shots. The E-S method assumes the existence of a correlation between the observed seismic entities, such as trace or shot gathers, making it possible to estimate one of these entities from all others for interpolation or seismic quality control. It first transforms a trace (or shot) gather into an eigenspace using PCA. Then in the eigenspace, it treats every trace as a point with its loading scores of PCA as its coordinates. Simple linear, bilinear, or cubic spline 1 dimensional (1D) interpolation is used to determine PCA loading scores for any arbitrary coordinate in the eigenspace, which are then used to construct an interpolated trace for the desired position in physical space. This E-S method works with either regular or irregular sampling and, unlike various other published methods, it is well-suited for band-limited seismic records with curvilinear reflection events. We developed related algorithms and applied these to processed synthetic and offshore seismic survey data with or without simulated noises to demonstrate their performance. By comparing the interpolated and observed seismic traces, we find that the E-S method can effectively assess the quality of the trace, and restore poor quality data by interpolation. The successful processing of synthetic and real data using the E-S method presented in this approach will be widely applicable to seismic trace/shot interpolation and seismic quality control.
References
[1]
Yilmaz, O. (2001) Seismic Data Analysis: Processing, Inversion, and Interpretation of Seismic Data. Society of Exploration Geophysicists, Tulsa, OK, 1028 p. https://doi.org/10.1190/1.9781560801580
[2]
Sun, R. and Wang, A. (1998) Removing Bad Traces in Shallow Seismic Exploration. Tao, 9, 183-196. https://doi.org/10.3319/TAO.1998.9.2.183(T)
[3]
Abma, R. and Kabir, N. (2005) Comparison of Interpolation Algorithms. The Leading Edge, 24, 984-989. https://doi.org/10.1190/1.2112371
[4]
Fomel, S. (2003) Seismic Reflection Data Interpolation with Differential Offset and Shot Continuation. Geophysics, 68, 733-744. https://doi.org/10.1190/1.1567243
[5]
Bardan, V. (1987) Trace Interpolation in Seismic Data Processing. Geophysical Prospecting, 35, 343-358. https://doi.org/10.1111/j.1365-2478.1987.tb00822.x
[6]
Spitz, S. (1991) Seismic Trace Interpolation in the f-x Domain. Geophysics, 56, 785-796. https://doi.org/10.1190/1.1443096
[7]
Gulunay, N. (2003) Seismic Trace Interpolation in the Fourier Transform Domain. Geophysics, 68, 355-369. https://doi.org/10.1190/1.1543221
[8]
Naghizadeh, M. and Sacchi, D.M. (2009) f-x Adaptive Seismic-Trace Interpolation. Geophysics, 74, 9-16. https://doi.org/10.1190/1.3008547
[9]
Spitz, S. (1990) 3-D Seismic Interpolation in the f-x-y Domain. SEG Technical Program Expanded Abstracts, 1641-1643. https://doi.org/10.1190/1.1890083
[10]
Wang Y. (2002) Seismic Trace Interpolation in the f-x-y Domain. Geophysics, 67, 1232-1239. https://doi.org/10.1190/1.1500385
[11]
Liu, B. and Sacchi, D.M. (2004) Minimum Weighted Norm Interpolation of Seismic Records. Geophysics, 69, 1560-1568. https://doi.org/10.1190/1.1836829
[12]
Gulunay, N. and Chambers, R.E. (1996) Unaliased f-k Domain Trace Interpolation (UFKI). Society of Exploration Geophysicists, 1461-1464. https://doi.org/10.1190/1.1826390
[13]
Claerbout, J. (1998) Multidimensional Recursive Filters via a Helix. Geophysics, 63, 1532-1541. https://doi.org/10.1190/1.1444449
[14]
Porsani, M. (1999) Seismic Trace Interpolation Using Half-Step Prediction Filters. Geophysics, 64, 1461-1467. https://doi.org/10.1190/1.1444650
[15]
Naghizadeh, M. and Sacchi, D.M. (2007) Multistep Autoregressive Reconstruction of Seismic Records. Geophysics, 72, 111-118. https://doi.org/10.1190/1.2771685
[16]
Naghizadeh, M. and Innanen, A.K. (2011) Seismic Data Interpolation Using a Fast Generalized Fourier Transform. Geophysics, 76, 1-10. https://doi.org/10.1190/1.3511525
[17]
Duijndam, A.J.W., Schonewille, M.A. and Hindriks, C.O. (1999) Reconstruction of Band-Limited Signals, Irregularly Sampled along One Spatial Direction. Geophysics, 64, 524-538. https://doi.org/10.1190/1.1444559
[18]
Zwartjes, P.M. and Sacchi, M. (2007) Fourier Reconstruction of Nonuniformly Sampled, Aliased Seismic Data. Geophysics, 72, 21-32. https://doi.org/10.1190/1.2399442
[19]
Guo, Z. (2017) Seismic Data Interpolation: The Arbitrary Sampled Fourier Transform Method. Master’s Dissertation, University of Calgary, Calgary, 66.
[20]
Yu, Z., Ferguson, J., Mcmechan, G. and Anno, P. (2007) Wavelet-Radon Domain Dealiasing and Interpolation of Seismic Data. Geophysics, 72, 41-49. https://doi.org/10.1190/1.2422797
[21]
Wang, J.M. and Pert, M. (2010) Seismic Data Interpolation by Greedy Local Radon Transform. Geophysics, 75, WB225-WB234. https://doi.org/10.1190/1.3484195
[22]
Herrmann, F.J. and Hennenfent, G. (2008) Non-Parametric Seismic Data Recovery with Curvelet Frames. Geophysical Journal International, 173, 233-248. https://doi.org/10.1111/j.1365-246X.2007.03698.x
[23]
Chen, H.W. and Chang, C.W. (1999) Implicit Noise Reduction and Trace Interpolation in the Wavefield Depth Extrapolation. Geophysical Research Letters, 26, 3705-3708. https://doi.org/10.1029/1999GL005421
[24]
Fu, L. (2004) Wave Field Interpolation in the Fourier Wave Field Extrapolation. Geophysics, 69, 257-264. https://doi.org/10.1190/1.1649393
[25]
Oropeza, V. and Sacchi, M. (2011) Simultaneous Seismic Data Denoising and Reconstruction via Multichannel Singular Spectrum Analysis. Geophysics, 76, 25 https://doi.org/10.1190/1.3552706
[26]
Ronen, J. (1987) Wave-Equation Trace Interpolation. Geophysics, 52, 973-984. https://doi.org/10.1190/1.1442366
[27]
Trad, D. (2003) Interpolation and Multiple Attenuation with Migration Operators. Geophysics, 68, 2043-2054. https://doi.org/10.1190/1.1635058
[28]
Martingson, D. and Hopper, J. (1992) Nonlinear Seismic Trace Interpolation. Geophysics, 57, 136-145. https://doi.org/10.1190/1.1443177
[29]
Mandelli, S., Lipari, V., Bestagini, P. and Tubaro, S. (2019) Interpolation and Denoising of Seismic Data Using Convolutional Neural Networks. arXiv:1901.07927v1.
[30]
Mikhailiuk, A. and Faul, A. (2018) Deep Learning Applied to Seismic Data Interpolation. 80th EAGE Conference and Exhibition 2018, Extended Abstract. https://doi.org/10.3997/2214-4609.201800918
[31]
Li, Q. and Dehler, S.A. (2015) Inverse Principal Component Analysis for Geophysical Survey Data Interpolation. Journal of Applied Geophysics, 115, 79-91. https://doi.org/10.1016/j.jappgeo.2015.02.010
[32]
Huang, W., Wang, R., Zhou, Y., Chen, Y. and Yang, R. (2016) Improved Principal Component Analysis for 3D Seismic Data Simultaneous Reconstruction and Denoising. SEG Technical Program Expanded Abstracts, 4102-4106. https://doi.org/10.1190/segam2016-13858769.1
[33]
Li, Q. and Dehler, S.A. (2012) E-x Method for Seismic Trace Interpolation, EGU General Assembly, Vienna, 22-27 April 2012, 6129.
[34]
Ulrych, T.J., Freire, S. and Siston, P. (1988) Eigenimage Processing of Seismic Sections. SEG Technical Program Expanded Abstracts, 1261-1265. https://doi.org/10.1190/1.1892508
[35]
Wang, B. and Lu, W. (2017) Accurate and Efficient Seismic Data Interpolation in the Principal Frequency Wavenumber Domain. Journal of Geophysics and Engineering, 14, 1475-1488. https://doi.org/10.1088/1742-2140/aa82dc
[36]
Wu, J. and Bai, M., (2018) Fast Principal Component Analysis for Stacking Seismic Data. Journal of Geophysics and Engineering, 15, 295-306. https://doi.org/10.1088/1742-2140/aa9f80
[37]
Nicholson, W.K. (2003) Linear Algebra with Applications. 4th Edition, McGraw-Hill Ryerson Limited, Canada, 481.
[38]
Jolliffe, I.T. (2002) Principal Component Analysis. 2nd Edition, Springer-Verlag, New York, 414.
[39]
Golub, G.H. and Kahan, W. (1965) Calculating the Singular Values and Pseudo-Inverse of a Matrix. Journal of the Society for Industrial and Applied Mathematics, Series B, Numerical Analysis, 2, 205-224. https://doi.org/10.1137/0702016
[40]
Li, Q., Dehler, S.A., Shimeld, J. and Chian, D. (2013) Constructing a 3D Sediment Velocity Cube and Extracting Its Features in the Deep Water Areas of the Arctic Ocean. Atlantic Geology, 49, 45.
