Enhancement of seismic data

12196,904

17/630654

July 31, 2019

Methods, systems, and computer-readable medium to perform operations including: generating a first time-frequency spectrum of a first seismic trace from an original seismic dataset; generating a second time-frequency spectrum of a second seismic trace from an enhanced seismic dataset, where the second seismic trace corresponds to the first seismic trace; and re-combining an amplitude spectrum of the first time-frequency spectrum and a phase spectrum of the second time-frequency spectrum to generate a third time-frequency spectrum of an output trace that corresponds to the first and second seismic traces.

Saudi Arabian Oil Company (Dhahran, SA); Trofimuk Institute of Petroleum Geology and Geophysics (Novosibirsk, RU)

Saudi Arabian Oil Company (Dhahran, SA), Trofimuk Institute of Petroleum Geology and Geophysics (Novosibirsk, RU)

1. A method comprising: generating a first time-frequency spectrum of a first seismic trace from an original seismic dataset received from a plurality of sensors, the original seismic dataset comprising a plurality of seismic signals of a wellsite subterranean region; generating a second time-frequency spectrum of a second seismic trace from a signal to noise ratio filtered seismic dataset, wherein the second seismic trace suppresses a noise component of the original seismic dataset; combining an amplitude spectrum of the first time-frequency spectrum and a phase spectrum of the second time-frequency spectrum, for each discrete frequency bin, and frame index, to generate a third time-frequency spectrum of an output trace, the third time-frequency spectrum of the output trace comprising time shifts between neighboring traces differentiating the plurality of seismic signals that integrates noise masked data from the first and second seismic traces; and controlling wellsite systems, using the third time-frequency spectrum, to perform field operations.

2. The method of claim 1 , further comprising: generating a time-based output trace from the third time-frequency spectrum.

3. The method of claim 2 , wherein generating the first time-frequency spectrum comprises using short-term Fourier transform (STFT) to generate the first time-frequency spectrum, and wherein generating the time-based output trace comprises using inverse short-term Fourier transform (ISTFT) to generate the time-based output trace.

4. The method of claim 1 , wherein the signal to noise ratio filtered seismic dataset is generated by performing a signal-to-noise ratio (SNR) enhancement procedure on the original seismic dataset.

5. The method of claim 4 , wherein the SNR enhancement procedure comprises: common-reflection surface method (CRS), multi-focusing (MF), supergrouping, or non-linear beamforming.

6. The method of claim 1 , wherein the signal to noise ratio filtered seismic dataset has an identical acquisition geometry to the original seismic dataset, and wherein the second seismic trace is located in the same position in the signal to noise ratio filtered seismic dataset as the first seismic trace in the original seismic dataset.

7. The method of claim 1 , wherein the original seismic dataset comprises a plurality of seismic traces, the method further comprising: for each seismic trace of the plurality, generating a respective output trace corresponding to the seismic trace, wherein the respective output traces collectively form an output seismic dataset.

8. A device comprising: one or more processors; and a non-transitory computer-readable storage medium coupled to the one or more processors and storing programming instructions for execution by the one or more processors, the programming instructions instructing the one or more processors to perform operations comprising: generating a first time-frequency spectrum of a first seismic trace from an original seismic dataset received from a plurality of sensors, the original seismic dataset comprising a plurality of seismic signals of a wellsite subterranean region; generating a second time-frequency spectrum of a second seismic trace from a signal to noise ratio filtered seismic dataset, wherein the second seismic trace suppresses a noise component of the original seismic dataset; combining an amplitude spectrum of the first time-frequency spectrum and a phase spectrum of the second time-frequency spectrum, for each discrete frequency bin, and frame index, to generate a third time-frequency spectrum of an output trace, the third time-frequency spectrum of the output trace comprising time shifts between neighboring traces differentiating the plurality of seismic signals that integrates noise masked data from the first and second seismic traces; and controlling wellsite systems, using the third time-frequency spectrum, to perform field operations.

9. The device of claim 8 , the operations further comprising: generating a time-based output trace from the third time-frequency spectrum.

10. The device of claim 9 , wherein generating the first time-frequency spectrum comprises using short-term Fourier transform (STFT) to generate the first time-frequency spectrum, and wherein generating the time-based output trace comprises using inverse short-term Fourier transform (ISTFT) to generate the time-based output trace.

