Next Article in Journal
Molecular Detection of Acarapis woodi Using Hive Debris as Innovative and Non-Invasive Matrix
Next Article in Special Issue
Stability Analysis and Fracture Mechanism of Surrounding Rock in a Hydraulic Tunnel Based on Microseismic Monitoring
Previous Article in Journal
Multi-Modal Long-Term Person Re-Identification Using Physical Soft Bio-Metrics and Body Figure
Previous Article in Special Issue
Influence of Super Absorbent Polymer on Root Characteristics and Anchorage of Amorpha fruticosa on Rocky Slope
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Lp-Order Time–Frequency Spectra and Its Application in Edge Detection

1
School of Geophysics, Chengdu University of Technology, Chengdu 610059, China
2
State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Chengdu 610059, China
3
School of Geoscience and Technology, Southwest Petroleum University, Chengdu 610500, China
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2022, 12(6), 2836; https://doi.org/10.3390/app12062836
Submission received: 24 December 2021 / Revised: 21 February 2022 / Accepted: 28 February 2022 / Published: 10 March 2022
(This article belongs to the Special Issue Structural Geology, Rock Mechanics and Their On-Site Testing Analysis)

Abstract

:
Generalized Hilbert transform (GHT) has been widely and successfully used in signal processing and interpretation for many years. However, the order parameter in GHT can only be chosen in a certain range, and the GHT result is not stable with respect to the chosen order parameter, which usually leads to poor-stability and low-resolution outputs. To improve the stability of GHT and broaden its order parameter chosen range, and finally achieve high-efficiency and high-resolution results, we extend GHT by introducing the sign function, the two frequency integration parameters, and propose the Lp-order time–frequency spectra (LTFS). In LTFS, the sign function aims to broaden the chosen range of the order parameter into any positive numbers, and then leads to high-stability and high-resolution LTFS results; the two frequency integration parameters aim to optimize the computation efficiency and output resolution of LTFS. Synthetic one-dimensional (1D) channel data testing, synthetic two-dimensional (2D) seismic data, and actual three-dimensional (3D) seismic data examples demonstrate favorable capability of LTFS for edge detection.

1. Introduction

The generalized Hilbert transform (GHT) [1] improves the noise robustness and resolution of the Hilbert transform (HT) [2], by introducing the order parameter (a real number) to the time–frequency spectra, which are obtained by time–frequency analysis methods such as the short-time Fourier transform (STFT) [3], wavelet transform [4,5,6], S transform [7,8], generalized S transform [9], EMD [10,11,12], VMD [13], modified short-time Fourier transform [14,15,16], synchrosqueezing transform [17,18,19,20], and so on.
Due to the good performance of noise robustness and resolution, GHT has been widely used in signal processing and interpretation. However, GHT results depend on the time–frequency spectra and the order parameter; and for a chosen order parameter, the output resolution of GHT is unique, which limits its application. In addition, the order parameter in GHT can only be chosen in a certain range; and the GHT result is not stable with respect to the chosen order parameter, which usually leads to poor-stability and low-resolution output.
In order to improve the stability of GHT and broaden its order parameter chosen range, thus obtain high-stability and high-resolution results; we first extend GHT by introducing the sign function, and then propose the high-stability and high-resolution GHT (SGHT). In SGHT, due to the sign function, the order parameter can be chosen as any positive numbers. Thus, we can obtain high-stability and high-resolution GHT results. In addition, the frequency integral interval in GHT and SGHT is fixed, which limits the corresponding computation efficiency and output resolution. In order to improve the computation efficiency and output resolution of GHT and SGHT, and finally obtain high-efficiency and high-resolution results, we extend SGHT by introducing the two frequency integration parameters, and then propose the Lp-order time–frequency spectra (LTFS). In LTFS, we can flexibly choose the order parameter and the two frequency integration parameters, and finally optimize the computation efficiency and output resolution of LTFS for more general application.
The remainder of this paper is organized as follows. First, we introduce the principle of GHT and LTFS. Then, we use synthetic two-dimensional (2D) seismic data to test the capability of LTFS for edge detection, with the comparison of HT and GHT. Finally, we apply LTFS to actual three-dimensional (3D) seismic data, with comparison of HT, GHT, curvature [21,22], and coherence [22,23].

2. Generalized Hilbert Transform (GHT)

GHT extends HT (Appendix A) by introducing the order parameter in the time–frequency domain to enhance the noise robustness and resolution of HT, which make GHT easier for new applications than HT.

2.1. Theory of GHT

