我有一些简单的地震跟踪如下所示。
当试图执行一个分散分析使用描述的技术纸如下图,我得到。
不幸的是,我的生活不能找出它是什么我做错了。所以我觉得我一定遗漏了一些关键点概念上。下面是我的python脚本处理离散图像。
# #。执行色散分析# #。定义参数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。0 (ng nfreq,一步法,ndata), dtype =复杂)# #。计算色散不同idx的相速度和频率,跟踪在列举(seismic_gather [0] .trace): # #的傅里叶变换和归一化波形U = np.fft.fft (trace.data) #傅里叶变换N = U / np.abs (U) # # #正常波形应用动态线性时差P = 2 * np.pi * X *痕迹。抵消/ Y #分散产权S [idx] = np.exp (1 j * P) [::, np。newaxis] * N # #。。堆栈跟踪S = np。(年代,0)# #。。 Extract amplitudes A = np.abs(S)/ng描述代码,我定义一个meshgrid空间200个测试点的相速度和频率。然后我把我的地震数据并将其转换成频域和规范化。一旦规范化,我对每个跟踪应用动态线性时差,然后堆栈的痕迹每一对沿着meshgrid相速度和频率点。
不幸的是,由此产生的图像不似乎是我所期望的,我希望是类似于纸联系在一起。我缺少一个步骤吗?