11. The device of claim 8 , wherein the signal to noise ratio filtered seismic dataset is generated by performing a signal-to-noise ratio (SNR) enhancement procedure on the original seismic dataset.

12. The device of claim 11 , wherein the SNR enhancement procedure comprises: common-reflection surface method (CRS), multi-focusing (MF), supergrouping, or non-linear beamforming.

13. The device of claim 8 , wherein the signal to noise ratio filtered seismic dataset has an identical acquisition geometry to the original seismic dataset, and wherein the second seismic trace is located in the same position in the signal to noise ratio filtered seismic dataset as the first seismic trace in the original seismic dataset.

14. The device of claim 8 , wherein the original seismic dataset comprises a plurality of seismic traces, and wherein the operations further comprise: for each seismic trace of the plurality, generating a respective output trace corresponding to the seismic trace, wherein the respective output traces collectively form an output seismic dataset.

15. A non-transitory computer-readable medium storing instructions executable by a computer system to perform operations comprising: generating a first time-frequency spectrum of a first seismic trace from an original seismic dataset received from a plurality of sensors, the original seismic dataset comprising a plurality of seismic signals of a wellsite subterranean region; generating a second time-frequency spectrum of a second seismic trace from a signal to noise ratio filtered seismic dataset, wherein the second seismic trace suppresses a noise component of the original seismic dataset; combining an amplitude spectrum of the first time-frequency spectrum and a phase spectrum of the second time-frequency spectrum, for each discrete frequency bin, and frame index, to generate a third time-frequency spectrum of an output trace, the third time-frequency spectrum of the output trace comprising time shifts between neighboring traces differentiating the plurality of seismic signals that integrates noise masked data from the first and second seismic traces; and controlling wellsite systems, using the third time-frequency spectrum, to perform field operations.

16. The non-transitory computer-readable medium of claim 15 , the operations further comprising: generating a time-based output trace from the third time-frequency spectrum.

17. The non-transitory computer-readable medium of claim 16 , wherein generating the first time-frequency spectrum comprises using short-term Fourier transform (STFT) to generate the first time-frequency spectrum, and wherein generating the time-based output trace comprises using inverse short-term Fourier transform (ISTFT) to generate the time-based output trace.

18. The non-transitory computer-readable medium of claim 15 , wherein the signal to noise ratio filtered seismic dataset is generated by performing a signal-to-noise ratio (SNR) enhancement procedure on the original seismic dataset.

19. The non-transitory computer-readable medium of claim 18 , wherein the SNR enhancement procedure comprises: common-reflection surface method (CRS), multi-focusing (MF), supergrouping, or non-linear beamforming.

20. The non-transitory computer-readable medium of claim 15 , wherein the signal to noise ratio filtered seismic dataset has an identical acquisition geometry to the original seismic dataset, and wherein the second seismic trace is located in the same position in the signal to noise ratio filtered seismic dataset as the first seismic trace in the original seismic dataset.

3597727 August 1971 Judson et al.
3806863 April 1974 Brown et al.
5850622 December 1998 Vassiliou et al.
6161076 December 2000 Barr et al.
6876599 April 2005 Combee
6904368 June 2005 Reshef et al.
7379386 May 2008 Muyzert et al.
7453764 November 2008 Cavalca et al.
9470811 October 2016 Sun et al.
10947789 March 2021 Al-Abdulrahman et al.
2003/0055567 March 2003 Kelly
2004/0233782 November 2004 Herkenhoff et al.
2005/0075790 April 2005 Taner
2005/0090989 April 2005 Kelly
2007/0055465 March 2007 Lecerf
2008/0049551 February 2008 Muyzert et al.
2009/0276159 November 2009 Strobbia
2010/0027376 February 2010 Washbourne et al.
2011/0013481 January 2011 Clark
2011/0295510 December 2011 Gulati
2017/0269246 September 2017 Anderson et al.
2020/0264327 August 2020 Gadylshin et al.
2020/0400847 December 2020 Bakulin et al.
2023/0266492 August 2023 Bakulin et al.
102109612 June 2011
102478671 May 2012
102854530 January 2013
104268412 January 2015
105388519 March 2016
106443787 February 2017
106646613 May 2017
106842297 June 2017
107340541 November 2017
2762926 August 2014
WO 2005017564 February 2005
WO 2011034870 March 2011
WO 2013190512 December 2013
WO 2019156586 August 2019
WO 2021020982 February 2021
WO 2021020983 February 2021
WO 2021020984 February 2021