If complex numbers S ( t , f ) denote the time–frequency spectra of the inputted real signal x ( t ) , which are obtained by the short-time Fourier transform (STFT) [3], wavelet transform [4,5,6], S transform [7,8], generalized S transform [9], EMD [10,11,12], VMD [13], modified short-time Fourier transform [14,15,16], synchrosqueezing transform [17,18,19,20], and so on; then, S ( t , f ) can be modified as
R ( t , f ) = 2 S ( t , f ) S ( t , f ) f > 0 f = 0 0 f < 0
and further considering two real numbers S 1 ( t , f ) and S 2 ( t , f ) as the real and imaginary parts of R ( t , f ) , respectively, then, we have
S 1 ( t , f ) = Re R ( t , f ) S 2 ( t , f ) = Im R ( t , f )
where Re   and Im   takes the real and imaginary parts of a complex number, respectively.
Thus, the GHT result of x ( t ) can be expressed as
g ( t ) = g 3 ( t ) + i g 4 ( t )
and
g 1 ( t ) = 0 f N S 1 ( t , f ) p d f g 2 ( t ) = 0 f N S 2 ( t , f ) p d f g 3 ( t ) = g 1 ( t ) p 1 g 4 ( t ) = g 2 ( t ) p 1
where f N is the Nyquist frequency [24,25,26]; p is the order parameter and p > 0 .
Therefore, the instantaneous amplitudes, phases, and frequencies of the GHT result g ( t ) can be obtained by
tan θ 1 ( t ) = g 3 1 ( t ) g 4 ( t ) f 1 ( t ) = 1 2 π θ 1 ( t ) t A 1 ( t ) = g 3 2 ( t ) + g 4 2 ( t )
and A 1 ( t ) , θ 1 ( t ) , and f 1 ( t ) are referred to as the instantaneous amplitudes, phases, and frequencies of x ( t ) obtained by GHT, respectively.

2.2. Analysis of GHT

According to Equations (1)–(4), the GHT result g ( t ) depends on the time–frequency spectra S ( t , f ) and the order parameter p . For a chosen order parameter, the output resolution of GHT is unique, which limits its application. As S ( t , f ) is related to the time–frequency analysis methods [3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20], here we only analyze the impact of the chosen order parameter p on the GHT result g ( t ) .
If c denotes the collection of positive numbers, then it can be expressed as
c = 2 k 1 1 2 k 2 1 , 2 k 1 2 k 2 1 1 , 2 k 1 1 2 k 2 1 1
where { } takes the number collection; k 1 and k 2 are positive integer numbers; if A and B both are real numbers, and
B = A c = A ( 2 k 1 1 ) ( 2 k 2 ) A ( 2 k 1 ) ( 2 k 2 1 ) A ( 2 k 1 1 ) ( 2 k 2 1 )
then, the following four conclusions are drawn:
(i) For A < 0 , then B = A ( 2 k 1 1 ) ( 2 k 2 ) does not exist. Therefore, when S n ( t , f ) < 0 or g n ( t ) < 0 ( n = 1 , 2 ), p cannot be chosen as 2 k 1 1 2 k 2 1 .
(ii) For A < 0 , then B = A ( 2 k 1 ) ( 2 k 2 1 ) > 0 . Therefore, when S n ( t , f ) < 0 or g n ( t ) < 0 and p is chosen as 2 k 1 2 k 2 1 1 , the signs of S n ( t , f ) and g n ( t ) are changed and GHT obtains instability result (as shown in Figure 1d).
(iii) For A < 0 , then B = A ( 2 k 1 1 ) ( 2 k 2 1 ) < 0 . Therefore, when S n ( t , f ) < 0 and g n ( t ) < 0 , p can be chosen as 2 k 1 1 2 k 2 1 1 ; and correspondingly, GHT preserves the signs of S n ( t , f ) and g n ( t ) .
(iv) For A 0 , then B = A c 0 . Therefore, when S n ( t , f ) 0 and g n ( t ) 0 , p can be chosen as any number; and correspondingly, GHT preserves the signs of S n ( t , f ) and g n ( t ) .
Above all, the order parameter p can be only chosen as 2 k 1 1 2 k 2 1 1 and 2 k 1 2 k 2 1 1 ; the GHT result with p = 2 k 1 2 k 2 1 1 usually has a serious instability problem, which limits its further application.

2.3. Synthetic 1D Channel Data Testing

