FWI model domain angle stacks with amplitude preservation
A method, including: obtaining a seismic dataset that is separated into subsets according to predetermined subsurface reflection angle ranges; performing, with a computer, an acoustic full wavefield inversion process on each of the subsets, respectively, to invert for density and generate respective density models; generating acoustic impedances for each of the subsets, as a function of reflection angle, using the respective density models; and transforming, using a computer, the acoustic impedances for each of the subsets into reflectivity sections, wherein the transforming includes normalizing the reflectivity sections by their respective bandwidth.
1. A method, comprising:
obtaining a seismic dataset that is separated into subsets according to predetermined subsurface reflection angle ranges;
performing, with a computer, an acoustic full wavefield inversion process on each of the subsets, respectively, to invert for density and generate respective density models;
generating acoustic impedances for each of the subsets, as a function of reflection angle, using the respective density models;
transforming, using a computer, the acoustic impedances for each of the subsets into reflectivity sections, wherein the transforming includes normalizing the reflectivity sections by their respective bandwidth; and
using, for each of the reflectivity sections, a Fourier transform, discrete Fourier transform, or a fast Fourier transform to calculate an average spectrum within at least one local window that is applied at a same location to all of the reflectivity sections, and determining a bandwidth for each average spectrum.
2. The method of claim 1 , wherein each of the full wavefield inversion processes start from a same velocity model.
3. The method of claim 1 , wherein each of full wavefield inversion processes are independently applied to the subsets.
4. The method of claim 1 , wherein the obtaining includes dividing a shot gather into the subsets by using a data mask that includes information of reflector dipping angles and P-wave velocity.
5. The method of claim 1 , wherein the determining the bandwidth is based on a distance between 10-dB points.
6. The method of claim 1 , wherein the determining the bandwidth is based on a distance between points with steepest slope.
7. The method of claim 1 , wherein the average spectrum is calculated within a plurality of local windows, and is averaged.
8. The method of claim 1 , further comprising determining reflectivity values at a plurality of angles and constructing an angle-vs-amplitude curve by interpolation.
9. The method of claim 1 , further comprising managing hydrocarbon production using the reflectivity sections.
10. The method of claim 1 , wherein the managing hydrocarbon production includes drilling a well at a location determined at least in part by the reflectivity sections.
11. A non-transitory computer readable storage medium encoded with instructions, which when executed by a computer cause the computer to implement a method comprising:
obtaining a seismic dataset that is separated into subsets according to predetermined subsurface reflection angle ranges;
performing, with a computer, an acoustic full wavefield inversion process on each of the subsets, respectively, to invert for density and generate respective density models;
generating acoustic impedances for each of the subsets, as a function of reflection angle, using the respective density models;
transforming, using a computer, the acoustic impedances for each of the subsets into reflectivity sections, wherein the transforming includes normalizing the reflectivity sections by their respective bandwidth; and
using, for each of the reflectivity sections, a Fourier transform, discrete Fourier transform, or a fast Fourier transform to calculate an average spectrum within at least one local window that is applied at a same location to all of the reflectivity sections, and determining a bandwidth for each average spectrum.