地震资料处理中的反褶积处理.ppt

地震资料处理中的反褶积处理 地震资料处理中的反褶积处理 反褶积概述（2001年已讲） 预备知识 预测反褶积的基本原理和计算方法 Omega系统的主要反褶积处理模块 反褶积概述 褶 积 地震记录的褶积模型 地震记录的分辨率 反褶积的定义 反褶积的类型 预备知识 信号 信号的频谱 信号的离散化 Z变换 信号的相位特征 相关分析 两种特殊信号 反信号 信号 什么叫信号？ 随时间变化的物理量称为信号。即，信号是以时间为自变量的函数。 物理可实现信号： 信号x(t),当t0时，x(t)=0 地震记录是物理可实现信号 信号的频谱 在某个时间区间内的地震记录通常是如下图所是的一个比较复杂的信号。 为了说明频谱的意义，先从最简单的信号说起。在物理学中，最简单的波是描述简谐振动的简谐波，它可用正弦函数表示： Asin(2πft+φ) 。 A：振幅，f: 频率， φ：初始相位（简称相位） 在某个区间内最简单的简谐波是Asin（2πf0（t-t0)+ φ）。其中， f0 =1/T 该简谐波的周期正好是时间区间长度T。 稍复杂一点的是频率为fn=nf0的正弦波： Asin(2πnf0 （t-t0) +φ) ，这时正弦波的周期是T/n f0 ：基频， Asin(2πnf0 （t-t0) +φ) ：n次谐波。一次谐波称为基波。 频谱的计算 用付氏变换公式： 正变换 反变换 X(f)是复数,称为x(t)的复频谱,它可表示为： X(f)=u(f)+iv(f) 由此可得： 振幅谱: 相位谱: 信号的离散化 实际地震记录是连续信号，数字仪记录时，要间隔一定的时间间隔Δ记录一个值，由此将地震记录x(t)变成时间序列 x(nΔ) (n=1,2,…N) Δ 称为采样间隔。 这就是信号的离散化。 对于离散化有以下采样定理： 若连续信号x(t)有截止频率fc，则当 时，离散x(nΔ) 可完全确定X(t): 奈魁斯特频率 如果x(t)不存在截止频率fc，或 时， x(nΔ) 不能完全恢复x(t)，但x(t)的频谱X(f)与x(nΔ) 的频谱xΔ(f)之间有以下关系： 该式表明： 1、 xΔ(f) 是一个周期函数（周期为1/Δ） 2、它在一个周期的值等于将X(f)以为 基础分为若干 小段，每段长1/ Δ，然后将各段的X(f)值相加。 由此可见，当采样率为Δ 时，离散序列的最大频率为1/2Δ，这就是奈魁斯特频率，也称折叠频率。 Z变换 序列(a0,a1,a2,…an)的Z变换定义为 A(z)= a0+a1z,+a2z2+…anzn (z是复数) 使A(z)= 0的z值称为Z变换的根，该序列的Z变换有n个根。 信号的相位特征 设一两项信号 a=(a1,a2),则 1、若a1a2,称a是最小相位延迟信号 2、若a1a2,称a是最大相位延迟信号 3、若a1=a2,称a是等延迟信号 任一n+1项信号 b=(b0,b1,…,bn)可分解为n个两项信号 的褶积。 如果 1、所有两项信号 都是最小相位延迟信号，则b是最小相位 2、所有两项信号 都是最大相位延迟信号，则b是最大相位 3、既有最大相位延迟也有最小相位延迟，则b是混合相位 信号的相位特征也可用其z变换来定义： 1、 z 变换的根都在单位圆外，信号是最小相位 2、 z 变换的根都在单位圆内，信号是最大相位 3、单位圆内外都有根，信号是混合相位 相关分析 相关函数的定义 相关与褶积的关系 相关函数的频谱 相关函数的定义 1、互相关函数 2、自相关函数 （Xn、yn为离散信号） 相关与褶积的关系 信号xn与gn的褶积为：xn*yn= 信号xn与yn的相关函数： 两个信号的互相关函数等于将后一个信号的翻转信号与前一信号的褶积： 相关函数的频谱 由： 有： 或 特别地，对于自相关函数有：Rxx(f)=|X(f)|2 以上公式说明： 1、自相关函数的频谱是实数； 2、由信号的振幅谱可确定其自相关函数的频谱进而确定自相关 函数。反过来，由自相关函数也可求振幅谱。 两种特殊信号 1、单