Here, we used one set of synthetic 1D channel data (as shown in Figure 1a) to analyze the GHT, with comparison of HT. The time–sample interval of the synthetic 1D channel data is one millisecond, thus, the Nyquist frequency fN is 500 Hz.
Figure 1b shows the instantaneous amplitudes obtained by HT, which indicate the edges of the synthetic channel; however, the HT result shows lower resolution than we desired. Figure 1c shows the instantaneous amplitudes obtained by GHT with p = 1.0. Figure 1d shows the instantaneous amplitudes obtained by GHT with p = 0.8 (ratio of an even number 4 and an odd number 5). Figure 1e shows the instantaneous amplitudes obtained by GHT p = 0.6 (ratio of two odd numbers 3 and 5).
We can see that when p = 1.0 and 0.6, GHT results effectively detected the edges of the synthetic channel and have a higher resolution than the HT result; whereas, GHT obtains a lower resolution or stability result with p = 0.8, due to the corresponding sign changing. Comparison of Figure 1c–e shows that GHT results are unstable with changing order parameter p.

3. Lp-Order Time–Frequency Spectra (LTFS)

LTFS aims to extend GHT by introducing the sign function, and two frequency integration parameters to produce high-stability and high-resolution results with high efficiency.

3.1. Theory of High-Stability GHT (SGHT)

In order to improve the stability of GHT and broaden its order parameter chosen range, and finally achieve high-stability and high-resolution GHT results, we can extend Equation (4) by introducing the sign function and propose the SGHT given by
s ( t ) = s 3 ( t ) + i s 4 ( t )
and
s 1 ( t ) = 0 f N Sign S 1 ( t , f ) S 1 ( t , f ) p d f s 2 ( t ) = 0 f N Sign S 2 ( t , f ) S 2 ( t , f ) p d f s 3 ( t ) = Sign s 1 ( t ) s 1 ( t ) p 1 s 4 ( t ) = Sign s 1 ( t ) s 1 ( t ) p 1
where || takes the absolute value of a real number; Sign   denotes the sign function.
Thus, instantaneous amplitudes, phases, and frequencies of s ( t ) can be obtained by
tan θ 2 ( t ) = s 3 1 ( t ) s 4 ( t ) f 2 ( t ) = 1 2 π θ 2 ( t ) t A 2 ( t ) = s 3 2 ( t ) + s 4 2 ( t )
and A 2 ( t ) , θ 2 ( t ) , and f 2 ( t ) are referred to as the instantaneous amplitudes, phases, and frequencies of x ( t ) obtained by SGHT, respectively.
According to Equations (8) and (9), as
Sign Sign S n ( t , f ) S n ( t , f ) p = Sign S n ( t , f ) Sign Sign s n ( t ) s n ( t ) p = Sign s n ( t )
thus, the order parameter p in SGHT can be chosen as any positive real number, which means that
p = c = 2 k 1 1 2 k 2 1 , 2 k 1 2 k 2 1 1 , 2 k 1 1 2 k 2 1 1

3.2. Theory of LTFS

Considering the amplitude spectra dynamic range [16,27,28], the actual signal belongs to the frequency band limited signals; and it means that amplitude spectra of actual signal ranges only in a certain frequency range, not from zero to the Nyquist frequency. Thus, in order to improve the computation efficiency and output resolution of SGHT, and finally obtain high-efficiency and high-resolution results; we extend SGHT by introducing the two frequency integration parameters, and then propose the LTFS given by
h f 1 , f 2 ( t ) = h 3 ( t ) + i h 4 ( t )
and
h 1 ( t ) = f 1 f 2 Sign S 1 ( t , f ) S 1 ( t , f ) p d f h 2 ( t ) = f 1 f 2 Sign S 2 ( t , f ) S 2 ( t , f ) p d f h 3 ( t ) = Sign s 1 ( t ) s 1 ( t ) p 1 h 4 ( t ) = Sign s 1 ( t ) s 1 ( t ) p 1
where f 1 and f 2 both are the frequency integration parameters, and 0 f 1 f 2 f N ; for actual seismic data, usually f 1 = 0 Hz and f 2 120 Hz.
Correspondingly, instantaneous amplitudes, phases, and frequencies of h f 1 , f 2 ( t ) can be obtained by
tan θ 3 ( t ) = h 3 1 ( t ) h 4 ( t ) f 3 ( t ) = 1 2 π θ 3 ( t ) t A 3 ( t ) = h 3 2 ( t ) + h 4 2 ( t )
and A 3 ( t ) , θ 3 ( t ) , and f 3 ( t ) are referred to as the instantaneous amplitudes, phases, and frequencies of x ( t ) obtained by LTFS, respectively.

