SeismicLab地震数据处理包帮助文档归纳翻译.docVIP

SeismicLab地震数据处理包 帮助文档归纳翻译 一．源代码Codes （一）BP-Filter——Apply a band-pass filter to a group of traces. freq. domain filtering （带通-频率域滤波） function [o] = bp_filter(d,dt,f1,f2,f3,f4); 函数注释： IN d: data (columns are traces) %数据 dt: sampling interval in sec %采样间隔（秒） f1: freq. in Hz %最小频率 f2: freq. in Hz %f2-f3表示频带宽度 f3: freq. in Hz %f2-f3表示频带宽度 f4: freq. in Hz %最大频率 | ___________ | / \ Amplitude spectrum 振幅谱 | / \ %镶边法用于避免吉布斯现象造成的影响 | / \ |------------------------------- f1 f2 f3 f4 OUT o: output (columns are traces) %输出数据（地震道） Example: d=linear_events; dout = bp_filter(d,0.004,1,3,30,40); wigb([d,dout]); %wigb画图的函数，[d,dout]表示横向画图在一张figure上。 （二）decon——提供的反褶积方法有8种 time domain deconvolution and wavelet estimation（时域反卷积与子波评价） 【一】kolmog.m KOLMOG: Kolmogoroff spectral factorization。 % Given a wavelet, this function retrieves the minimum phase wavelet using Kolmogoroff factorization.If the input is a trace, the spectral factorization is applied to the autocorrelation after smoothing. %给定一个小波，该函数检索采用Kolmogoroff分解的最小相位子波。如果输入是一个地震道，在平滑后谱分解再应用于自相关。 Function [w_min] = kolmog(w,type_of_input,L); 函数注释： IN w: a wavelet of arbitrary phase if type_of_input = w %任意相位的子波 or a seismic trace if type_of_input = t %一个地震道 L: lenght of wavelet if type_of_input=t %子波长度 OUT w_min: a min phase wavelet %输出一个最小相位子波 Example: w = [1;2;-1;0;0;0;0] wmin = kolmog(w,w); subplot(221); stem(w); subplot(222); stem(wmin); 【二】ls_inv_filter.m LS_INV_FILTER: Least-squares inverse filter of a wavelet. %一种子波的最小二乘法反滤波。 function [f,o] = ls_inv_filter(w,NF,Lag,mu); 函数注释： IN w: the wavelet NF: lenght of the inverse filter %反滤波的长度 Lag: the position of the spike in the desired output %在期望输出中尖峰的位置 Lag=1 for minimum phase w