IP Library Granted Patent US 12,360,271
Granted Patent B2
US 12,360,271 · App. 17/728,718 · Granted Jul 15, 2025

Methods and systems for generating an image of a subterranean formation based on low frequency reconstructed seismic data

Inventor: Maiza Bekara (Kingston Upon Thames, GB)
Assignee: PGS Geophysical AS
G01V1/301G01V1/32G01V1/34G01V1/3808G01V1/3843G01V20/00G01V2210/65G01V2210/74
View Patent ↗
Loading inventors, assignments & file history…
Monitor This Case
Get email alerts when status or documents change.
Order Certified Copies
Most orders are placed with the USPTO same day — all within 24 business hours.
Order via The Patent Place →
Pre-filled with this patent's details
Quick Facts
Patent No.
US 12,360,271
App. No.
17/728,718
Granted
Jul 15, 2025
Kind
B2
Abstract

This disclosure presents processes and systems for generating an image of a subterranean formation from seismic data recorded in a seismic survey of the subterranean formation. The seismic data is contaminated with low frequency noise in a low frequency band. Processes and systems reconstruct seismic data in the low frequency band of the seismic data to obtain low frequency reconstructed seismic data that is free of the low frequency noise. The low frequency reconstructed seismic data is used to construct a velocity model of the subterranean formation. The velocity model and the low frequency reconstructed seismic data are used to generate an image of the subterranean formation that reveals structures of the subterranean formation without contamination from the low frequency noise.

Claims (98)

1. In a computer-implemented process for generating an image of a subterranean formation from seismic data recorded in a seismic survey of the subterranean formation, the improvement comprising:

using an autoregressive model to reconstruct seismic data in a low frequency band of the seismic data from signals in the seismic data with relatively higher frequencies than the low frequency band to obtain low frequency reconstructed seismic data that is free of low frequency noise in the low frequency band;

constructing a velocity model of the subterranean formation based on the low frequency reconstructed seismic data; and

generating an image of the subterranean formation based on the velocity model and the low frequency reconstructed seismic data, the image revealing structures of the subterranean formation without contamination from the low frequency noise.

2. The process of claim 1 wherein using the autoregressive model to reconstruct seismic data in the low frequency band comprises:

for each location of a sliding space-time window of the seismic data,

transforming seismic data in the space-time window from a space-time domain to an intermediate domain, wherein the intermediate domain is one of a frequency-slowness domain, a frequency-velocity domain, and a frequency-wavenumber domain,

using the autoregressive model to reconstruct seismic data in the low frequency band of each trace in the intermediate domain, and

transforming the low frequency reconstructed traces to the space-time domain to obtain low frequency reconstructed seismic data in the space-time window.

3. The process of claim 1 wherein using the autoregressive model to reconstruct seismic data in the low frequency band comprises:

for each location of a sliding space-time window of the seismic data,

transforming seismic data in the space-time window from the space-time domain to an intermediate domain, wherein the intermediate domain is one of a frequency-slowness domain, a frequency-velocity domain, and a frequency-wavenumber domain, and

for each trace in the intermediate domain,

detecting amplitudes with frequencies in the low frequency band,

reconstructing the amplitudes in the low frequency band to obtain a low frequency reconstructed trace using the autoregressive model, and

replacing the trace with the low frequency reconstructed trace.

4. The process of claim 1 wherein using the autoregressive model to reconstruct seismic data in the low frequency band comprises:

transforming the seismic data from the time domain to obtain transformed seismic data in an intermediate domain, wherein the intermediate domain is one of a frequency-slowness domain, a frequency-velocity domain, and a frequency-wavenumber domain;

detecting amplitudes of the transformed seismic data with frequencies in a low frequency band;

reconstructing amplitudes of traces in the low frequency band of the intermediate domain using the autoregressive model to obtain low frequency reconstructed seismic data in the intermediate domain; and

inverse transforming the low frequency reconstructed seismic data in the space time domain.

5. The process of claim 1 wherein constructing the velocity model of the subterranean formation based on the low frequency reconstructed seismic data comprises performing full-waveform inversion with the low frequency reconstructed seismic data to obtain the velocity model free of contamination from the low frequency noise.

6. The process of claim 1 wherein the autoregressive model is one of a linear autoregressive model and a non-linear autoregressive model.

7. A computer system for generating an image of a subterranean formation from seismic data recorded in a seismic survey of the subterranean formation, the system comprising:

one or more processors;

one or more data-storage devices; and

machine-readable instructions stored in the one or more data-storage devices that when executed using the one or more processors controls the system to perform operations comprising:

using an autoregressive model to reconstruct seismic data in a low frequency band of the seismic data from signals in the seismic data with relatively higher frequencies than the low frequency band to obtain low frequency reconstructed seismic data that is free of low frequency noise in the low frequency band;

constructing a velocity model of the subterranean formation based on the low frequency reconstructed seismic data; and

generating an image of the subterranean formation based on using the velocity model and the low frequency reconstructed seismic data, thereby the image revealing structures of the subterranean formation without contamination from the low frequency noise.

8. The computer system of claim 7 wherein using the autoregressive model to reconstruct seismic data in the low frequency band comprises:

for each location of a sliding space-time window of the seismic data,

transforming seismic data in the space-time window from the space-time domain to an intermediate domain, wherein the intermediate domain is one of a frequency-slowness domain, a frequency-velocity domain, and a frequency-wavenumber domain,

using the autoregressive model to reconstruct seismic data in the low frequency band of each trace in the intermediate domain, and

transforming the low frequency reconstructed traces to the space-time domain to obtain low frequency reconstructed seismic data in the space-time window.

9. The computer system of claim 7 wherein using the autoregressive model to reconstruct seismic data in the low frequency band comprises:

for each location of a sliding space-time window of the seismic data,

transforming seismic data in the space-time window from the space-time domain to an intermediate domain, wherein the intermediate domain is one of a frequency-slowness domain, a frequency-velocity domain, and a frequency-wavenumber domain, and

for each trace in the intermediate domain,

detecting amplitudes with frequencies in the low frequency band,

reconstructing the amplitudes in the low frequency band to obtain a low frequency reconstructed trace using the autoregressive model, and

replacing the trace with the low frequency reconstructed trace.

10. The computer system of claim 7 wherein using the autoregressive model to reconstruct seismic data in the low frequency band comprises:

transforming the seismic data from the time domain to obtain transformed seismic data an intermediate domain, wherein the intermediate domain is one of an intercept time-slowness domain, a frequency-slowness domain, a frequency-velocity domain, and a frequency-wavenumber domain;

detecting amplitudes of the transformed seismic data with frequencies in a low frequency band;

reconstructing amplitudes of traces in the low frequency band of the intermediate domain using the autoregressive model to obtain low frequency reconstructed seismic data in the intermediate domain; and

inverse transforming the low frequency reconstructed seismic data in the space time domain.

11. The computer system of claim 7 wherein constructing the velocity model of the subterranean formation based on the low frequency reconstructed seismic data comprises performing full-waveform inversion with the low frequency reconstructed seismic data to obtain the velocity model free of contamination from the low frequency noise.

12. The computer system of claim 7 wherein the autoregressive model is one of a linear autoregressive model and a non-linear autoregressive model.

13. A non-transitory computer-readable medium having machine-readable instructions encoded thereon for enabling one or more processors of a computer system to perform operations comprising:

reconstructing seismic data in a low frequency band of the seismic data from signals in the seismic data with relatively higher frequencies than the low frequency band to obtain low frequency reconstructed seismic data that is free of low frequency noise in the low frequency band;

constructing a velocity model of a subterranean formation based on the low frequency reconstructed seismic data; and

generating an image of the subterranean formation based on the velocity model and the low frequency reconstructed seismic data.

14. The medium of claim 13 wherein reconstructing seismic data in the low frequency band comprises:

for each location of a sliding space-time window of the seismic data,

transforming seismic data in the space-time window from the space-time domain to an intermediate domain, wherein the intermediate domain is one of a frequency-slowness domain, a frequency-velocity domain, and a frequency-wavenumber domain,

using an autoregressive model to reconstruct seismic data in the low frequency band of each trace in the intermediate domain, and

transforming the low frequency reconstructed traces to the space-time domain to obtain low frequency reconstructed seismic data in the space-time window.

15. The medium of claim 13 wherein reconstructing seismic data in the low frequency band comprises:

for each location of a sliding space-time window of the seismic data,

transforming seismic data in the space-time window from the space-time domain to an intermediate domain, wherein the intermediate domain is one of a frequency-slowness domain, a frequency-velocity domain, and a frequency-wavenumber domain, and

for each trace in the intermediate domain,

detecting amplitudes with frequencies in the low frequency band,

reconstructing the amplitudes in the low frequency band to obtain a low frequency reconstructed trace using an autoregressive model, and

replacing the trace with the low frequency reconstructed trace.

16. The medium of claim 13 wherein reconstructing seismic data in the low frequency band comprises:

transforming the seismic data from the time domain to obtain transformed seismic data an intermediate domain, wherein the intermediate domain is one of a frequency-slowness domain, a frequency-velocity domain, and a frequency-wavenumber domain;

detecting amplitudes of the transformed seismic data with frequencies in a low frequency band;

reconstructing amplitudes of traces in the low frequency band of the intermediate domain using an autoregressive model to obtain low frequency reconstructed seismic data in the intermediate domain; and

inverse transforming the low frequency reconstructed seismic data in the space time domain.

17. The medium of claim 13 wherein constructing the velocity model of the subterranean formation based on the low frequency reconstructed seismic data comprises performing full-waveform inversion with the low frequency reconstructed seismic data to obtain the velocity model free of contamination from the low frequency noise.

18. Apparatus for generating an image of a subterranean formation from seismic data recorded in a seismic survey of the subterranean formation, the apparatus comprising:

means for reconstructing seismic data in a low frequency band of the seismic data from signals in the seismic data with relatively higher frequencies than the low frequency band to obtain low frequency reconstructed seismic data that is free of low frequency noise in the low frequency band;

means for constructing a velocity model of the subterranean formation based on the low frequency reconstructed seismic data; and

