This is the quick test of F-J method.
https://doi.org/10.1029/2018JB016595
F-J method is a kind of seismic method for extracting surface wave dispersion.
- freq_bessel_trans.m is the matlab program for F-J transform.
- PhaseShiftOfSW.m is the matlab program for phase shift method to extract surface wave dispersion.
- seis-1 / seis-2 or other files like that are the synthetic data generated by QSEIS.
- Uf.dat / Ur.dat or other files like that are the synthetic data generated by H.M Zhang's program.
There are some paper about F-J method.
- 基于矢量波数变换法_VWTM_的多道Rayleigh波分析方法_杨振涛
- Frequency-Bessel Transform Method for Effective Imaging of Higher-Mode Rayleigh Dispersion Curves From Ambient Seismic Noise Data_Jiannan Wang
We choose four methods to generate the surface wave data. In order to replicate the results of paper above, we choose same model which is presented below.
- QSEIS
- H.M. Zhang's Code
- SOFI2D
- SPECFEM2D
- complete all the description of program and data file
- check the validity of F-J method
- find out a reliable surface wave data, whatever synthtic or real, and treat it as benchmark.
- the speed of calculation is too slow. Nanqiao provide a python version of F-J method program. With the help of scipy, it runs much faster than before.
- I still do NOT figure out how to conduct the F-J transform. There are much different version of formula, and I am sure that there must be some typo in the paper especially in chinese version.