US 11,852,767 B2
Robust full waveform inversion of seismic data method and device
Bingbing Sun, Thuwal (SA); and Tariq Alkhalifah, Thuwal (SA)
Assigned to KING ABDULLAH UNIVERSITY OF SCIENCE AND TECHNOLOGY, Thuwal (SA)
Appl. No. 16/981,058
Filed by KING ABDULLAH UNIVERSITY OF SCIENCE AND TECHNOLOGY, Thuwal (SA)
PCT Filed Jan. 25, 2019, PCT No. PCT/IB2019/050634
§ 371(c)(1), (2) Date Sep. 15, 2020,
PCT Pub. No. WO2019/186286, PCT Pub. Date Oct. 3, 2019.
Claims priority of provisional application 62/713,246, filed on Aug. 1, 2018.
Claims priority of provisional application 62/648,569, filed on Mar. 27, 2018.
Prior Publication US 2021/0018639 A1, Jan. 21, 2021
Int. Cl. G01V 1/30 (2006.01)
CPC G01V 1/303 (2013.01) [G01V 2210/614 (2013.01); G01V 2210/6161 (2013.01); G01V 2210/6222 (2013.01)] 11 Claims
OG exemplary drawing
 
1. A method for calculating a velocity model for a subsurface of the earth, the method comprising:
receiving measured seismic data d;
calculating predicted seismic data p;
selecting a matching filter w such that a convolution operation between the matching filter w and one of the measured seismic data d or the predicted seismic data p is equal to the other one of the measured seismic data d or the predicted seismic data p;
selecting a misfit function J that calculates (1) a distance between the matching filter w and a Dirac Delta function or (2) a travel time shift associated with the measured seismic data;
calculating a new velocity model using the misfit function J, the measured seismic data d, and the predicted seismic data p; and
generating an image of the subsurface of the earth, based on the new velocity model, for identifying resources in the subsurface,
wherein the measured seismic data d includes wavefields generated by a seismic source and the wavefields propagate through the subsurface where they are attenuated and reflected, and the attenuated and reflected wavefields are recorded by plural seismic receivers.