3.3. Analysis of LTFS

According to Equations (13) and (14), the output of LTFS depends on the time–frequency spectra, the sign function, the order parameter, and the two frequency integration parameters. A comparison of SGHT and LTFS shows that LTFS is a general representation of SGHT, and SGHT is the special case of LTFS. Thus, we can flexibly choose the order parameter and the two frequency integration parameters, to optimize the output resolution and computation efficiency of LTFS, respectively. Finally, high-stability and high-resolution results can be obtained with high efficiency.
If we compare LTFS results with R ( t , f ) , GHT, and SGHT results, then
h f 1 , f 2 ( t ) = R ( t , f ) , if f 1 = f 2 = f h f 1 , f 2 ( t ) = g ( t ) , if p = 2 k 1 1 2 k 2 1 1 , f 1 = 0 , f 2 = f N h f 1 , f 2 ( t ) = s ( t ) , if f 1 = 0 , f 2 = f N
and if we consider we consider T 1 , T 2 , and T 3 as the time costs of GHT, SGHT, and LTFS, respectively; then, according to Equations (4), (9), and (14), we have
T 2 = T 1 T 1 T 3 1 = f N f 2 f 1 1
Above all, SGHT has higher stability than GHT without computation efficiency decreasing, while LTFS has higher computation efficiency than SGHT without stability decreasing. Moreover, LTFS have higher stability and computation efficiency than GHT, and finally obtain high-stability and high-resolution results with high efficiency.

3.4. Synthetic 1D Channel Data Testing

Figure 2a shows the instantaneous amplitudes of the synthetic 1D channel data in Figure 1a, which are obtained by LTFS with p = 1.0, f1 = 0 Hz, and f2 = fN = 500 Hz. Figure 2b shows the instantaneous amplitudes of the synthetic 1D channel data in Figure 1a, which are obtained by LTFS with p = 0.8, f1 = 0 Hz, and f2 = fN = 500 Hz. Figure 2c shows the instantaneous amplitudes of the synthetic 1D channel data in Figure 1a, which are obtained by LTFS with p = 0.6, f1 = 0 Hz, and f2 = fN = 500 Hz. LTFS has successfully detected the abrupt and gradual changes no matter which order parameter is chosen, and the LTFS results become more and more sharp with decreasing order parameter p.
According to Figure 1 and Figure 2, we can see that GHT and LTFS results with p = 1.0 are the same, while GHT and LTFS results with p = 0.6 are the same; in addition, LTFS results with p = 1.0 have much higher resolutions than GHT results with p = 1.0, which benefits from the sign function used in LTFS. Comparison of Figure 2a–c shows that the LTFS resolutions are increased with decreasing order parameter p.

4. Examples

In the three examples shown in this paper, the time–frequency spectra are generated by the modified short-time Fourier transform [14] given by
S ( t , f ) = + z ( τ + t ) w ( τ ) exp ( i 2 π f τ ) d τ
and
w ( τ ) = a f b 2 π exp a 2 f 2 b 2 τ 2
where a > 0 and b 0 .

4.1. Synthetic 2D Seismic Data Example

For a comparative study of the LTFS with HT and GHT in seismic edge detection, one set of synthetic 2D seismic data (as shown in Figure 3a) with three horizontal reflections and one gradient reflection (the gradient is 1 millisecond per trace) is generated by a 20 Hz dominant frequency zero-phase Ricker wavelet. The time–sample interval is one millisecond, thus the Nyquist frequency is 500 Hz and the reflection coefficients of the four reflections are 1.0, −1.0, 1.0, and −1.0, respectively.
Figure 3b shows the instantaneous amplitudes obtained by HT; Figure 3c shows the instantaneous amplitudes obtained by GHT, with p = 0.4 (ratio of an even number 2 and an odd number 5). Figure 3d shows the instantaneous amplitudes obtained by LTFS, with p = 0.4, f1 = 0 Hz, and f2 = fN = 500 Hz. Figure 3e shows the instantaneous amplitudes obtained by GHT with p = 0.2 (ratio of two odd numbers 1 and 5), which is the same as the LTFS result with p = 0.2, f1 = 0 Hz, and f2 = fN = 500 Hz. Figure 3f shows the instantaneous amplitudes obtained by LTFS, with p = 0.2, f1 = 10 Hz, and f2 = 100 Hz.
For more clear comparison, we extracted the corresponding results in the yellow rectangles in Figure 3. Figure 4a–f shows the zoomed-in views of yellow rectangles in Figure 3a–f, respectively; and we can see that Figure 4a–f clearly detected the edges between the second horizontal layer and the gradient layer when trace numbers are greater than 40, 42, 34, 24, 21, and 21, respectively. In contrast, the GHT result in Figure 4c shows strong artifacts (or instabilities) due to the sign changing during GHT. A comparison of Figure 4b–f indicates that Figure 4e,f obtains higher-resolution results than Figure 4b–d and Figure 4e,f has no obvious differences.
Table 1 shows the time costs of GHT and LTFS results in Figure 3e,f. Figure 4f has higher computation efficiency than Figure 4e due to the chosen frequency integration parameters f1 = 10 Hz and f2 = 100 Hz.
Above all, LTFS have higher resolution and computation efficiency than GHT, and finally obtain high-efficiency and high-resolution results.

