地震滤波.pptVIP

一、一维频率域滤波 2、频率域滤波的实现 1）实现步骤 首先要对地震记录x(t)作傅氏变换，得到其频谱X(ω)，进行频谱分析。根据有效波的频带范围，设计合适的滤波器H(ω) ，在频率域进行滤波。然后对其输出 做傅氏反变换，得到滤波后的 。整个过程可归结为下列运算过程(见图2-15): §2.2 一维滤波 所用的傅氏变换有离散傅氏变换DFT和快速傅氏变换，由于DF}I速度慢，实际应用中采用FFT。 一、一维频率域滤波 2、频率域滤波的实现 3）使用FFT应注意的问题 (1)输入数据：输入数据点数NFFT应是2k个点，当实际点数不够时，应补零。当输入数据是实数时、直接将子程序中复型数组CX的实部赋值，虚部赋零。 2)输出数据：计算出的频谱共NFFT个点(第一个正周期内的值)。从第1点开始，以NFFT/2+1处为对称中心，与后面的点有共轭关系：实 部相同，虚部差一符号。图2-l6示意地画出了振幅谱，其对称点在 NFFT/2+1处，它实际上是下一个周期的第一个值，即相当于图中第一个点的值。 §2.2 一维滤波 一、一维频率域滤波 2、频率域滤波的实现 3）使用FFT应注意的问题 (3)输入与输出数据采样间隔间的关系： 从（2-30）式中也可以看出振幅谱的对称点正好在折叠频率处，而第NFFT+1个点，则重复第1个点的值，对应于采样频率处。 §2.2 一维滤波 一、一维频率域滤波 2、一个实际应用例子 设原信号采样问隔为4us.样点数为300个，如图2-17所示。图2-18是它的频谱，其中分布有大量的高频噪声。现在设计一低通滤波器，高截频为30Hz，如图2-19，注意它是对称的。图2-20是它的时域响应，他是一抽样函数。对上述信号进行频率域滤波，时域响应如图2-21所示。 §2.2 一维滤波 二、一维时间域滤波 在地震勘探中，地震仪记录下来的是地震波随时间变化的信号。在前一节里介绍的频率滤波是先将信号通过FFT转换到频率域。频率域滤波后，再进行IFFT变为时间序列。那么，能否不进行变换直接在时间域进行处理呢？这就是时间域滤波问题。 时间域滤波有两种常用的方法，即褶积滤波和递归滤波。由于前面已经介绍了频率滤波，由傅氏变换与时域褶积的关系，很容易过渡到褶积滤波，下面介绍时域褶积滤波： §2.2 一维滤波 二、一维时间域滤波 1、褶积滤波 1）地震记录的褶积模型 设有一地震子波b(t),经过地下几个反射界面(反射系数函数考ξ(t))中的时间已用双层旅行时表示)，具体的褶积过程可以从图2-22看到，该子波每遇到一个界面，就产生一个反射波，其幅度、极性由反射系数决定，各输出的响应叠加起来，就得到了褶积结果(地震记录x(t) )，它实际上就是无噪声时的合成地震记录。需要指出的是，时域褶积除滤彼外，另一个重要应用就是制作合成地震记录。 以上过程可用下式表示 §2.2 一维滤波 如果考虑噪声存在，则褶积模型变为 二、一维时间域滤波 1、褶积滤波 1）地震记录的褶积模型 §2.2 一维滤波 二、一维时间域滤波 1、褶积滤波 1）地震记录的褶积模型 §2.2 一维滤波 二、一维时间域滤波 1、褶积滤波 3）褶积滤波应注意的问题 (1)滤波因子的设计要求 考虑到滤波器频响H(ω)是一非负的实偶函数，滤波因子h(t)应当取一偶函数，即h(t)关于t=0对称。 (2)褶积结果的长度 设输入x(t)是M个点，滤波因子h(t)为N个点，则滤波结果为L=M+N-1个点。一般的地震信号褶积滤波后，结果只需要原来长度就行了，这时可各舍去结果两端(N-1)/2个点(因为h(t)为偶函数，N取奇数)。 §2.2 一维滤波 二、一维时间域滤波 1、褶积滤波 3）褶积滤波应注意的问题 (3) 褶积运算的速度问题 褶积滤波和频率滤波都能实现一维滤波问题。究竟选用哪种方法，主要应考虑运算速度问题。对于较长地震记录，最好用频率域滤波，而且在多道处理时运算速度优势将更加明显。对于短记录或少量的数据，时域的褶积滤波是可取的。 §2.2 一维滤波 二、一维时间域滤波 2、递归滤波 除褶积滤波外，有一种时间域滤波方法——递归滤波，它比褶积滤波和频率滤波有更高的运算速度，这对人批量的地震数据处理工作，优势是不言而喻的。另外，从原理上来说，递归方法是一种常用的数学方法，不光是在地震数据处理中，在其它很多领域中都要用到它。因此，学习其基本原理和实现方法，是很有必要的。 1）基本原理 §2.2 一维滤波 对于褶积运 二、一维时间域滤波 2、递归滤波 1）基本原理 如果滤波因子h(n)的长度为N个点，那么每计算