System and method for seismic imaging
A method is described for seismic imaging including receiving a full waveform inversion (FWI) image of a subsurface volume of interest including complex geology; transforming the FWI image to the Fourier domain to generate a Fourier domain image; performing a multi-dimensional low-cut filter of the Fourier domain image to generate a low-cut filter image; phase rotating by 90 degrees the low-cut filter image to generate a phase-rotated image; performing an inverse Fourier transform on the phase-rotated image to generate a transformed image; and displaying the transformed image on a graphical display. The method is executed by a computer system.
1 . A computer-implemented method of seismic imaging, comprising:
a. receiving, at a computer processor, a full waveform inversion (FWI) image of a subsurface volume of interest including complex geology with two dips at subsurface points;
b. transforming, via the computer processor, the FWI image to the Fourier domain to generate a Fourier domain image, wherein the Fourier domain image is expressed in wavenumber space as Image FWI (k x , k y , k z );
c. performing, via the computer processor, a multi-dimensional low-cut filter of the Fourier domain image to generate a low-cut filter image, wherein the multi-dimensional low-cut filter is applied to k r , where
k
r
=
k
x
2
+
k
y
2
+
k
z
2
to
generate the low-cut filter image Low_cut FWI (k x , k y , k z ) and wherein there is no direction discrimination so all dips in the low-cut filter image are preserved;
d. phase rotating by 90 degrees the low-cut filter image, via the computer processor, to generate a phase-rotated image in order to match a phase of a seismic migration image;
e. performing an inverse Fourier transform on the phase-rotated image, via the computer processor, to generate a transformed image, wherein the transformed image includes conflicting dips of reflectors representative of the complex geology; and
f. using the transformed image for improved hydrocarbon exploration and hydrocarbon production by identifying traps for subsurface hydrocarbon deposits or hazards for drilling.
2 . The method of claim 1 wherein the phase rotating is performed as Rotated FWI (k x , k y , k z )=i*Low_cut FWI (k x , k y , k z ).
3 . A computer system, comprising:
one or more processors;
memory; and
one or more programs, wherein the one or more programs are stored in the memory and configured to be executed by the one or more processors, the one or more programs including instructions that when executed by the one or more processors cause the system to:
a. receive, at the one or more processors, a full waveform inversion (FWI) image of a subsurface volume of interest including complex geology with two dips at subsurface points;
b. transform the FWI image to the Fourier domain to generate a Fourier domain image, wherein the Fourier domain image is expressed in wavenumber space as Image FWI (k x , k y , k z );
c. perform a multi-dimensional low-cut filter of the Fourier domain image to generate a low-cut filter image, wherein the multi-dimensional low-cut filter is applied to
k
r
=
k
x
2
+
k
y
2
+
k
z
2
to generate the low-cut filter image Low_cut FWI (k x , k y , k z ) and wherein there is no direction discrimination so all dips in the low-cut filter image are preserved;
d. phase rotate by 90 degrees the low-cut filter image to generate a phase-rotated image in order to match a phase of a seismic migration image;
e. perform an inverse Fourier transform on the phase-rotated image to generate a transformed image, wherein the transformed image includes conflicting dips of reflectors representative of the complex geology; and
f. using the transformed image for improved hydrocarbon exploration and hydrocarbon production by identifying traps for subsurface hydrocarbon deposits or hazards for drilling.
4 . The computer system of claim 3 wherein the instructions include that the phase rotating is performed as Rotated FWI (k x , k y , k z )=i*Low_cut FWI (k x , k y , k z ).
5 . A non-transitory computer readable storage medium storing one or more programs, the one or more programs comprising instructions, which when executed by an electronic device with one or more processors and memory, cause the device to
a. receive, at the one or more processors, a full waveform inversion (FWI) image of a subsurface volume of interest including complex geology with two dips at subsurface points;
b. transform the FWI image to the Fourier domain to generate a Fourier domain image, wherein the Fourier domain image is expressed in wavenumber space as Image FWI (k x , k y , k z );
c. perform a multi-dimensional low-cut filter of the Fourier domain image to generate a low-cut filter image, wherein the multi-dimensional low-cut filter is applied to k r , where
k
r
=
k
x
2
+
k
y
2
+
k
z
2
to generate the low-cut filter image Low_cut FWI (k x , k y , k z ) and wherein there is no direction discrimination so all dips in the low-cut filter image are preserved;
d. phase rotate by 90 degrees the low-cut filter image to generate a phase-rotated image in order to match a phase of a seismic migration image;
e. perform an inverse Fourier transform on the phase-rotated image to generate a transformed image, wherein the transformed image includes conflicting dips of reflectors representative of the complex geology; and
f. using the transformed image for improved hydrocarbon exploration and hydrocarbon production by identifying traps for subsurface hydrocarbon deposits or hazards for drilling.
6 . The electronic device of claim 5 wherein the instructions include that the phase rotating is performed as Rotated FWI (k x , k y , k z )=i*Low_cut FWI (k y , k y , k z ).