4.2. Actual 3D Seismic Data Example

In order to better use the actual 3D seismic data for geologic body detection by LTFS, a four-step procedure was followed. First, HT was applied to the input actual 3D seismic data in the time (or depth) direction to generate its analytical representation. Second, MSTFT was applied to the analytical representation of the input actual 3D seismic data in the time direction to generate the time–frequency spectra. Third, the 3D instantaneous amplitudes obtained by HT, GHT, and LTFS were calculated. Finally, time slices of the instantaneous amplitudes obtained by HT, GHT, and LTFS were extracted for comparison.
Figure 5a,b shows the seismic data of inline 160 and crossline 96 extracted from the input actual 3D seismic data, respectively. The time–sample interval is one millisecond, thus the Nyquist frequency is 500 Hz.
Figure 6a shows a time slice of 1.7 s extracted from the input actual 3D seismic data. Figure 6b shows the same time slice of instantaneous amplitudes obtained by HT. Figure 6c shows the same time slice of instantaneous amplitudes obtained by GHT with p = 6 (an even number). Figure 6d shows the same time slice of the instantaneous amplitudes obtained by LTFS with p = 6, f1 = 0 Hz, and f2 = 90 Hz. For comparison, the same time slices of curvature and third-generation coherence (C3) attributes used for characterizing the channels are shown as Figure 6e,f, respectively.
Figure 7a–f shows the time slices in Figure 6a–f, respectively, and with two interpreted channels (shown as the two yellow curves). Comparison of Figure 7b with Figure 7a,c–f indicates that HT shows the main features of the two channels, but lacks the details obtained by GHT, LTFS, most positive curvatures, and third-generation coherences. A comparison of Figure 7d with Figure 7a–c,e,f indicates that Figure 7d is consistent with the two interpreted channels, and exhibits higher resolution and more detail than Figure 7a–c,e,f.

5. Conclusions

In this paper, we propose the LTFS to extend GHT by introducing the sign function, the two frequency integration parameters. The sign function broadens the chosen range of the order parameter into any positive numbers and leads to high-stability and high-resolution results. The two frequency integration parameters not only control the computation efficiency of LTFS, but also control its output resolution. Thus, we can flexibly choose the order parameter and two frequency integration parameters to optimize the computation efficiency and output resolution of LTFS. Synthetic 1D channel data testing, synthetic 2D seismic data, and actual 3D seismic data examples reveal that LTFS has higher resolution than GHT and can be widely used in signal processing and interpretation.

Author Contributions

Writing—original draft preparation, H.L.; methodology, H.L. and G.Z.; writing—review and editing, J.C. and Z.H.; software, Y.L.; data processing, K.L. All authors have read and agreed to the published version of the manuscript.

Funding

The National Natural Science Foundation of China (Grant Nos. 41874168 and 41974160), the Sichuan Science and Technology Program for Distinguished Young Scholars (Grant No. 2019JDJQ0053).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Theory of HT

