地震勘探资料处理第一章第六章复习要点总结第一章 地震数据处理基础一维谱分析 数字地震记录中，每个地震道是一个按一定时间采样间隔排列的时间序列，每一个地震 道都可以用一系列具有不同频率、不同振幅、相位的简谐曲线叠加而成。应用一维傅里叶变换可以得到地震道的各个简谐成分； 应用一维傅里叶反变换可以将各个简谐成分合并为原来的地震道序列。 连续函数正反变换公式：gX () = J x(t)e-iotdt正变换gx(t)=丄 J X (o)eiotdo2兀g反变换通常由傅里叶变换得到的频谱为一个复函数，称为复数谱。它可以写成指数形式X （）=1 X （）I e（3）= A（）ei（）式中A（3）为复数的模，称为振幅谱；申（3）为复数的幅角，称为相位谱。A（3）=#阿战阿，（3）=tan-1（弧度也可换算为角度）r离散情况下和这个差不多（看PPT和书P2-3）一维傅里叶变换频谱特征：1、一维傅里叶变换的几个基本性质（推导）线性 翻转 共轭 时移 褶积 相关（功率谱）， P3-72、Z 变换（推导）13、采样定理假频尼奎斯特频率，f厂N 2At二维谱分析二维傅里叶变换X （o, k）称为二维函数X （x, t）的频波谱。其模量I X （3, k）1称为函 数X（x,t）的振幅谱。由X（,k）这些频率f与波数k的简谐成分叠加即可恢复原来的波场 函数X （x, t）（二维傅里叶反变换）。如果有效波和干扰波的在f-k平面上有差异，就可以利用二维频率一波数域滤波将它们分开，达到压制干扰波，提高性噪比的目的。二维频谱产生空间假频的原因数字滤波在地震勘探中，用数字仪器记录地震波时，为了保持更多的波的特征，通常利用宽频带 进行记录，因此在宽频带范围内记录了各种反射波的同时，也记录了各种干扰波。有效波和 干扰波的差异表现在多个方面（频谱、传播方向、能量）。利用频谱特征的不同来压制 干扰波，以突出有效波的方法就是数字滤波。滤波器的频率特性：滤波器时间函数或滤波因子h（t）的频谱方仙）称为滤波器的频率特性。 滤波器的时间特性（单位脉冲响应）：在时间域的表示方法中，令一个单位脉冲通过一个滤 波器，然后观察滤波器的输出，这个滤波器输出的自然过程曲线称为滤波器的时间特性，也 称滤波器的单位脉冲响应。一维频率域滤波机理:啟）!图频率域滤波原理图文（）二方（）文（）（对输出进行傅里叶反变换）输出信号的振幅谱等于输入信号的振幅谱与滤波器的振幅频率特性的乘积, 输出信号的相位谱等于输入信号的相位谱与滤波器的相位谱之和。I X （）1=1方（）1 -1文（）1，九（）= （）+ （）xxH对H（w ）的要求（看PPT和书P13）一维频率域滤波的实现步骤(1) 对已知地震记录道进行频谱分析，确定有效波和干扰波频带范围。(2) 设计合适的滤波器，要能压制干扰波，突出有效波。(3) 进行频率域滤波运算，文()二方()文()，计算输出信号的频谱X(o)。(4) 对输出信号的频谱X(o )进行傅里叶反变换，求滤波后的地震记录X(t)。最小相位、最大相位、混合相位离散傅里叶变换(DFT)滤波的特殊问题：周期性、对称性。由连续傅里叶变换过渡到离散傅里叶变换利用：AtAf = N一维时间域滤波机理：时间域滤波就是褶积滤波。x(t) = x(t) * h(t)褶积滤波的物理意义：它相当于把地震信息x(t)分解为起始时间、极性、幅度各不相同的脉冲序列，令这些脉冲序列按时间顺序依次通过滤波器，这样在滤波器的输出端就得到对输入脉冲序列的脉冲响应，这些脉冲响应有不同的起始时间，不同的极性和不同的幅度(这个 幅度是与引起它的输入脉冲幅度成正比的)，将它们叠加起来就得到滤波后的输出X(t)。对h(t)的要求(看PPT和书P22)一维时间域滤波(褶积滤波)的实现步骤(1) 对地震记录进行频谱分析，确定通频带中心频率f0和带宽2Af。(2) 确定滤波因子长度N, N通常为101。(3) 通过傅里叶反变换求滤波因子h(nAt)。(4) 通过时间域褶积计算：X(t) = x(t) * h(t)，在计算机上实现时要将其离散化。数字滤波的两个特殊性：1、数字滤波的离散性造成伪门现象；2、数字滤波的有限性造成吉卜斯现象。伪门产生的过程和克服方法：对连续的滤波因子h(t)用时间采样间隔At离散采样后，得到 h(nAt)。如果，再按h(nAt)计算出与它相应的滤波器的频率特性，这时在频率特性的图形 上,除了有同原来的H(o)对应的“门”外，还会周期性地重复出现很多个“门”，这些门 称为“伪门”。*j t _ L,i . i n i T| t | i r i八/V1IMIJ裁圈i at伪门jtlmI为了避免干扰波通过“伪门”而对有效波造成影响，可以适当地选取采样间隔At，使第一 个“伪门”出现在干扰波的频谱范围之外。吉卜斯现象产生的过程和克服方法有限长度的h（t）对应的方（f）不再是一个门式滤波，而是有波动的曲线（图1-31），曲线由间断点向远处波动衰减，在间断点波动最大，这种现象叫作非连续函数频率响应的吉卜 斯现象。为了避免吉卜斯现象对有效波造成畸变的影响，可以采用镶边法，即在频率特性曲线的 不连续点附近，镶上一条连续的边，这种做法克服了频率特性曲线的波动问题（吉卜斯现象）。03 L-3L非违续罄数的吉k期现鏡滤波器类型：低通、高通、带通（时域和频域） 1,1 l Aw高通：H（w）=1,1 w | Aw0,1 w | Aw带通：低通和高通中间。第二章预处理及真振幅恢复时序记录：按照采样时间顺序排列的多路传输地震记录，称为时序记录。 道序记录：以地震道为顺序排列的记录，称为道序记录。数据解编目的：通过数据解编将地震数据的记录顺序由时序转化为道序，它相当于对地震数 据进行矩阵转置。真振幅恢复的目的和主要内容：其目的是尽量对地震波能量的衰减和畸变进行补偿和校正, 主要处理内容包括波前扩散能量补偿、地层吸收能量补偿和地表一致性能量调整。vv波前扩散因子：D(均匀介质)，D二(层状介质)d V 2td V 2 trms品质因子： 实际地震资料处理中常用品质因子Q来描述地震波的衰减，其意义是地震波在传播一个波 长九距离后，原来储存的能量E与所消耗能量AE之比，即A20A2 A20 尢品质因子近似表达式为：Q二沁吸收衰减因子：-ftD 二 ear 二 ePr 二 e qa 可见在非完全弹性介质中，地震波的高频成分比低频成分衰减得要快。地层吸收和波前扩散对地震波振幅影响的区别：(为什么说大地是一个低通滤波器，地层吸 收具有低通滤波作用)地震波振幅的衰减与频率有关，频率越高，振幅衰减越严重。地层吸收不仅造成地震波振幅 的衰减，而且对地震波产生低通滤波作用。第三章反褶积1、什么是反褶积，其主要作用是什么？反褶积定义：反褶积处理是褶积处理的反过程，因而称为反褶积。具体说来如果我们设计一 个滤波器，其滤波因子w(t)具有与滤波器w(t)有恰好相反的性质，即当输入为地震道记录 x(t)时，其输出为地层脉冲响应e(t)，称这个反过程为反滤波或反褶积。 Q图3亠2图3-1(b)的反榕 枳滤波过捲反褶积的作用：其主要作用是压缩地震记录中的地震子波，提高地震的垂直分辨率；同时， 可以压制海上鸣震和多次波，提高地震资料信噪比(利用预测反滤波)。2、求取地震子波的方法主要有哪些，其优缺点及适用条件是什么？P62-69地震子波：由震源发出的尖脉冲经过大地滤波器的滤波作用后，变成一个具有一定时间延续度的波形w(t)，该波形通常叫作地震子波(1) 直接观测法这种方法是用专门布置在震源附近的检波器直接记录地震子波w(t)，此方法只适用于 海上地震勘探(海水含盐量不同，具有明显的分层性)。(2) 自相关法掌握：如果地震子波w(t)是零相位的，且反射系数r(t)为白噪声，则w(t) = - f W()eietd2兀eitd=丄f|文)|2兀_兀(3) 利用测井资料求子波的方法 这种方法要求有良好的声波测井和密度测井资料，并且在井旁有质量较高的地震记录 但该方法不必预先知道地震子波的相位特性，也不必假设反射系数是白噪声。(4) 对数分解法 这种方法也不需要假设地震子波是最小相位和反射系数是白噪声，但对地震记录质量要 求较高，用对数谱序列平均法求取地震子波效果较好。实际地震记录模型：除了一系列反射波S(t) = w(t) * r(t)而外，还存在着干扰波n(t)，因此，地震记录x(t) 的一般模型可以写为：x(t) = S(t) + n(t) = w(t) * r(t) + n(t) = w(T)r(t -t) + n(t)T=03、试推导求最小平方滤波因子或反滤波因子的最小平方滤波方程？ P70-75 或试推导求最小平方反滤波因子的最小平方反滤波方程？最小平方滤波概念：它是在滤波器实际输出与期望输出的误差平方和为最小的情况下，确定 滤波器的滤波因子的，因而称为最小平方滤波。最小平方滤波两个假设条件:(1)反射系数序列r(t)为随机白噪序列。2)随机噪声干扰n(t) 为白噪声。托布里兹(Toeplitz)矩阵：它不但以主对角线为对称，也以次对角线为对称，而且主对角线 及与主对角线平行的直线上的元素均相同最小平方滤波方程：r (0)r bbbbr (1)r (0)bbbbr (m)bbr (m -1)bba(0)a(1)(r (0)dbr (1)dbr (m)r (m -1)bbbbr (0)人a(m)丿bbr (m)丿db最小平方反滤波方程(不考虑噪声时)：(r (0)wwr (1)wwr (1)wwr (0)wwr (m)wwr (m -1)wwa(0)a(1)(r (0)dwr (1)dw、r (m)r (m -1)wwwwr (0)人a(m)丿wwI rdw (m)丿4、简述预测反滤波的基本原理：预测滤波就是要设计一个预测因子 c(t) ， 对 输入 x(t) 已知的过去值 x(t m), x(t m +1),x(t 2), x(t 1)和现在值x(t)进行滤波处理，获得未来某个时刻 t+Q时的预测值(海上鸣震和多次波干扰)X(t + a) = c (t) * x(t) = c (t )x(t T)T=0通过这种预测的方法，根据地震记录一次反射波和干扰波的信息预测出纯干扰波部分，再从 包括一次波和干扰波的地震记录X(t + a)中减去纯干扰波部分X(t + a)，就得到消除干扰后 的一次反射波信号S(t + )，以消除一次反射波后面的海上鸣震和多次波干扰。S (t + a) = x(t + a) x (t + a)预测反滤波参数选取方法及其对预测反滤波结果的影响(P88-90)( 1 )预测步长 a( 2)预测因子长度 m(3) 预白化量第四章 动校正和