油气勘探方法程序设计课程设计(论文)-最小平方反滤波.docxVIP

油气勘探方法程序设计课程设计(论文) 设计（论文）题目 最小平方反滤波 学院名称 地球物理学院 专业名称 勘查技术与工程（石油物探） 学生姓名 学生学号 任课教师 田仁飞 设计（论文）成绩 教务处 制 2015年7月1日  PAGE \* MERGEFORMAT 22 最小平方反滤波 一、方法原理 由地震波的传播理论可知，在粘弹性介质中地震波是以地震子波的形式在地下传播，地面接收到的反射波地震记录是地层反射系数与地震子波的褶积，因此，地层相当于一个滤波器，是反射系数序列变成了由子波组成的地震记录，降低了地震勘探的纵向分辨率。反滤波的目的就是要设计一个反滤波器，来对地震记录滤波，消除地层滤波的作用，提高地震记录的纵向分辨率。 最小平方反滤波的基本思想在于设计一个滤波算子，用它把已知的输入信号转换为与给定的期望输出信号在最小平方误差的意义下是最佳输出。 1）用最小平方法求反滤波因子 对输入子波b(t)反滤波后的期望输出为d(t),实际输出为y(t)，按最小平方原理，使二者的误差平方和Q为最小时求得的反滤波因子称为最小平方反滤波因子，用它对地震记录x(t)进行的反滤波为最小平方反滤波。 设输入离散信号为地震子波b(n)={b(1),b(2), …,b(m)},待求的反滤波因子a(n)={a(0),a(1),a(2), …,a(m)}, a(t)的起始时间为0，（m+1）为a(t)的延续长度，b(n)与a(n)的褶积为实际输出y(n),即 yn=an*bn=τ=0maτbn-τ 1 实际输出与期望输出的误差平方和为 Q=n=0m[yn-d(n)]2 =n=0m[τ=0maτbn-τ-d(n)]2 （2） 要使Q为最小，数学上就是求Q的极值问题，即求满足 ?Q?a(l)=0 l=0,1,…,m （3） 的滤波因子a(t)。 即满足方程 τ=0maτrbbl-τ=rbdl l=0,1,…,m （4） 其中rbbl-τ=n=0mbn-τb(n-l)为地震子波的自相关函数，rbdl=n=0md(n)b(n-l),写成矩阵形式 rbb(0)rbb(1)?rbb(m)rbb(1)…rbb(m)rbb(0)…rbb(m-1)?…?rbb(m-1)…rbb(0)a0a1?am=adb0adb1?adbm （5） 上式系数矩阵称为托布利兹方程。 当期望输出是单位脉冲δt时，即 dt=δt=1,t=00,t≠0 （6） 则 rbdl=b-l （7） 由于b(t)是地震子波，b(t)=0,当t0时，上述矩阵可写成： rbb(0)rbb(1)?rbb(m)rbb(1)…rbb(m)rbb(0)…rbb(m-1)?…?rbb(m-1)…rbb(0)a0a1?am=b00?0=10?0 （8） 用地震记录自相关系数rxx(τ)代替ybb(τ) 2）已知地震记录x(t) x(t)=b(t)* ζ(t) （9） 其中x(t)为地震记录，ζ(t)是反射系数。根据自相关函数的性质，可证明：地震记录自相关rxx(n)是地震子波自相关rbb(n)和反射系数自相关rζζ(n)的褶积,