HT [2] is derived from the Fourier transform (FT) [24,25,26] and has been used in signal processing and interpretation for many years. A well-known application of HT is seismic complex-trace analysis using instantaneous amplitude, phase, and frequency.
For a real signal x ( t ) in the time domain, its FT results can be expressed as
X ( f ) = x ( t ) exp i 2 π f t d t
where t is the time; f is the frequency; i = 1 ; and the HT result (a real signal) of real signal x ( t ) can be expressed as
y ( t ) = Im Z ( f ) exp i 2 π f t d f
and
Z ( f ) = 2 X ( f ) X ( f ) f > 0 f = 0 0 f < 0
where Im [   ] takes the imaginary part of a complex number.
Correspondingly, the analytical signal (or analytical representation) of the real signal x ( t ) obtained by HT can be expressed as
z ( t ) = x ( t ) + i y ( t ) = Z ( f ) exp i 2 π f t d f
And if A ( t ) , θ ( t ) , and f ( t ) denote the instantaneous amplitudes, phases, and frequencies of x ( t ) or z ( t ) , respectively, then
tan θ ( t ) = x 1 ( t ) y ( t ) f ( t ) = 1 2 π θ ( t ) t A ( t ) = x 2 ( t ) + y 2 ( t )
and all can be used for edge detection, hydrocarbon detection, and other applications, but with lower noise robustness and resolution.

References

  1. Luo, Y.; Saleh, A.; Marhoon, M. Generalized Hilbert transform and its applications in Geophysics. Lead. Edge 2003, 22, 198–202. [Google Scholar] [CrossRef]
  2. Taner, M.T.; Koehler, F.; Sheriff, R.E. Complex seismic trace analysis. Geophysics 1979, 44, 1041–1063. [Google Scholar] [CrossRef]
  3. Gabor, D. Theory of communication. J. Inst. Electr. Eng. 1946, 93, 429–457. [Google Scholar] [CrossRef]
  4. Kronland-Martinet, R.; Morlet, J.; Grossmann, A. Analysis of sound patterns through wavelet transforms. Int. J. Pattern Recognit. Artif. Intell. 1987, 1, 273–302. [Google Scholar] [CrossRef]
  5. Mallat, S.; Hwang, W.L. Singularity detection and processing with wavelets. IEEE Trans. Inf. Theory 1992, 38, 617–643. [Google Scholar] [CrossRef]
  6. Daubechies, I.; Maes, S. A nonlinear squeezing of the continuous wavelet transform based on auditory nerve models. In Wavelets in Medicine and Biology; CRC Press: Boca Raton, FL, USA, 1996; pp. 527–546. [Google Scholar]
  7. Stockwell, R.G.; Mansinha, L.; Lowe, R.P. Localization of the complex spectrum: The S-transform. IEEE Trans. Signal Process. 1996, 44, 998–1001. [Google Scholar] [CrossRef]
  8. Pinnegar, C.R.; Mansinha, L. The S-transform with windows of arbitrary and varying shape. Geophysics 2003, 68, 381–385. [Google Scholar] [CrossRef]
  9. McFadden, P.D.; Cook, J.G.; Forster, L.M. Decomposition of gear vibration signals by the generalized s transform. Mech. Syst. Signal Process. 1999, 13, 691–707. [Google Scholar] [CrossRef]
  10. Huang, N.E.; Shen, Z.; Long, S.R.; Wu, M.C.; Shih, H.H.; Zheng, Q.; Yen, N.C.; Tung, C.C.; Liu, H.H. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc. R. Soc. A Math. Phys. Eng.Sci. 1998, 454, 903–995. [Google Scholar] [CrossRef]
  11. Huang, N.E.; Wu, Z. A review on Hilbert-Huang transform: Method and its applications to geophysical studies. Rev. Geophys. 2008, 46, RG2006. [Google Scholar] [CrossRef] [Green Version]
  12. Daubechies, I.; Lu, J.; Wu, H.T. Synchrosqueezed wavelet transforms: An empirical mode decomposition-like tool. Appl. Comput. Harmon. Anal. 2011, 30, 243–261. [Google Scholar] [CrossRef] [Green Version]
  13. Dragomiretskiy, K.; Zosso, D. Variational mode decomposition. IEEE Trans. Signal Process. 2014, 62, 531–544. [Google Scholar] [CrossRef]
  14. Zhang, G.L. Time-phase amplitude spectra based on a modified short-time Fourier transform. Geophys. Prospect 2018, 66, 34–46. [Google Scholar] [CrossRef] [Green Version]
  15. Zhang, G.L.; Li, Y.; Wang, T.Y.; Du, H.; Luo, F.; Zhan, Y.Z. Extended Hilbert transform and application for seismic attributes. Earth Space Sci. 2019, 6, 873–886. [Google Scholar] [CrossRef]
  16. Li, Y.; Zhang, G.L.; Duan, J.; He, C.J.; Du, H.; Luo, F.; Zhan, Y.Z.; Wang, J. Extended stable factor method for the inverse Q-filter. Geophysics 2020, 85, T155–T163. [Google Scholar] [CrossRef]
  17. Auger, F.; Flandrin, P.; Lin, Y.T.; McLaughlin, S.; Meignen, S.; Oberlin, T.; Wu, H.T. Time-frequency reassignment and synchrosqueezing: An overview. IEEE Signal Process. Mag. 2013, 30, 32–41. [Google Scholar] [CrossRef] [Green Version]
  18. Herrera, R.H.; Han, J.J.; Baan, M.V.D. Applications of the synchrosqueezing transform in seismic time-frequency analysis. Geophysics 2014, 79, V55–V64. [Google Scholar] [CrossRef]
  19. Zhang, G.L.; Duan, J.; Li, Y.; He, C.J.; Zhan, Y.Z.; Wang, J. Adaptive time-resample high-resolution synchrosqueezing transform and its application in seismic data. IEEE Trans. Geosci. Remote Sens. 2020, 58, 6691–6698. [Google Scholar] [CrossRef]
  20. Li, Y.; Zhang, G.L. High-precision and high-resolution synchrosqueezing transform via time-frequency instantaneous phases. Appl. Sci. 2021, 11, 11760. [Google Scholar] [CrossRef]
  21. Marfurt, K.J. Robust estimates of 3D reflector dip and azimuth. Geophysics 2006, 71, 29–40. [Google Scholar] [CrossRef]
  22. Marfurt, K.J. Techniques and best practices in multiattribute display. Interpretation 2015, 3, B1–B23. [Google Scholar] [CrossRef]
  23. Robert, A. Curvature attributes and their application to 3D interpreted horizons. First Break 2001, 19, 85–100. [Google Scholar] [CrossRef]
  24. Nyquist, H. Certain topics in telegraph transmission theory. Trans. Am. Inst. Electr. Eng 1928, 47, 617–644. [Google Scholar] [CrossRef]
  25. Shannon, C.E. Communication in the presence of noise. Proc. IRE. 1949, 37, 10–21. [Google Scholar] [CrossRef]
  26. Bracewell, R.N. The Fourier Transform and Its Applications; McGraw-Hill Book Co., Inc.: New York, NY, USA, 1965; pp. 268–271. [Google Scholar]
  27. Li, B.; Xu, N.W.; Dai, F.; Zhang, G.L.; Xiao, P.W. Dynamic analysis of rock mass deformation in large underground caverns considering microseismic data. Int. J. Rock Mech. Min. Sci. 2019, 122, 104078. [Google Scholar] [CrossRef]
  28. Li, B.; Ding, Q.F.; Xu, N.W.; Dai, F.; Xu, Y.; Qu, H.L. Characteristics of microseismic b-value associated with rock mass large deformation in underground powerhouse caverns at different stress levels. J. Cent. South Univ. 2022, 29, 693–711. [Google Scholar]
