Difference between revisions of "F-k filtering"

From SEG Wiki
Jump to: navigation, search
(added noise removal example image)
(added another external link and fixed some grammar issues)
(2 intermediate revisions by the same user not shown)
Line 1: Line 1:
[https://wiki.seg.org/wiki/Dictionary:F-k_filtering F-k filtering] is when seismic data, that is traditionally in the time and displacement domain, is converted into the frequency and wave number (F-k) domain as seen in figure 1 below.  The seismic data is then filtered to remove certain values, and then converted back to the time displacement domain.  The process is a two dimensional Fourier transformation and must be sampled according to the  [[Nyquist criterion]].<ref>“BASIC DEFINITIONS.” ''Basic Definitions'', www.xsgeo.com/course/basic.htm.</ref>
+
[https://wiki.seg.org/wiki/Dictionary:F-k_filtering F-k filtering] is when seismic data, that is traditionally in the time and displacement domain, is converted into the frequency and wave number (F-k) domain as seen in [[:File:Fk domain.png|figure 1]] below.  The [[Seismic data acquisition|seismic data]] is then filtered to remove unwanted frequencies higher and/or lower than the seismic signal band, and then converted back to the time-displacement domain.  The process is a two dimensional Fourier transformation and must be sampled according to the  Nyquist criterion to avoid aliasing.<ref>“BASIC DEFINITIONS.” ''Basic Definitions'', www.xsgeo.com/course/basic.htm.</ref>  
 
[[File:Fk domain.png|none|frame|Figure 1: Surface and reflected waves in the time–space domain before and after phase-shift processing in the FK<sub><font size="2">X</font></sub>K<sub><font size="2">Y</font></sub> domain: (a) surface and reflected waves in the time–space domain; (b) phase of the surface and reflected waves in the FK<sub><font size="2">X</font></sub>K<sub><font size="2">Y</font></sub> domain; (c) phase of the surface and reflected waves after phase-shift processing in the FK<sub><font size="2">X</font></sub>K<sub><font size="2">Y</font></sub> domain<ref>Wang, D., & Ling, Y. (2016). Phase-shift- and phase-filtering-based surface-wave suppression method. ''Applied Geophysics,'' ''13''(4), 614-620.</ref>]]
 
[[File:Fk domain.png|none|frame|Figure 1: Surface and reflected waves in the time–space domain before and after phase-shift processing in the FK<sub><font size="2">X</font></sub>K<sub><font size="2">Y</font></sub> domain: (a) surface and reflected waves in the time–space domain; (b) phase of the surface and reflected waves in the FK<sub><font size="2">X</font></sub>K<sub><font size="2">Y</font></sub> domain; (c) phase of the surface and reflected waves after phase-shift processing in the FK<sub><font size="2">X</font></sub>K<sub><font size="2">Y</font></sub> domain<ref>Wang, D., & Ling, Y. (2016). Phase-shift- and phase-filtering-based surface-wave suppression method. ''Applied Geophysics,'' ''13''(4), 614-620.</ref>]]
  
 
== Background ==
 
== Background ==
[[Dictionary:Acoustic|Acoustic]] signals that are not reflections appear in shot records as [[Dictionary:Noise (seismic)|noise]]<ref name=":1">Seismic processing basics. (n.d.). Retrieved from <nowiki>https://wiki.aapg.org/Seismic_processing_basics#f.E2.80.93k_or_apparent_velocity_filter</nowiki>.</ref>. These signals have a constant “apparent velocity” as they travel along the receiver cable. This simple organization allows them to be isolated from the reflection signal and to be removed from the record.  To accomplish this you can use the F-k, or pie slice, filter in which the range of apparent velocities can be selected and removed to eliminate this linear [[Dictionary:Noise (seismic)|noise]].
+
[[Dictionary:Acoustic|Acoustic]] signals that are not reflections appear in shot records as [[Dictionary:Noise (seismic)|noise]]<ref name=":1">Seismic processing basics. (n.d.). Retrieved from <nowiki>https://wiki.aapg.org/Seismic_processing_basics#f.E2.80.93k_or_apparent_velocity_filter</nowiki>.</ref>. These signals have a constant “apparent velocity” as they travel along the receiver cable. This simple organization allows them to be isolated from the reflection signal and to be removed from the record.  To accomplish this you can use the F-k, or pie slice, filter in which the range of apparent velocities can be selected and removed to eliminate this linear noise.
 
== Process ==
 
== Process ==
 
[[File:Eq.png|Figure 2: FFT<sub>xy</sub>( ) denotes the two-dimensional Fourier transform with x, y as variables, e<sup>iθd( f, kx, ky)</sup> , e<sup>iθs( f, kx, ky)</sup> , and e<sup>iθg( f, kx, ky)</sup> denote the two-dimensional Fourier transform of e<sup>iθd( f, kx, ky)</sup> ,e<sup>iθs( f, kx, ky)</sup> and e<sup>iθg( f, kx, ky)</sup> with x, y as variables<ref>Wang, D., & Ling, Y. (2016). Phase-shift- and phase-filtering-based surface-wave suppression method. ''Applied Geophysics,'' ''13''(4), 614-620.
 
[[File:Eq.png|Figure 2: FFT<sub>xy</sub>( ) denotes the two-dimensional Fourier transform with x, y as variables, e<sup>iθd( f, kx, ky)</sup> , e<sup>iθs( f, kx, ky)</sup> , and e<sup>iθg( f, kx, ky)</sup> denote the two-dimensional Fourier transform of e<sup>iθd( f, kx, ky)</sup> ,e<sup>iθs( f, kx, ky)</sup> and e<sup>iθg( f, kx, ky)</sup> with x, y as variables<ref>Wang, D., & Ling, Y. (2016). Phase-shift- and phase-filtering-based surface-wave suppression method. ''Applied Geophysics,'' ''13''(4), 614-620.
 
</ref>. |frame]]
 
</ref>. |frame]]
Without considering random noise the seismic record is noted as d(t,x,y) = s(t,x,y) + g(t,x,y)  where d is the seismic record and is a sum of the reflected wave (s) and the surface wave (g).  In the F-k domain this equation becomes D(ƒ,x,y)e<sup>iΘd(ƒ,x,y)</sup> = S(ƒ,x,y)e<sup>iΘs(ƒ,x,y)</sup> + G(ƒ,x,y)e<sup>iΘg(ƒ,x,y)</sup>.  To obtain the results in the F-k domain a two-dimensional [[wikipedia:Fourier_transform|Fourier transformation]] is done along the x–y plane using the equation in the figure 2.<ref name=":0" />
+
Without considering random noise the seismic record is noted as d(t,x,y) = s(t,x,y) + g(t,x,y)  where d is the seismic record and is a sum of the reflected wave (s) and the surface wave (g).  In the F-k domain this equation becomes D(ƒ,x,y)e<sup>iΘd(ƒ,x,y)</sup> = S(ƒ,x,y)e<sup>iΘs(ƒ,x,y)</sup> + G(ƒ,x,y)e<sup>iΘg(ƒ,x,y)</sup>.  To obtain the results in the F-k domain a two-dimensional [[wikipedia:Fourier_transform|Fourier transformation]] is done along the x–y plane using the equation in [[:File:Eq.png|figure 2]].<ref name=":0" />
  
 
According to the [[wikipedia:Dispersion_relation|dispersion relation]] the equation ƒ = vk must be satisfied<ref name=":0" />
 
According to the [[wikipedia:Dispersion_relation|dispersion relation]] the equation ƒ = vk must be satisfied<ref name=":0" />
Line 14: Line 14:
 
[[File:Seismograms-of-shot-19-a-raw-data-b-after-mute-and-trace-edition-F-k-filter-After.png|frame|Figure 3: Seismograms of shot 19: (a) raw data; (b) after mute and trace edition. F-k filter After editing and mute of trace amplitudes, commonly in seismic processing flow chart, it recommends some kind of filtering. The first filter used was f-k filter, which is a two- dimensional filter defined in the frequency domain obtained by Double Fourier Transform: <ref>Comparison of FK and SVD filtering in the processing of a land seismic data - Scientific Figure on ResearchGate. Available from: <nowiki>https://www.researchgate.net/figure/Seismograms-of-shot-19-a-raw-data-b-after-mute-and-trace-edition-F-k-filter-After_fig1_280521231</nowiki> [accessed 4 Dec, 2019]</ref>]]
 
[[File:Seismograms-of-shot-19-a-raw-data-b-after-mute-and-trace-edition-F-k-filter-After.png|frame|Figure 3: Seismograms of shot 19: (a) raw data; (b) after mute and trace edition. F-k filter After editing and mute of trace amplitudes, commonly in seismic processing flow chart, it recommends some kind of filtering. The first filter used was f-k filter, which is a two- dimensional filter defined in the frequency domain obtained by Double Fourier Transform: <ref>Comparison of FK and SVD filtering in the processing of a land seismic data - Scientific Figure on ResearchGate. Available from: <nowiki>https://www.researchgate.net/figure/Seismograms-of-shot-19-a-raw-data-b-after-mute-and-trace-edition-F-k-filter-After_fig1_280521231</nowiki> [accessed 4 Dec, 2019]</ref>]]
 
F-k filtering is used for things like eliminating the aliasing [[Dictionary:Surface wave|surface-wave]] energy and maintain the low frequency information of the reflected waves,<ref name=":0">Wang, D., & Ling, Y. (2016). Phase-shift- and phase-filtering-based surface-wave suppression method. ''Applied Geophysics,'' ''13''(4), 614-620.
 
F-k filtering is used for things like eliminating the aliasing [[Dictionary:Surface wave|surface-wave]] energy and maintain the low frequency information of the reflected waves,<ref name=":0">Wang, D., & Ling, Y. (2016). Phase-shift- and phase-filtering-based surface-wave suppression method. ''Applied Geophysics,'' ''13''(4), 614-620.
</ref> and attenuating the residual weak energy of ground roll<ref>Chen, H., Li, X., Qian, Z., & Zhao, G. (2013). Robust adaptive polarization analysis method for eliminating ground roll in 3C land seismics. ''Applied Geophysics,'' ''10''(3), 295-304.</ref>.  The main application for F-k filtering is to eliminate coherent noise in seismic data as seen in figure 3.<ref>Herman, M., Syahmi Hashim, H., Abdul Latif, A.H., Ghosh, D.P., 2017. Application of FK Filtering for Coherent Noise Removal in High Frequency Shallow Marine Data. IOP Conference Series: Earth and Environmental Science 88, 012010.</ref>
+
</ref> and [[Attenuation of air waves|attenuating]] the residual weak energy of [[ground roll]]<ref>Chen, H., Li, X., Qian, Z., & Zhao, G. (2013). Robust adaptive polarization analysis method for eliminating ground roll in 3C land seismics. ''Applied Geophysics,'' ''10''(3), 295-304.</ref>.  The main application for F-k filtering is to eliminate coherent noise in seismic data as exemplified by [[:File:Seismograms-of-shot-19-a-raw-data-b-after-mute-and-trace-edition-F-k-filter-After.png|figure 3]].<ref>Herman, M., Syahmi Hashim, H., Abdul Latif, A.H., Ghosh, D.P., 2017. Application of FK Filtering for Coherent Noise Removal in High Frequency Shallow Marine Data. IOP Conference Series: Earth and Environmental Science 88, 012010.</ref>
  
 
== Issues ==
 
== Issues ==
 
Traditional F-k filtering cannot eliminate all of [[spatial aliasing]] of surface waves. <ref name=":0" />
 
Traditional F-k filtering cannot eliminate all of [[spatial aliasing]] of surface waves. <ref name=":0" />
  
A filter that covers too wide of a range can remove too great an amount of information and make interpretation correctly difficult.<ref name=":1" />
+
A filter that covers too wide of a range can remove too great an amount of information and make interpreting correctly difficult.<ref name=":1" />
  
 
== External Links ==
 
== External Links ==
 
Domain Transformations - http://www.xsgeo.com/course/basic.htm
 
Domain Transformations - http://www.xsgeo.com/course/basic.htm
 +
 +
Filtering Techniques - https://gpg.geosci.xyz/content/seismic/seismic_reflection_filtering.html
  
 
Fourier Transformation - https://en.wikipedia.org/wiki/Fourier_transform
 
Fourier Transformation - https://en.wikipedia.org/wiki/Fourier_transform

Revision as of 21:32, 3 December 2019

F-k filtering is when seismic data, that is traditionally in the time and displacement domain, is converted into the frequency and wave number (F-k) domain as seen in figure 1 below. The seismic data is then filtered to remove unwanted frequencies higher and/or lower than the seismic signal band, and then converted back to the time-displacement domain. The process is a two dimensional Fourier transformation and must be sampled according to the Nyquist criterion to avoid aliasing.[1]

Figure 1: Surface and reflected waves in the time–space domain before and after phase-shift processing in the FKXKY domain: (a) surface and reflected waves in the time–space domain; (b) phase of the surface and reflected waves in the FKXKY domain; (c) phase of the surface and reflected waves after phase-shift processing in the FKXKY domain[2]

Background

Acoustic signals that are not reflections appear in shot records as noise[3]. These signals have a constant “apparent velocity” as they travel along the receiver cable. This simple organization allows them to be isolated from the reflection signal and to be removed from the record. To accomplish this you can use the F-k, or pie slice, filter in which the range of apparent velocities can be selected and removed to eliminate this linear noise.

Process

Figure 2: FFTxy( ) denotes the two-dimensional Fourier transform with x, y as variables, eiθd( f, kx, ky) , eiθs( f, kx, ky) , and eiθg( f, kx, ky) denote the two-dimensional Fourier transform of eiθd( f, kx, ky) ,eiθs( f, kx, ky) and eiθg( f, kx, ky) with x, y as variables[4].

Without considering random noise the seismic record is noted as d(t,x,y) = s(t,x,y) + g(t,x,y) where d is the seismic record and is a sum of the reflected wave (s) and the surface wave (g). In the F-k domain this equation becomes D(ƒ,x,y)eiΘd(ƒ,x,y) = S(ƒ,x,y)eiΘs(ƒ,x,y) + G(ƒ,x,y)eiΘg(ƒ,x,y). To obtain the results in the F-k domain a two-dimensional Fourier transformation is done along the x–y plane using the equation in figure 2.[5]

According to the dispersion relation the equation ƒ = vk must be satisfied[5]

Applications

Figure 3: Seismograms of shot 19: (a) raw data; (b) after mute and trace edition. F-k filter After editing and mute of trace amplitudes, commonly in seismic processing flow chart, it recommends some kind of filtering. The first filter used was f-k filter, which is a two- dimensional filter defined in the frequency domain obtained by Double Fourier Transform: [6]

F-k filtering is used for things like eliminating the aliasing surface-wave energy and maintain the low frequency information of the reflected waves,[5] and attenuating the residual weak energy of ground roll[7]. The main application for F-k filtering is to eliminate coherent noise in seismic data as exemplified by figure 3.[8]

Issues

Traditional F-k filtering cannot eliminate all of spatial aliasing of surface waves. [5]

A filter that covers too wide of a range can remove too great an amount of information and make interpreting correctly difficult.[3]

External Links

Domain Transformations - http://www.xsgeo.com/course/basic.htm

Filtering Techniques - https://gpg.geosci.xyz/content/seismic/seismic_reflection_filtering.html

Fourier Transformation - https://en.wikipedia.org/wiki/Fourier_transform

Dispersion Relationship - https://en.wikipedia.org/wiki/Dispersion_relation

Sources

  1. “BASIC DEFINITIONS.” Basic Definitions, www.xsgeo.com/course/basic.htm.
  2. Wang, D., & Ling, Y. (2016). Phase-shift- and phase-filtering-based surface-wave suppression method. Applied Geophysics, 13(4), 614-620.
  3. 3.0 3.1 Seismic processing basics. (n.d.). Retrieved from https://wiki.aapg.org/Seismic_processing_basics#f.E2.80.93k_or_apparent_velocity_filter.
  4. Wang, D., & Ling, Y. (2016). Phase-shift- and phase-filtering-based surface-wave suppression method. Applied Geophysics, 13(4), 614-620.
  5. 5.0 5.1 5.2 5.3 Wang, D., & Ling, Y. (2016). Phase-shift- and phase-filtering-based surface-wave suppression method. Applied Geophysics, 13(4), 614-620.
  6. Comparison of FK and SVD filtering in the processing of a land seismic data - Scientific Figure on ResearchGate. Available from: https://www.researchgate.net/figure/Seismograms-of-shot-19-a-raw-data-b-after-mute-and-trace-edition-F-k-filter-After_fig1_280521231 [accessed 4 Dec, 2019]
  7. Chen, H., Li, X., Qian, Z., & Zhao, G. (2013). Robust adaptive polarization analysis method for eliminating ground roll in 3C land seismics. Applied Geophysics, 10(3), 295-304.
  8. Herman, M., Syahmi Hashim, H., Abdul Latif, A.H., Ghosh, D.P., 2017. Application of FK Filtering for Coherent Noise Removal in High Frequency Shallow Marine Data. IOP Conference Series: Earth and Environmental Science 88, 012010.