地震资料处理复习总结文档格式.docx
《地震资料处理复习总结文档格式.docx》由会员分享,可在线阅读,更多相关《地震资料处理复习总结文档格式.docx(11页珍藏版)》请在冰豆网上搜索。
其模量为函数X(x,t)的振幅谱。
如果有效波和干扰波得平面简谐波成分有差异,有效波的平面简谐波成分与干扰波的平面简谐波成分不同的视速度传播,则可以用二维视速度滤波将他们分开,达到压制干扰,提高性噪比的目的。
二维频谱的特征:
空间假频
数字滤波
在地震勘探中,用数字仪器记录地震波时,为了保持更多的波得特征,,通常利用宽频带进行记录,因此在宽频带范围内记录了各种反射波的同时,也记录了各种干扰波。
有效波和干扰波得差异表现在多个方面(频谱、传播方向、能量……)。
利用频谱特征的不同来压制干扰波,以突出有效波的方法就是数字滤波。
滤波器的响应特性:
对滤波器能力的最普遍度量是其响应特性
滤波器的频率特性:
其滤波器时间函数或滤波因子的频谱称为滤波器的频率特性,
滤波器的时间特性(单位脉冲响应):
在时间域的表示方法中,令一个单位脉冲通过一个滤波器,然后观侧滤波器的输出,这个滤波器输出的自然过程曲线称为滤波器的时间特性。
也称滤波器的“脉冲响应”
频率响应函数应该就是
时间和频率响应函数合起来应该就是就是响应特征
滤波机理:
输出信号的振幅谱等于输入信号的振幅谱与滤波器的振幅频率特性的乘积,
输出信号的相位谱等于输入信号的相位潜与滤波器相位特性之和。
(频率)
时间域上就是褶积
褶积滤波的物理意义:
它相当于把地震信息分解为起始时间、极性、幅度各不相同的脉冲序列,令这些脉冲按时间顺序依次通过滤波器,这样在滤波器的输出端就得到对输入脉冲序列的脉冲响应,这些脉冲响应有不同的起始时间,不同的极性和不同的幅度(这个幅度是与引起它的输入脉冲幅度成正比的),将它们叠加起来就得到滤波后的输出分。
频率域滤波的步骤
(1)对已知地震记录道进行频谱分析。
(2)设计合适的滤波器
(3)进行滤波运算
(4)对输出信号谱进行傅里叶反变换
褶积滤波的具体计算
褶积滤波的具体计算步骤如下:
(1)对地震记录进行频谱分析,确定通频带中心频率和带宽。
(2)确定滤波因子长度N。
(3)求滤波因子。
(4)将式
写成离散形式
数字滤波的离散型造成伪门现象,
数字滤波的有限性造成吉布斯现象。
伪门:
对连续的滤波因子用时间采样间隔离散采样后,得到。
如果,再按计算出与它相应的滤波器的频率特性,这时在频率特性的图形上,除了有同原来的对应的“门”外,还会周期性地重复出现很多“门”,这些门称为“伪门”。
为了避免“伪门”造成的影响,可以适当地采样间隔使第一个“伪门”出现在干扰波的频谱范围之外。
吉卜斯现象
有限长度的对应的不再是一个门式滤波,而是有波动的曲线(图1-31),曲线由间断点向远处波动衰减,在间断点波动最大,这种现象叫作非连续函数频率响应的吉卜斯现象。
实际计算中脉冲响应函数只能取有限长,截断后的脉冲响应所对应的频率响应函数不再是一个理想的“门”,而是接近于这个门得一条幅值有波动的光滑曲线。
镶边法解决。
类型:
低、带、高通
低通:
高通
带通:
低通和高通中间。
滤波器的分类:
按滤波器的性质(即响应函数)划分
1)无畸变滤波器。
振幅特性为常数,相位特性是线性的滤波器。
这种滤波器不改变输入信号的波形,
频率响应为、均为常数
2)相位畸变滤波器(纯相位滤波器,全啊通滤波器)
只改变输入信号的相位谱,振幅形状不改变。
振幅特性常数,但相位特性不是线性。
3)振幅畸变滤波器
振幅特性不是常数,实际工作中希望滤波时不是信号产生相位畸变或相位移,为零相位滤波器。
滤波器的稳定性:
当输入信号为有限,其输出信号也为有限,这种滤波器就是稳定的。
充要条件:
物理可实现性:
满足因果律(输入之前不产生输出)的滤波器称为物理可实现的。
当t<
0时
对于z变换为多项式的滤波器来说,
对于z变换为有理分式的滤波器来说,
当所有的根均不在单位圆(|z|=1)上时。
这个滤波器是稳定的,
当所有根都在圆外时,这个滤波器是物理可实现的。
反滤波作用
反褶积的基本作用是压缩地震记录中的地震子波,同时,可以压制鸣震和多次波,因而反褶积可以明显提高地震的垂直分辨率。
反褶积处理是褶积处理的反过程,因而称为反褶积。
如果我们设计一个滤波器,其滤波因子具有与滤波器有恰好相反的性质,即当输入为地震道记录约时,其输出为地层脉冲响应。
称这个反过程为反滤波或反褶积。
由震源发出的尖脉冲经过大地滤波器的滤波作用后,变成一个具有一定时间延续的波形,通常叫作地震子波
实际地震记录:
除了所表示的一系列反射波而外,还存在着干扰波,因此,地震记录的一般模型可以写为
最小平方反滤波两个假设条件:
1.反射系数为白噪声
2.假设随机噪声干扰为白噪声
子波未知情况下的假设条件:
1.反射系数为随机的白噪序列
2.假设地震子波是最小相位
(预测误差滤波也是这两个假设条件)
预测反滤波
预测滤波就是要设计一个预测因子,对输入已知的过去值和现在值进行滤波处理,获得未来某个时刻时的预测值
用预测的方法,根据地震记录一次反射和干扰的信息预测出纯干扰部分,再由包括一次波和干扰的地震记录中减去纯干扰部分,得到消除干扰后的一次反射信号,以消除一次反射后面的海上鸣震等多次波干扰。
叠加的目的是压制干扰,提高地震数据的信噪比。
动校正的目的是消除炮检距对反射波旅行时的影响
动校正:
将不同炮检距的反射时间校正到零炮检距反射时间的过程称为动校正
动校正量随炮检距递增,随反射深度和速度递减
动校正拉伸:
数字动校正造成的波形拉伸称为动校正拉伸
由于浅层的动校正时差大于深层的动校正时差,所以
浅层大炮检距的拉伸最为严重。
校正后反射波的形状不同,发生畸变,畸变的一般规律是反射波得波形被拉长,周期加大,频谱向低频方向移动。
拉伸系数:
反射深度越浅,炮检距越大(越大),动校正拉伸越严重,子波的主频向低频转移也随之严重。
静校正:
静校正是校正以及消除由于地表高程和地下低、降速带变化对反射波旅行时的影响。
野外一次静校正又称为基准面静校正,包括井深校正,地形校正,低速带校正。
人为选定一个海拔高程作为基准面,利用野外实测得到的各点高程、低速带厚度、速度、或井口时间等资料,将所有的炮点和检波点都校正到此面上,用基岩速度代替低速带速度,从而去掉表层因素的影响。
定义为:
用于补偿由于地表高程变化、风化层的厚度和速度变化对地震资料的影响。
其目的是获得在一个平面上进行采集,且没有风化层或低速介质存在时的反射波到达时间。
静校正的目的是获得在一个平面上进行采集,且没有风化层或低速介质存在时的反射波到达时间.
静校正信息来自于两个方面:
一是野外测量和观测的数据,包括地面高程数据、井口检波器记录时间、微测井和小折射数据等;
称为基准面校止或野外静校正.
二是根据初至波时间和地下反射信息求取静校正量,称为初至折射静校正和反射波地表一致性剩余静校正。
地震数据被校止到参考基准面上,消除了地表起伏和风化层横向变化的影响,后续地震处理工作就好像地震数据是在基准而上采集的。
有时我们不是将地震数据一次校正到参考基准面或最终基准面上,而是首先将地震数据校正到一个中间基准面上,这个基准面有时也称为浮动基准面或CMP叠加基准面。
速度分析、剩余静校正、动校正、叠加都在这个基准面上进行。
叠加之后,再将地震数据由浮动基准面校正到参考基准面或最终基准面上,在最终基准面上完成偏移处理。
静校正概念中“静”的含义
是相对动校正中“动”的含义而言的,我们知道,地震道的动校正时差是反射时间的函数,而地震道的静校正时差与地震道的时间无关,无论是浅层反射,还是深层反射,整个地震道只有一个静校正量。
基准面静校正
也称为野外静校正,顾名思义,就是将在地表采集的地震记录校正到基准面上,消除地表高程和风化层对地震记录旅行时的影响。
第一步是剥去风化层的影响
第二步,再将地震记录时间由风化层的底界校正到参考基准面上,这种消除高程影响的校正称为高程校正。
剩余静校正
地表一致性剩余静校正是在应用了前面介绍的野外静校正或者折射静校正以后进行的。
由于多种因素,一个CMP道集中的各个地震道,经过上面的静校正之后,仍然存在着剩余静校正量。
一类是基于地表一致性时差分解的方法;
另一类是基于互相关(或称叠加能量最大)的剩余静校正方法。
1.最大叠加能量法剩余静校正
2.相邻叠加道相关法静校正
见书P124-128
叠加速度普原理
见书P133-134
速度谱用途
速度谱的速度分析是主要的速度分析手段,准确可靠的速度分析是地震数据处理的基础,速度参数不仅关系到地震数据处理诸多环节的质量,其本身也提供了关于地下构造和岩性的重要环节。
为水平叠加、偏移等处理提供处理参数。
叠加速度普含义:
叠加公式求的是平均振幅能量就是平均振幅能量
叠加公式求的是平均振幅就是平均振幅
平均速度:
均方根速度:
从第一层到第n层的均方根速度为
狄克斯(Dix)公式;
利用DIX公式不但可以有均方根速度(或叠加速度)求层速度,也可以求平均速度。
偏移
为什么要进行偏移归位?
(1)从几何反射的观点讨论
(2)从广义反射的观点讨论
(3)从波长分析的角度讨论
射线法偏移思想
叠加剖面的同相轴可以看作由若干脉冲振幅组成,通常将输入剖面中的一个脉冲振幅,经过偏移后得到的地下空间图形称为输入剖面的偏移脉冲响应。
如图7一4(a)所示,假设在水平坐标处有一地震记录,其t时刻有一脉冲振幅,波的传播速度是一个常数,则偏移脉冲响应的轨迹如图7-4(b)所示,是以点为圆心,以vt/2为半径的一个半圆,半圆轨迹上的振幅与输入脉冲振幅成正比。
也就是说,如果地下反射界面是图7-4(b)所示的一个半圆,在点作自激自收观测,则来自半圆界面上的反射波将汇聚成一个脉冲出现在地震记录上。
叠后偏移——叠前部分时间偏移(DMO)——叠前时间偏移(Kirchhoff积分偏移和波动方程偏移)——叠前深度偏移(Kirchhoff积分偏移和波动方程偏移)
炮域波动方程叠前深度偏移
不同偏移方法的优缺点
I绕射扫描叠加偏移
优点:
不受倾角大小限制、计算量小
缺点:
不适用于横向变速介质、几何偏移,动力学特点受到畸变,偏移噪声比较大。
IIFK偏移
适用于大倾角,效率高,运行时间小(FFT)
产生空间假频,偏移噪声比较大,不适用于横向变速介质。
ⅢKirchhoff积分
适用于大倾角,效率高
不适用于横向变速介质,偏移噪声比较大。
Ⅳ有限差分法波动方程偏移
偏移的倾角达到15°
、45°
、60°
等,偏移噪声小,抗干扰强,在信噪比比较低时也有效,能适应于横向