资料处理中的反褶积

《资料处理中的反褶积》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、地震资料处理中的反褶积处理黄大云2002年5月地震资料处理中的反褶积处理反褶积概述（2001年已讲）预备知识预测反褶积的基本原理和计算方法Omega系统的主要反褶积处理模块反褶积概述褶积地震记录的褶积模型地震记录的分辨率反褶积的定义反褶积的类型预备知识信号信号的频谱信号的离散化Z变换信号的相位特征相关分析两种特殊信号反信号信号什么叫信号？随时间变化的物理量称为信号。即，信号是以时间为自变量的函数。物理可实现信号：信号x(t),当t<0时，x(t)=0地震记录是物理可实现信号00tt信号的频谱在某个时间区间内的地震记录通常是如下图所是的一个比

2、较复杂的信号。为了说明频谱的意义，先从最简单的信号说起。在物理学中，最简单的波是描述简谐振动的简谐波，它可用正弦函数表示：Asin(2πft+φ)。A：振幅，f:频率，φ：初始相位（简称相位）在某个区间内最简单的简谐波是Asin（2πf0（t-t0)+φ）。其中，f0=1/T该简谐波的周期正好是时间区间长度T。稍复杂一点的是频率为fn=nf0的正弦波：Asin(2πnf0（t-t0)+φ)，这时正弦波的周期是T/nf0：基频，Asin(2πnf0（t-t0)+φ)：n次谐波。一次谐波称为基波。t0t0+Tt0t0+TΦ=0t0t0+TΦ≠0多

3、个简谐信号相加可得到一个复杂信号：一个复杂信号可分解为多个简谐信号的叠加：各简谐信号的振幅构成信号的振幅谱A(f)，各简谐信号的初始相位构成信号的相位谱Φ(f)，振幅谱和相位谱合称频谱。频谱可用复数表示：频谱表示成复数后可做数学运算X(t)=x1(t)+x2(t)频谱的计算用付氏变换公式：正变换反变换X(f)是复数,称为x(t)的复频谱,它可表示为：X(f)=u(f)+iv(f)由此可得：振幅谱:相位谱:信号的离散化实际地震记录是连续信号，数字仪记录时，要间隔一定的时间间隔Δ记录一个值，由此将地震记录x(t)变成时间序列x(nΔ)(n=1,2

4、,…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Δ，这就是奈魁斯特频率，也称折叠频率。频率折叠示

5、意图Z变换序列(a0,a1,a2,…an)的Z变换定义为A(z)=a0+a1z,+a2z2+…anzn(z是复数)使A(z)=0的z值称为Z变换的根，该序列的Z变换有n个根。信号的相位特征设一两项信号a=(a1,a2),则1、若a1>a2,称a是最小相位延迟信号2、若a1

6、延迟，则b是混合相位信号的相位特征也可用其z变换来定义：1、z变换的根都在单位圆外，信号是最小相位2、z变换的根都在单位圆内，信号是最大相位3、单位圆内外都有根，信号是混合相位相关分析相关函数的定义相关与褶积的关系相关函数的频谱相关函数的定义1、互相关函数2、自相关函数（Xn、yn为离散信号）相关与褶积的关系信号xn与gn的褶积为：xn*yn=信号xn与yn的相关函数：两个信号的互相关函数等于将后一个信号的翻转信号与前一信号的褶积：相关函数的频谱由：有：或特别地，对于自相关函数有：Rxx(f)=

7、X(f)

8、2以上公式说明：1、自相关函数的频谱

9、是实数；2、由信号的振幅谱可确定其自相关函数的频谱进而确定自相关函数。反过来，由自相关函数也可求振幅谱。两种特殊信号1、单位脉冲δ(t)(狄拉克函数）（当t=0时）（当t≠0时）δ(t)频谱是Δ(f)=12、白噪声b(t)∑b(t)=0Rbb(t)=δ(t)反信号对信号x(t)，如果有信号a(t)，使x(t)*a(t)=δ(t)，则称a(t)是x(t)的反信号。由于写成指数形式：所以，反信号的频谱与原信号的频谱有以下关系：1、2、φx=-φa并非任何信号都有反信号。如在某些频率点f，使，则反信号不存在。预测反褶积的基本原理和计算方法脉冲反褶

10、积预测反褶积的基本原理和计算方法脉冲反褶积1、脉冲反褶积的假设条件2、脉冲反褶积的基本原理3、脉冲反褶积的计算1、脉冲反褶积的假设条件两个假设条件（1）反射系数函数