《2022年地震资料处理复习总结》由会员分享,可在线阅读,更多相关《2022年地震资料处理复习总结(10页珍藏版)》请在金锄头文库上搜索。
1、学习必备欢迎下载一维谱分析数字地震记录中,每个地震到是按一个按一定时间采样间隔排列的时间序列数字滤波,每一个地震道都可以用一系列具有不同频率和不同振幅、相位的简谐曲线叠加而成。应用一维傅里叶变换可以得到每个地震道德简谐成分;应用傅里叶反变换可以将简谐成分合成为原来的地震道的时间序列函数。通常由傅里叶变换得到的频谱为一个复函数,称为复数谱。它可以写成指数形式式中复数的模,称为振幅谱; 复数的幅角,称为相位谱。离散情况下和这个差不多一维频谱的特征:1 傅里叶变换的几个基本性质线性 翻转共轭时移褶积相关(功率谱)2 假频 尼奎斯特频率二维谱分析二维波场函数X(x,t)的二维傅里叶变换X(,)表明了二
2、维波场函数X(x,t)的各个频率f一波数简谐成分的频一波谱。由X(,)这些频率f 一波数的简谐成分叠加即可恢复原来的波场函数X(x,t)。二维傅里叶变换X(w,k)称为二维函数X(x,t)的频一波谱。其模量为函数 X(x,t)的振幅谱。如果有效波和干扰波得平面简谐波成分有差异,有效波的平面简谐波成分与干扰波的平面简谐波成分不同的视速度传播,则可以用二维视速度滤波将他们分开,达到压制干扰, 提高性噪比的目的。二维频谱的特征:空间假频()()()()( )iwiwX wX w eA w e)(A()1()()tan()irx wwxw22( )( )( )riA wxwxwtf21NokokX(,
3、 )k精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 1 页,共 10 页学习必备欢迎下载数字滤波在地震勘探中,用数字仪器记录地震波时,为了保持更多的波得特征,通常利用宽频带进行记录, 因此在宽频带范围内记录了各种反射波的同时,也记录了各种干扰波。有效波和干扰波得差异表现在多个方面(频谱、传播方向、能量)。利用频谱特征的不同来压制干扰波,以突出有效波的方法就是数字滤波。滤波器的响应特性:对滤波器能力的最普遍度量是其响应特性滤波器的频率特性: 其滤波器时间函数或滤波因子的频谱称为滤波器的频率特性,滤波器的时间特性(单位脉冲响应):在时间域的表示方法中
4、,令一个单位脉冲通过一个滤波器, 然后观侧滤波器的输出,这个滤波器输出的自然过程曲线称为滤波器的时间特性。也称滤波器的“脉冲响应”频率响应函数应该就是时间和频率响应函数合起来应该就是就是响应特征滤波机理:输出信号的振幅谱等于输入信号的振幅谱与滤波器的振幅频率特性的乘积,输出信号的相位谱等于输入信号的相位潜与滤波器相位特性之和。(频率)时间域上就是褶积褶积滤波的物理意义:它相当于把地震信息分解为起始时间、极性、 幅度各不相同的脉冲序列, 令这些脉冲按时间顺序依次通过滤波器,这样在滤波器的输出端就得到对输入脉冲序列的脉冲响应,这些脉冲响应有不同的起始时间,不同的极性和不同的幅度(这个幅度是与引起它
5、的输入脉冲幅度成正比的),将它们叠加起来就得到滤波后的输出分。频率域滤波的步骤(1)对已知地震记录道进行频谱分析。(2)设计合适的滤波器(3)进行滤波运算(4)对输出信号谱进行傅里叶反变换褶积滤波的具体计算褶积滤波的具体计算步骤如下:(1)对地震记录进行频谱分析,确定通频带中心频率和带宽。(2)确定滤波因子长度N。)()()(wHwXwX)( tx)(wH)(th)(wH)()()(wwwHxx)()()(wHwXwX)(wX)(tx精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 2 页,共 10 页学习必备欢迎下载(3)求滤波因子。(4)将式写成
6、离散形式数字滤波的 离散型 造成伪门现象,数字滤波的 有限性 造成吉布斯现象。伪门: 对连续的滤波因子用时间采样间隔离散采样后,得到。如果,再按计算出与它相应的滤波器的频率特性,这时在频率特性的图形上,除了有同原来的对应的 “ 门” 外,还会周期性地重复出现很多“ 门” ,这些门称为 “ 伪门 ” 。为了避免“伪门”造成的影响,可以适当地采样间隔使第一个“伪门”出现在干扰波的频谱范围之外。吉卜斯现象有限长度的( )h t对应的( )Hf不再是一个门式滤波,而是有波动的曲线(图 1-31),曲线由间断点向远处波动衰减, 在间断点波动最大, 这种现象叫作非连续函数频率响应的吉卜斯现象。实际计算中脉
7、冲响应函数只能取有限长,截断后的脉冲响应所对应的频率响应函数不再是一个理想的“门” ,而是接近于这个门得一条幅值有波动的光滑曲线。镶边法解决。dtxytx)()()()()()(2121mhmtnxtnxNNm)(wH)(tht)(tnht精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 3 页,共 10 页学习必备欢迎下载类型:低、带、高通低通:高通1,0,()wwwwH w带通:低通和高通中间。滤波器的分类:按滤波器的性质(即响应函数)划分1) 无畸变滤波器。振幅特性为常数,相位特性是线性的滤波器。这种滤波器不改变输入信号的波形,频率响应为00(
8、)itHa e0a、0t均为常数00( )()y ta x tt00()()itYa Xe2) 相位畸变滤波器(纯相位滤波器,全啊通滤波器)只改变输入信号的相位谱,振幅形状不改变。振幅特性常数0|( ) |Ha,但相位特性不是线性。3) 振幅畸变滤波器振幅特性不是常数,实际工作中希望滤波时不是信号产生相位畸变或相位移,为零相位滤波器。滤波器的稳定性:当输入信号为有限,其输出信号也为有限,这种滤波器就是稳定的。充要条件:|( ) |h tdt物理可实现性:满足因果律(输入之前不产生输出)的滤波器称为物理可实现的。充要条件:( )0h t当 t0 时对于 z 变换为多项式的滤波器来说,对于 z 变
9、换为有理分式的滤波器来说,当所有的根均不在单位圆(|z|=1 )上时。这个滤波器是稳定的,当所有根都在圆外时,这个滤波器是物理可实现的。反滤波作用反褶积的 基本作用 是压缩地震记录中的地震子波,同时, 可以压制鸣震和多次波,因而反褶积可以明显提高地震的垂直分辨率。反褶积处理是褶积处理的反过程,因而称为反褶积。如果我们设计一个滤波器,其滤波因子( )w t具有与滤波器( )w t有恰好相反的性质,即当输入为地震道记录( )x t约时,其输出为地层脉冲响应( )e t。 称这个反过程为反滤波或反褶积。由震源发出的尖脉冲经过大地滤波器的滤波作用后,变成一个具有一定时间延续的波形( )w t,通常叫作
10、 地震子波ww, 1)(wHww, 0精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 4 页,共 10 页学习必备欢迎下载实际地震记录:除了所表示的一系列反射波( )S t而外,还存在着干扰波( )n t,因此,地震记录( )x t的一般模型可以写为0( )( )( )( ) ()( )x tS tn twr tn t最小平方反滤波两个假设条件:1.反射系数为白噪声2.假设随机噪声干扰( )n t为白噪声子波未知情况下的假设条件:1.反射系数为随机的白噪序列2.假设地震子波是最小相位(预测误差滤波也是这两个假设条件)预测反滤波预 测 滤 波 就 是
11、 要 设 计 一 个 预 测 因 子( )c t, 对 输 入( )x t已 知 的 过 去 值(), (1), , (2), (1)x tm x tmx tx tL和现在值( )x t进行滤波处理,获得未来某个时刻ta时的预测值0?()( )( )( ) ()mx tac tx tcx t用预测的方法, 根据地震记录一次反射和干扰的信息预测出纯干扰部分,再由包括一次波和干扰的地震记录中减去纯干扰部分,得到消除干扰后的一次反射信号,以消除一次反射后面的海上鸣震等多次波干扰。叠加的目的是压制干扰,提高地震数据的信噪比。动校正的目的是消除炮检距对反射波旅行时的影响动校正 :将不同炮检距的反射时间校
12、正到零炮检距反射时间的过程称为动校正22220( )( )(0)(0)( )(0)2NMOxxt xt xtttvVt动校正量随炮检距递增,随反射深度和速度递减动校正拉伸 : 数字动校正造成的波形拉伸称为动校正拉伸( )r t( )r t212211(T)()12()T精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 5 页,共 10 页学习必备欢迎下载由于浅层的动校正时差大于深层的动校正时差,所以浅层大炮检距的拉伸最为严重。校正后反射波的形状不同,发生畸变,畸变的一般规律是反射波得波形被拉长,周期加大,频谱向低频方向移动。拉伸系数 : 0mNMOm
13、ftTTft反射深度越浅,炮检距越大( 越大 ),动校正拉伸越严重,子波的主频向低频转移也随之严重。静校正 : 静校正是校正以及消除由于地表高程和地下低、降速带变化对反射波旅行时的影响。野外一次静校正又称为基准面静校正,包括井深校正,地形校正,低速带校正。人为选定一个海拔高程作为基准面,利用野外实测得到的各点高程、低速带厚度、 速度、或井口时间等资料,将所有的炮点和检波点都校正到此面上,用基岩速度代替低速带速度,从而去掉表层因素的影响。定义 为:用于补偿由于地表高程变化、风化层的厚度和速度变化对地震资料的影响。其目的是获得在一个平面上进行采集,且没有风化层或低速介质存在时的反射波到达时间。静校
14、正的目的是获得在一个平面上进行采集,且没有风化层或低速介质存在时的反射波到达时间 . 静校正信息来自于两个方面: 一是野外测量和观测的数据,包括地面高程数据、井口检波器记录时间、微测井和小折射数据等;称为基准面校止或野外静校正. 二是根据初至波时间和地下反射信息求取静校正量,称为初至折射静校正和反射波地表一致性剩余静校正。地震数据被校止到参考基准面 上,消除了地表起伏和风化层横向变化的影响,后续地震处理工作就好像地震数据是在基准而上采集的。有时我们不是将地震数据一次校正到参考基准面或最终基准面上,而是首先将地震数据校正到一个中间基准面上,这个基准面有时也称为浮动基准面或CMP 叠加基准面 。速
15、度分析、剩余静校正、动校正、叠加都在这个基准面上进行。叠加之后,再将地震数据由浮动基准面校正到参考基准面或最终基准面上,在最终基准面上完成偏移处理。静校正概念中“ 静” 的含义TTTTTTTNMOt精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 6 页,共 10 页学习必备欢迎下载是相对动校正中“ 动” 的含义而言的, 我们知道, 地震道的动校正时差是反射时间的函数,而地震道的静校正时差与地震道的时间无关,无论是浅层反射,还是深层反射, 整个地震道只有一个静校正量。基准面静校正也称为野外静校正,顾名思义, 就是将在地表采集的地震记录校正到基准面上,
16、消除地表高程和风化层对地震记录旅行时的影响。第一步是剥去风化层的影响第二步, 再将地震记录时间由风化层的底界校正到参考基准面上,这种消除高程影响的校正称为高程校正。剩余静校正地表一致性剩余静校正是在应用了前面介绍的野外静校正或者折射静校正以后进行的。由于多种因素,一个CMP 道集中的各个地震道,经过上面的静校正之后,仍然存在着剩余静校正量。一类是基于地表一致性时差分解的方法;另一类是基于互相关(或称叠加能量最大)的剩余静校正方法。1.最大叠加能量法剩余静校正2.相邻叠加道相关法静校正见书 P124-128 叠加速度普原理见书 P133-134 速度谱用途速度谱的速度分析是主要的速度分析手段,准
17、确可靠的速度分析是地震数据处理的基础,速度参数不仅关系到地震数据处理诸多环节的质量,其本身也提供了关于地下构造和岩性的重要环节。为水平叠加、偏移等处理提供处理参数。叠加速度普含义:叠加公式求的是平均振幅能量就是平均振幅能量叠加公式求的是平均振幅就是平均振幅平均速度:11( )nkkkankkv tVvtv精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 7 页,共 10 页学习必备欢迎下载均方根速度:从第一层到第n 层的均方根速度为22110,1nnk kk kkkrmsnnkkv tv tVtt狄克斯 (Dix)公式 ; 22,0,10,120,0
18、,1()()rms nnrms nnnnnvtvtvtt利用 DIX公式不但可以有均方根速度(或叠加速度)求层速度nv,也可以求平均速度。偏移为什么要进行偏移归位?(1)从几何反射的观点讨论(2)从广义反射的观点讨论(3)从波长分析的角度讨论射线法偏移思想叠加剖面的同相轴可以看作由若干脉冲振幅组成,通常将输入剖面中的一个脉冲振幅,经过偏移后得到的地下空间图形称为输入剖面的偏移脉冲响应。如图 7 一 4(a)所示, 假设在水平坐标dx处有一地震记录,其 t 时刻有一脉冲振幅,波的传播速度是一个常数,则偏移脉冲响应的轨迹如图7-4(b)所示,是以dx点为圆心,以vt/2为半径的一个半圆,半圆轨迹上
19、的振幅与输入脉冲振幅成正比。也就是说,如果地下反射界面是图7-4(b)所示的一个半圆,在dx点作自激自收观测,则来自半圆界面上的反射波将汇聚成一个脉冲出现在地震记录上。叠后偏移 叠前部分时间偏移(DMO)叠前时间偏移 (Kirchhoff积分偏移和波动方程偏移) 叠前深度偏移(Kirchhoff积分偏移和波动方程偏移)精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 8 页,共 10 页学习必备欢迎下载炮域波动方程叠前深度偏移不同偏移方法的优缺点I 绕射扫描叠加偏移优点:不受倾角大小限制、计算量小缺点:不适用于横向变速介质、几何偏移,动力学特点受到畸
20、变,偏移噪声比较大。II FK偏移优点:适用于大倾角,效率高,运行时间小(FFT )缺点:产生空间假频,偏移噪声比较大,不适用于横向变速介质。 Kirchhoff积分优点:适用于大倾角,效率高缺点:不适用于横向变速介质,偏移噪声比较大。 有限差分法波动方程偏移优点:偏移的倾角达到15 、45 、60 等,偏移噪声小,抗干扰强,在信噪比比较低时也有效,能适应于横向变速缺点:计算量大道内平衡基本思想:将一道上的能量强的波减弱,能量弱的波增强。所使用的方法是加权,即能量强的波乘以一个相对较小的权值, 而能量弱的波乘以一个相对较大的权值。权值的确定由计算机自动进行。道间平衡基本思想:与道内平衡相似,区
21、别:其对象为地震道而不是一个单纯的反射波。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 9 页,共 10 页学习必备欢迎下载采集Acquisitim 处理processing 解释interpretation 信噪比Signal/Noise 分辨率Revolution 参数parameters 褶积convolution 反褶积De convolution 时间域time-domain 带通Band-pass 低通Low-pass 高通High-pass 子波wavelet地震道记录trace 采样间隔intervade 正演Forward Modeling 反演inversion (Winener )least-squares-deconvolution Toepliz matrix 海上鸣震干扰Water reverberation 正常时差Normal movement 低速带LVL low velocity level (layer)倾角时差校正DipMoveOut CMP 共中心点反射CDP 共深度反射CRP 共反射点反射DMO 叠前部分时间偏移精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 10 页,共 10 页