数字处理第一章.ppt

在图1-24上理想低通滤波器I 的时间特性为 理想低通滤波器I I的时间特性为 理想带通滤波器的时间特性为：h(t)=h1(t)-h2(t) 所以 最后得到 理想带通滤波器的时间特性 见图1-26 5）褶积滤波的具体计算 褶积滤波的具体计算步骤如下： （1）对地震记录进行频谱分析，确定通频带中心频率f0和带宽2 △f。 （2）确定滤波因子长度N。理论上，滤波因子是无限长的。实际上，要在计算机上计算，不可能去无限长，而是取某一长度N。滤波因子长度N的选取很重要，即使滤波效果较好，又要节省计算工作量。一般根据实验工作来确定，现在多用N=101. (3)求滤波因子，如果是带通滤波器，则将 右式写成离散形式： （4）将 写成离散形式 用左式编制程序，对地震记录道可进行单道滤波。在实际运算时，应根据滤波因子的长度对滤波钱的地震信号的前后补适当个数的零值。 Interactive Horizon RMO Analysis View CRP Gathers MVA Xsection Window MVA Xsection Window displays stacked data, velocity fields, and allows you to pick horizons. View CRP Gathers allows you to look at migrated CRP gathers by clicking on the desired CRP location in the MVA Xsection Window. IHVA allows you to pick RMO values that can be used to update a velocity field. View Ray Density Difference Files from MVA Tomography Inversion Small inversion aggressiveness factor Large inversion aggressiveness factor VSP Wave Field Separation o Flattening the downgoing using the first arrivals o Wavefield Separation using a Median Filter o Wavefield Separation using an F-K filter o Wavefield Separation using an EigenVector Filter. o Wavefield Separation Test Comparisons o Saving the Upgoing to Disk o Saving the Downgoing to Disk VSP Wave Field Separation Hodogram of Two Horizontal Traces Hodogram of Oriented Horizontal and Original Vertical 频率域滤波原理 设有一滤波器如图1-8所示，其滤波器时间函数或滤波因子h(t)的频谱H(w)称为滤波器的频率特性，输入信号x(t)的频谱X(w)，输出信号x*(t)的频谱X*(w)，根据滤波理论，在线性滤波条件下，滤波器输出信号的频谱X*(w)为输入信号的频谱X(w)与滤波器频率特性H(w)的乘积，即 X*(w)= X(w) H(w) 频率域滤波就是利用这个方程进行运算，已知两个量就可求另一个量。一般是输入记录已知，滤波器是根据需要设计的，而要求的则是经过滤波后的输出。为此将输入信号的频谱X(w)和滤波器的频率特性H(w)相乘即可求出。 1）对滤波器振幅频率特性和相位频率特性的要求 设输入信号的频谱为X(w)=|X(w)|exp[iφx (w)] 设输出信号的频谱为X*(w)=|X*(w)|exp[iφ*x (w)] 滤波器的频谱为H(w)=|H(w)|exp[iφH (w)] 由线性滤波方程有：|X*(w)|=|X(w)||H(w)| φ*x (w)= φx (w)+ φH (w)] 上式表明，输出信号的振幅谱等于输入信号的振幅谱与滤波器的振幅频率特性的乘积，输出信号的相位谱等于输入信号的相位谱与滤波器相位特性之和。 （2）对滤波器H(w)的要求 如果有效波的频谱是在w1—w2的频率范围，则滤波器的振幅频率特性在w1—w2范围内应为常数（一般取为1），即 |H(W)|