Figure 1. (a) Synthetic 1D channel data. Instantaneous amplitudes obtained by: (b) HT; (c) GHT with p = 1.0; (d) GHT with p = 0.8 (ratio of an even number 4 and an odd number 5); (e) GHT with p = 0.6 (ratio of two odd numbers 3 and 5).
Figure 1. (a) Synthetic 1D channel data. Instantaneous amplitudes obtained by: (b) HT; (c) GHT with p = 1.0; (d) GHT with p = 0.8 (ratio of an even number 4 and an odd number 5); (e) GHT with p = 0.6 (ratio of two odd numbers 3 and 5).
Applsci 12 02836 g001
Figure 2. Instantaneous amplitudes of the synthetic 1D channel data in Figure 1a, which are obtained by LTFS with: (a) p =1.0, f1 = 0 Hz and f2 = fN = 500 Hz; (b) p = 0.8, f1 = 0 Hz and f2 = fN = 500 Hz; (c) p = 0.6, f1 = 0 Hz and f2 = fN = 500 Hz.
Figure 2. Instantaneous amplitudes of the synthetic 1D channel data in Figure 1a, which are obtained by LTFS with: (a) p =1.0, f1 = 0 Hz and f2 = fN = 500 Hz; (b) p = 0.8, f1 = 0 Hz and f2 = fN = 500 Hz; (c) p = 0.6, f1 = 0 Hz and f2 = fN = 500 Hz.
Applsci 12 02836 g002
Figure 3. (a) Synthetic 2D seismic data. Instantaneous amplitudes obtained by: (b) HT; (c) GHT with p = 0.4; (d) LTFS with p = 0.4, f1 = 0 Hz, and f2 = fN = 500 Hz; (e) GHT with p = 0.2 (or LTFS with p = 0.2, f1 = 0 Hz, and f2 = fN = 500 Hz); (f) LTFS with p = 0.2, f1 = 10 Hz, and f2 = 100 Hz.
Figure 3. (a) Synthetic 2D seismic data. Instantaneous amplitudes obtained by: (b) HT; (c) GHT with p = 0.4; (d) LTFS with p = 0.4, f1 = 0 Hz, and f2 = fN = 500 Hz; (e) GHT with p = 0.2 (or LTFS with p = 0.2, f1 = 0 Hz, and f2 = fN = 500 Hz); (f) LTFS with p = 0.2, f1 = 10 Hz, and f2 = 100 Hz.
Applsci 12 02836 g003aApplsci 12 02836 g003b
Figure 4. Zoomed-in views of yellow rectangles in Figure 3. (af) show the zoomed-in views of yellow rectangles in Figure 3a–f, respectively.
Figure 4. Zoomed-in views of yellow rectangles in Figure 3. (af) show the zoomed-in views of yellow rectangles in Figure 3a–f, respectively.
Applsci 12 02836 g004
Figure 5. Two profiles seismic data extracted from the input actual 3D seismic data. (a) Inline 160; (b) Crossline 96.
Figure 5. Two profiles seismic data extracted from the input actual 3D seismic data. (a) Inline 160; (b) Crossline 96.
Applsci 12 02836 g005
Figure 6. Time slices at 1.7s: (a) amplitudes; (b) instantaneous amplitudes; (c) instantaneous amplitudes obtained by GHT with p = 6; (d) instantaneous amplitudes obtained by LTFS with p = 6, f1 = 0 Hz, and f2 = 90 Hz; (e) most positive curvatures; (f) third-generation coherences (C3).
Figure 6. Time slices at 1.7s: (a) amplitudes; (b) instantaneous amplitudes; (c) instantaneous amplitudes obtained by GHT with p = 6; (d) instantaneous amplitudes obtained by LTFS with p = 6, f1 = 0 Hz, and f2 = 90 Hz; (e) most positive curvatures; (f) third-generation coherences (C3).
Applsci 12 02836 g006
Figure 7. Time slices at 1.7s with two interpreted channels: (a) amplitudes; (b) instantaneous amplitudes; (c) instantaneous amplitudes obtained by GHT with p = 6; (d) instantaneous amplitudes obtained by LTFS with p = 6, f1 = 0 Hz, and f2 = 90 Hz; (e) most positive curvatures; (f) third-generation coherences (C3).
Figure 7. Time slices at 1.7s with two interpreted channels: (a) amplitudes; (b) instantaneous amplitudes; (c) instantaneous amplitudes obtained by GHT with p = 6; (d) instantaneous amplitudes obtained by LTFS with p = 6, f1 = 0 Hz, and f2 = 90 Hz; (e) most positive curvatures; (f) third-generation coherences (C3).
Applsci 12 02836 g007
Table 1. Time costs of GHT and LTFS results in Figure 3e,f.
Table 1. Time costs of GHT and LTFS results in Figure 3e,f.
MethodsGHT with p = 0.2 (or LTFS with p = 0.2, f1 = 0 Hz, and f2 = 500 Hz)LTFS with p = 0.2,
f1 = 10 Hz, and f2 = 100 Hz
Time costs (s)T1 = 4.21T3 = 0.77
Ratios of time costsT1/T3 = 5.47
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, H.; Cao, J.; He, Z.; Li, Y.; Lin, K.; Zhang, G. Lp-Order Time–Frequency Spectra and Its Application in Edge Detection. Appl. Sci. 2022, 12, 2836. https://doi.org/10.3390/app12062836

AMA Style

Liu H, Cao J, He Z, Li Y, Lin K, Zhang G. Lp-Order Time–Frequency Spectra and Its Application in Edge Detection. Applied Sciences. 2022; 12(6):2836. https://doi.org/10.3390/app12062836

Chicago/Turabian Style

Liu, Houjun, Junxing Cao, Zhenhua He, Yong Li, Kai Lin, and Gulan Zhang. 2022. "Lp-Order Time–Frequency Spectra and Its Application in Edge Detection" Applied Sciences 12, no. 6: 2836. https://doi.org/10.3390/app12062836

APA Style

Liu, H., Cao, J., He, Z., Li, Y., Lin, K., & Zhang, G. (2022). Lp-Order Time–Frequency Spectra and Its Application in Edge Detection. Applied Sciences, 12(6), 2836. https://doi.org/10.3390/app12062836

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop