US 11,656,378 B2
Seismic imaging by visco-acoustic reverse time migration
Hongwei Liu, Dhahran (SA); and Yi Luo, Dhahran (SA)
Assigned to Saudi Arabian Oil Company, Dhahran (SA)
Filed by Saudi Arabian Oil Company, Dhahran (SA)
Filed on Jun. 8, 2020, as Appl. No. 16/895,762.
Prior Publication US 2021/0382192 A1, Dec. 9, 2021
Int. Cl. G01V 1/30 (2006.01); G01V 1/18 (2006.01); G01V 1/28 (2006.01); G01V 1/34 (2006.01)
CPC G01V 1/301 (2013.01) [G01V 1/181 (2013.01); G01V 1/282 (2013.01); G01V 1/345 (2013.01); G01V 2210/14 (2013.01); G01V 2210/43 (2013.01); G01V 2210/642 (2013.01)] 17 Claims
OG exemplary drawing
 
1. A method for generating a seismic image representing a subsurface formation by compensating for distortions of seismic data, the method comprising:
receiving seismic data for the subsurface formation, including receiver wavelet data and source wavelet data;
generating source wavefield data based on a forward modeling of the source wavelet data;
generating receiver wavefield data that compensates for distortions in the seismic data by:
applying a dispersion-only model to the receiver wavelet data to generate a first reconstructed back-propagated receiver wavefield portion, wherein the dispersion-only model comprises applvina a wave eauation to the receiver wavelet, the wave equation comprising:

OG Complex Work Unit Math
wherein p represents an acoustic pressure value, wherein ν0 represents an initial velocity value, wherein ω0 represents an initial angular frequency value, wherein γ represents a constant value, wherein k represents a wavenumber, and wherein t represents time;
applying a dissipation-only model to the receiver wavelet data to generate a second reconstructed back-propagated receiver wavefield portion; and
combining the first back-propagated receiver wavefield portion and the second back-propagated receiver wavefield portion into the receiver wavefield data;
applying an imaging condition to the receiver wavefield data and the source wavefield data;
generating, based on applying the imaging condition, visco-acoustic reverse time migration (VARTM) result data;
based on the VARTM result data, generating a seismic image of the subsurface formation; and
based on the seismic image, drilling into the subsurface formation.