System and method for interpolating seismic data
AssigneeCgg Services SAS
InventorStewart Trickett
Life of this patent
Estimating the term…
Abstract
A system and method of interpolating seismic data is provided. The system and method form a plurality of pairwise Hankel tensors from acquired seismic data, and a respective pairwise Hankel tensor for each of a plurality of originally collected frequency slices, perform tensor completion on each of said pairwise Hankel tensors to recover a plurality of interpolated frequency slices, and combine said plurality of interpolated frequency slices with said originally collected frequency slices to form a set of trace data of a geographical area of interest.
Claims20
A seismic survey method that improves an image of an explored subsurface geological formation by mitigating missing or faulty data, the method comprising:
- generating seismic waves with seismic sources;
- recording, with seismic receivers, seismic data upon detecting the seismic waves including ones traveling through the explored subsurface geological formation;
- determining that one or more traces of the seismic data are missing, corrupt or noisy;
- forming, with a computing device, a plurality of pairwise Hankel tensors from the acquired seismic data, so that a pairwise Hankel tensor corresponds to each of a plurality of originally collected frequency slices and has two orders for each spatial dimension of the originally collected frequency slices, a number of spatial dimensions being maxim four;
- performing a tensor completion for at least one of the pairwise Hankel tensors, to obtain interpolated seismic data suitable to replace the one or more traces of the seismic data that are missing, corrupt or noisy;
- combining the seismic data with the interpolated seismic data to obtain a complete set of traces;
- and
- generating the image based on the complete set of traces, to estimate an oil and gas yielding potential of the explored subsurface geological formation.
The method of 1, wherein said performing of tensor completion comprises:
- determining a low rank tensor R that minimizes ∥Z(TH−R)∥F, which is a Frobenius norm of said pairwise Hankel tensor TH, wherein ∥·∥F is the Frobenius norm of the tensor elements, and Z(·) is an operator that zeroes out all elements that are unknown in TH.
A seismic survey system that improves an image of an explored subsurface geological formation by mitigating missing or faulty data, comprising:
- a receiver configured to acquire transmitted seismic data detecting seismic waves with seismic sources after the seismic waves travelled through the explored subsurface geological formation;
- and
- a processor configured to,
- determine that one or more traces of the seismic data are missing, corrupt or noisy;
- form a plurality of pairwise Hankel tensors from said acquired seismic data, a respective pairwise Hankel tensor being formed for each of a plurality of originally collected frequency slices and each of the pairwise Hankel tensors having two orders for each spatial dimension of the originally collected frequency slices, a number of spatial dimensions being maxim four,
- perform tensor completion at least one of said pairwise Hankel tensors to obtain interpolated seismic data suitable to replace the one or more traces of the seismic data that are missing, corrupt or noisy,
- combine the seismic data with the interpolated seismic data to obtain a complete set of traces, and
- generate the image of the explored subsurface geological formation based on the complete set of traces, for estimating an oil and gas yielding potential of the explored subsurface geological formation.
Dependencies are inferred from the claim text.