means for generating the image of the subterranean formation based on the velocity model and the low frequency reconstructed seismic data.

19. The apparatus of claim 18 wherein means for reconstructing seismic data in the low frequency band comprises:

for each location of a sliding space-time window of the seismic data,

transforms seismic data in the space-time window from the space-time domain to an intermediate domain, wherein the intermediate domain is one of a frequency-slowness domain, a frequency-velocity domain, and a frequency-wavenumber domain,

uses an autoregressive model to reconstruct seismic data in the low frequency band of each trace in the intermediate domain, and

transforms the low frequency reconstructed traces to the space-time domain to obtain low frequency reconstructed seismic data in the space-time window.

20. The apparatus of claim 18 wherein means for reconstructing seismic data in the low frequency band comprises:

for each location of a sliding space-time window of the seismic data,

transforms seismic data in the space-time window from the space-time domain to an intermediate domain, wherein the intermediate domain is one of a frequency-slowness domain, a frequency-velocity domain, and a frequency-wavenumber domain, and

for each trace in the intermediate domain,

detects amplitudes with frequencies in the low frequency band,

reconstructs the amplitudes in the low frequency band to obtain a low frequency reconstructed trace using an autoregressive model, and

replaces the trace with the low frequency reconstructed trace.

21. The apparatus of claim 18 wherein means for reconstructing seismic data in the low frequency band comprises:

transforms the seismic data from the time domain to obtain transformed seismic data an intermediate domain, wherein the intermediate domain is one of a frequency-slowness domain, a frequency-velocity domain, and a frequency-wavenumber domain;

detects amplitudes of the transformed seismic data with frequencies in a low frequency band;

reconstructs amplitudes of traces in the low frequency band of the intermediate domain using an autoregressive model to obtain low frequency reconstructed seismic data in the intermediate domain; and

inverse transforms the low frequency reconstructed seismic data in the space time domain.

22. The apparatus of claim 18 wherein constructing the velocity model of the subterranean formation based on the low frequency reconstructed seismic data comprises performing full-waveform inversion with the low frequency reconstructed seismic data to obtain the velocity model free of contamination from the low frequency noise.

23. A method for manufacturing a geophysical data product, the method comprising:

reconstructing seismic data in a low frequency band of the seismic data from signals in the seismic data with relatively higher frequencies than the low frequency band using an autoregression model to obtain low frequency reconstructed seismic data that is free of low frequency noise in the low frequency band;

constructing a velocity model of a subterranean formation based on the low frequency reconstructed seismic data;

generating an image of the subterranean formation based on the velocity model and the low frequency reconstructed seismic data; and

storing the velocity model and the image in a computer-readable medium.

Assignments (1)
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Apr 25, 2022
From: BEKARA, MAIZA
To: PGS GEOPHYSICAL AS
Reel/Frame 059701/0979 →
Continuity (2)
Provisional Application 63183367 · May 3, 2021
Related Publication 20220373703A1 · Nov 24, 2022
References Cited (19)
US 10224900B2 · Poole · 2019 [cited by examiner]
US 11035968B2 · Li · 2021 [cited by examiner]
US 20110103182A1 · Martin · 2011 [cited by examiner]
US 20120081999A1 · Halliday · 2012 [cited by examiner]
US 20120224454A1 · Soubaras · 2012 [cited by examiner]
US 20150346368A1 · Dellinger · 2015 [cited by examiner]
US 20210132242A1 · Sallas · 2021 [cited by examiner]
US 20210223424A1 · Valensi · 2021 [cited by examiner]
US 20210318458A1 · Baumstein · 2021 [cited by examiner]
CN 106291661A · 2017 [cited by applicant]
WO 2016154411A1 · 2016 [cited by applicant]
Bianchin et al. (Acoustic impedance estimation from combined harmonic reconstruction and interval velocity, 2019, Geophysics) (Year: 2019). [cited by examiner]
Chen et al. (Application of Envelope in Salt Structure Velocity Building: From Objective Function Construction to the Full-Band Seismic Data Reconstruction, IEEE,2020). (Year: 2020). [cited by examiner]
Mostafa et al. (Multistep autoregressive reconstruction of seismic records, 2007, Geophysics) (Year: 2007). [cited by examiner]
Multistep autoregressive reconstruction of seismic records Mostafa Naghizadeh1 and Mauricio D. Sacchi1 (Year: 2007). [cited by examiner]
Oleg Ovcharenko, Vladimir Kazei, Daniel Peter, and Tariq Alkhalifah. (2017). Neural network based low-frequency data extrapolation. 3 pages. [cited by applicant]
Scott Keating and Kris Innanen. CREWES Research Report—vol. 27 (2015). Extrapolation of low frequencies and application to physical modeling data. 14 pages. [cited by applicant]
Oleg Ovcharenko, Vladimir Kazei, Mahesh Kalita, and Daniel Peter. (2019). Deep learning for low-frequency extrapolation from multi-offset seismic data. 54 pages. [cited by applicant]
Hongyu Sun and Laurent Demanet. (2018). Low frequency extrapolation with deep learning. 5 pages. [cited by applicant]