Mosher, Charles C. “Generalized windowed transforms for seismic processing and imaging.” SEG Technical Program Expanded Abstracts 2012, Sep. 2012, https://doi.org/10.1190/segam2012-1196.1. (Year: 2012). cited by examiner
U.S. Appl. No. 17/595,365, Bakulin et al., Nov. 15, 2021. cited by applicant
U.S. Appl. No. 17/630,663, Bakulin et al., Jan. 27, 2022. cited by applicant
U.S. Appl. No. 17/630,665, Bakulin et al., Jan. 27, 2022. cited by applicant
Alsteris et al., “Short-time phase spectrum in speech processing: a review and some experimental results,” Digital Signal Processing, 2007, 17:578-616, 39 pages. cited by applicant
Bakulin et al., “Application of supergrouping to enhance 3D prestack seismic data from a desert environment,” The Leading Edge, Mar. 2018, 201-207, 8 pages. cited by applicant
Bakulin et al., “Application of supergrouping to land seismic data in desert environment,” SEG International Exposition and 86th Annual Meeting, 2016, 5 pages. cited by applicant
Bakulin et al., “Nonlinear beaforming for enhancing prestack seismic data with a challenging near surface or overburden,” Data Processing & Machine Learning, Dec. 2018, 6 pages. cited by applicant
Bakulin et al., “Nonlinear beamforming for enhancing pre-stack data with challenging near surface or overburden,” SEG Technical Program Expanded Abstracts, Aug. 2017, 5 pages. cited by applicant
Baykulov et al., “Prestack seismic data enhancement with partial common-reflection-surface (CRS) stack,” Geophysics, May-Jun. 2009, 74(3):49-58, 10 pages. cited by applicant
Berkovitch et al., “How non-hyperbolic MultiFocusing improves depth imaging,” First Break, Sep. 2011, 29:95-103, 9 pages. cited by applicant
Billete et al., “Practical aspects and applications of 2D stereotomography,” Geophysics, May 2003, 68(3):1008-1021, 14 pages. cited by applicant
Bishop et al., “Tomographic determination of velocity and depth in laterally varying media,” Geophysics, Jun. 1985, 50(6):903-923, 22 pages. cited by applicant
Buzlukov et al., “Imaging improvement by prestack signal enhancement,” Geophysical Prospecting, 2013, 61:1150-1158, 9 pages. cited by applicant
Buzlukov et al., “Prestack data enhancement using local traveltime approximation,” 72nd EAGE Conference & Exhibition incorporating EUROPEC, Barcelona, Spain, Jun. 14-17, 2010, 5 pages. cited by applicant
Curia et al., “The impact of multifocusing in the processing of land 3D seismic data in a fold and thrust belt setting: Ranquil Norte Block, Neuquén Basin, Argentina,” The Leading Edge, 2013, 36(9):770-774, 5 pages. cited by applicant
Digital Signal Processing, Blackledget, Horwood Publishing, 2006, 2nd ed., 840 pages. cited by applicant
Duveneck, “Tomographic determination of seismic velocity models with kinematic wavefield attributes,” Zur Erlangung des akademischen Grades eines, Doktors De Naturwissenschaften von der Fakultat fur Physik der Universitat Karlsruhe (TH) genehmigte, DISSERTATION von Dipl.-Geophys., Jun. 2004, 201 pages. cited by applicant
Duveneck, “Tomographic velocity model estimation with data-derived first and second travel time derivatives,” presented at the 8th International Congress of the Brazilian Geophysical Society, Brazil, Sep. 14-18, 2003, 4 pages. cited by applicant
Duveneck, “Velocity model estimation with data-derived wavefront attributes,” Geophysics, Jan. 2004, 69(1):265-274, 10 pages. cited by applicant
Erdogan et al., “Phase-sensitive and recognition-boosted speech separation using deep recurrent neural networks,” Proc. Int. Conf. Acoust., Speech, Signal Process., 2015, 708-712, 5 pages. cited by applicant
Fagin, “Model based depth imaging,” Society of Exploration Geophysics, Course Notes Series No. 10, Tulsa, USA, 1999, 184 pages. cited by applicant
Gjoystdal et al., “Inversion of reflection times in three dimensions,” Geophysics, Jul. 1981, 46(7):972-983, 12 pages. cited by applicant
Hoecht et al., “Operator-oriented CRS interpolation,” Geophysical Prospecting, 2009, 57:957-979, 23 pages. cited by applicant
Hubral et al., “Interval Velocities from Seismic Reflection Time Measurements,” Society of Exploration Geophysicists, 1980, 214 pages. cited by applicant
Jones, “Estimating subsurface parameter fields for seismic migration: Velocity model building,” Seismic velocity/statics, Oct. 2015, 24 pages. cited by applicant
Khaidukov et al., “Enhancement of seismic data gathered by floating ice acquisition: application of local kinematic attributes,” SEG International Exposition and 86th Annual Meeting, 2016, 5 pages. cited by applicant
Klüver, “Velocity model building using migration to residual time,” SEG Technical Program Expanded Abstracts, SEG/New Orleans 2006 Annual Meeting, 5 pages. cited by applicant
Lambaré et al., “Multi-layer non-linear slope tomography,” 75th Annual International Conference and Exhibition, EAGE, Extended Abstracts, London, Jun. 10-13, 2013, 5 pages. cited by applicant
Lambaré, “Stereotomography,” Geophysics, Sep. 2008, 73(5):VE25-VE34, 10 pages. cited by applicant
Lavaud et al., “Automatic robust velocity estimation by poststack stereotomography,” SEG Technical Program Expanded Abstracts 2004, SEG International Exposition and 74th Annual Meeting, Oct. 10-15, 2004, 4 pages. cited by applicant
Li et al., “A low rank approximation method for ice-break noise attenuation of simultaneous vibroseis data,” SEG Technical Program Expanded Abstracts, Aug. 2017, 5 pages. cited by applicant
Liang et al., The optimal ratio time-frequency mask for speech separation in terms of the signal-to-noise ratio, The Journal of the Acoustical Society of America, 2013, 134(5):EL452, 8 pages. cited by applicant
Lichman et al., “Automated Phase-based Moveout Correction,” 69th SEG Annual International Meeting, Expanded Abstracts, 2000, 5 pages. cited by applicant
Lichman et al., “Phase Inversion Deconvolution for Surface Consistent Processing and Multiple Attenuation,” SEG Annual International Meeting, Expanded Abstracts, 1999, 4 pages. cited by applicant
Martin et al., “Noise Power Spectral Density Estimation Based on Optimal Smoothing and Minimum Statistics,” IEEE Transactions on Speech and Audio Processing, Jul. 2001, 9(5):504-512, 9 pp. 512. cited by applicant
Mueller et al., “Improving Prestack Migration with CRS Techniques—a Case Study,” 72nd EAGE Conference and Exhibition, Expanded Abstract, Jun. 2010, 5 pages. cited by applicant
Neckludov et al., “Residual stereotomographic inversion,” Geophysics, Jul. 2006, 71(4):E35-E39, 6 pages. cited by applicant
Neklyudov et al., “Intra-array statics and phase corrections obtained by beamforming in the short-time fourier transform domain: application to supergrouping,” SEG International Exposition, 2017, 6 pages. cited by applicant
Robinson et al., “Geophysical signal analysis,” Society of Exploration Geophysicists, 2000, 481 pages. cited by applicant
Tanushev et al., “Fast, high-resolution beam tomography and velocity-model building,” The Leading Edge, Feb. 2017, 36(2):140-145, 6 pages. cited by applicant
Ulrych et al., “The essence of phase in seismic data processing and inversion,” 77th SEG Annual Meeting, Expanded Abstracts, 2007, 1765-1769, 5 pages. cited by applicant
Van Der Made et al., “Estimation of complex velocity models from stacking information,” SEG Technical program Expanded Abstracts, 1987, 3 pages. cited by applicant
Wang et al., “Time-frequency masking for speech separation and its potential for hearing aid design,” Trends in Amplificat., 2008, 12:332-353, 22 pages. cited by applicant
Wikipedia.com [online], “Median filter,” Last edited Sep. 10, 2021 retrieved Apr. 8, 2022, retrieved from URL , 4 pages. cited by applicant
Williamson et al., “Time-Frequency Masking in the Complex Domain for Speech Dereverberation and Denoising,” IEEE/ACM Transactions on Audio, Speech and Language Processing, May 2016, 13(9):1-10, 10 pages. cited by applicant
Woodward et al., “A decade of tomography,” Geophysics, Sep.-Oct. 2008, 73(5), 7 pages. cited by applicant
Yilmaz et al., “Blind separation of speech mixtures via time-frequency masking,” IEEE Trans. Signal Process, Jul. 2004, 52(7):1830-1846, 19 pages. cited by applicant
Zhang et al., “Common-reflection-surface (CRS) stack for common offset,” Geophysical Prospecting, 2001, 49:709-18, 10 pages. cited by applicant
PCT International Search Report and Written Opinion in International Appln. No. PCT/RU2019/000336, dated Feb. 26, 2020, 14 pages. cited by applicant
PCT International Search Report and Written Opinion in International Appln. No. PCT/RU2018/000079, dated Dec. 6, 2018, 12 pages. cited by applicant
PCT International Search Report and Written Opinion in International Appln. No. PCT/RU2019/00532, dated Apr. 3, 2020, 17 pages. cited by applicant
PCT International Search Report and Written Opinion in International Appln. No. PCT/RU2019/000534, dated Apr. 3, 2020, 17 pages. cited by applicant
PCT International Search Report and Written Opinion in International Appln. No. PCT/RU2019/000533, dated Apr. 8, 2020, 13 pages. cited by applicant
Bakulin et al., “3D data-domain reflection tomography for intial velocity model building using challenging 3D seismic data,” Journal of Seismic Exploration, 2021, 30:419-446, 28 pages. cited by applicant
Stork, “Reflection tomography in the postmigrated domain,” Geophysics, May 1992, 57(5):680-692, 13 pages. cited by applicant
RU Office Action in Russian Appln. No. 2021135462, dated Sep. 20, 2022, 12 pages, with English Translation. cited by applicant
RU Office Action in Russian Appln. No. 2022103676, dated Jan. 18, 2023, 10 pages, with English Translation. cited by applicant
RU Office Action in Russian Appln. No. 2022104308, dated Jan. 12, 2023, 11 pages, with English Translation. cited by applicant
RU Office Action in Russian Appln. No. 2022104454, dated Jan. 16, 2023, 10 pages, with English Translation. cited by applicant
Bakulin et al., “3D Prestack Land Data Enhancement Based on Nonlinear Beamforming in Cross-Spread Domain,” EAGE 81st Annual Conference and Exhibition, Jun. 6, 2019, 5 pages. cited by applicant
Bakulin et al., “Bootstrapping invisible signals: prestack land data enhancement using nonlinear beamforming with local waveform corrections,” SEG Expanded Abstracts, 2019, 3954-3958, 5 pages. cited by applicant
Bakulin et al., “Enhance—Estimate—Image: new processing approach for single-sensor and other seismic data with low pre-stack signal-to-noise ratio,” SEG Technical Program Expanded Abstracts, Aug. 17, 2017, 5 pages. cited by applicant
Bakulin et al., “Nonlinear beamforming for enhancement of 3D prestack land seismic data,” Geophysics, May 2020, 85(3):V283-V296, 14 pages. cited by applicant
Bakulin et al., “Prestack data enhancement with phase corrections in time-frequency domain guided by local multidimensional stacking,” Geophysical Prospecting, May 15, 2020, 68(6):1811-1818, 8 pages. cited by applicant
Billette et al., “Practical aspects and applications of 2D stereotomography,” Geophysics, May-Jun. 2003, 68(3):1008-1021, 14 pages. cited by applicant
Gadylshin et al., “Inpainting of local wavefront attributes using artificial intelligence,” SEG Technical Program Expanded Abstracts, Aug. 1, 2019, 2212-2216, 5 pages. cited by applicant
Neckludov et al., “Residual stereotomographic inversion,” Geophysics, Jul.-Aug. 2006, 71(4): E35-E39, 5 pages. cited by applicant
Özbek et al., “3-D filter design on a hexagonal grid with applications to point-receiver land acquisition,” SEG Int'l Exposition and 74th Annual Meeting, Oct. 10-15, 2004, 5 pages. cited by applicant
PCT International Search Report and Written Opinion in International Appln. No. PCT/RU2020/000469, dated May 28, 2021, 18 pages. cited by applicant
RU Office Action in Russian Appln. No. 2023106594, dated Nov. 30, 2023, 17 pages, with English Translation. cited by applicant
SAIP Examination Report in SAIP Appln. No. 521430857, dated Dec. 22, 2023, 11 pages with English translation. cited by applicant
SAIP Examination Report in SAIP Appln. No. 522431566, dated Dec. 27, 2023, 15 pages. cited by applicant

edspgr.12196904