< p >我有一些简单的地震跟踪如下所示。< a href = " https://i.stack.imgur.com/xDEFP.png " rel = " nofollow noreferrer " > < img src = " https://i.stack.imgur.com/xDEFP.png " alt = "地震采集" / > < / > < / p > < p >当试图使用这个描述的技术执行色散分析< a href = " https://expservices.ku.edu/sites/expservices.drupal.ku.edu/files/docs/1998/Park%20et%20al%20 (1998) _Imaging % 20色散曲线% 20 % 20 % 20的表面波% 20 % 20 % 20多通道% 20 record.pdf”rel = " nofollow noreferrer " > < / >,得到如下图。< / p > < p > < a href = " https://i.stack.imgur.com/L37Pz.png " rel = " nofollow noreferrer " > < img src = " https://i.stack.imgur.com/L37Pz.png " alt = "分散分析" / > < / > < / p > < p >不幸的是,我不能对我的生活算出我做错了什么。所以我觉得我一定遗漏了一些关键点概念上。下面是我的python脚本处理离散图像。< / p > < pre > <代码> # #。。执行色散分析# #。定义参数ng = 20 #铁板一步数量= 200 #的副总裁测试点nfreq = 200 #的频率测试点ndata = 600 #每个跟踪的数据点数量dt = 1/16000 # Time_Delta # #。。定义gridspace频率= np。linspace (0 40 nfreq) #频谱(Hz) Vp = np。一步法linspace(500) #相速度谱(m / s) X, Y = np。副总裁meshgrid(频率)# #。。设置gridspace倾斜叠加函数S = np。zeros((ng, nfreq, nVp, ndata), dtype=complex) ## . . Calculate dispersion with varying phase velocity and frequency for idx, trace in enumerate(seismic_gather[0].trace): ## Get the fourier transformation and normalization of the waveform U = np.fft.fft(trace.data) # Fourier transform N = U/np.abs(U) # Normalize waveform ## Apply dynamic linear moveout P = 2*np.pi*X*trace.offset/Y # Dispersion property S[idx] = np.exp(1j*P)[:,:,np.newaxis]*N ## . . Stack traces S = np.sum(S, 0) ## . . Extract amplitudes A = np.abs(S)/ng To describe the code a little, I am defining a meshgrid space of 200 testing points for the Phase Velocity and Frequency. I then take my seismic data and transform it into the frequency domain and normalize it. Once normalized, I apply the dynamic linear moveout to each trace, and then stack the the traces for each pair of phase velocity and frequency points along the meshgrid.
Unfortunately, the resulting image doesn't appear to be what I expect, which I expect to be similar to the paper linked. Am I missing a step here?