地震数据处理中关于地震子波相位特性的探讨
地震子波
盲解卷积算法1. 原理：几个重要概念：1.1 褶积模型假设：假设1：地层是由具有常速的水平层组成；假设2：震源产生一个平面压缩波(P 波)，法向入射到层边界上，在这种情况下，不产生剪切波(S 波)；假设3：震源波形在地下传播过程中不变，即它是稳定的； 假设4：噪音成分是零；假设5：震源波形是已知的；假设6：反射系数序列是一个随机过程。
这意味着地震记录具有地； 震子波的特征，即它们的自相关和振幅谱是相似的；假设7：地震子波是最小相位的，因此，它有一个最小相位的逆。
1.2 反滤波如果定义滤波算子为f (t),则f (t) 与已知地震记录x(t)的褶积得到一个对地层脉冲响应e(t)的估计e(t) = f (t)∗ x(t)； （1） x(t) = w(t)* f (t)* x(t)； （2） δ(t) = w(t)* f (t)； （3）1()()*()f t t w t δ= （4） 用流程图表示为：1.3 震源反子波计算震源反子波在数学上是利用z 变换来实现的。
例如，假设基本子波为两点时间序列（1，－0.5)1()12w z z =-2111()1 (12412)F z z z z ==+++-； （5） 2111()1 (12412)F z z z z ==+++- （6） ()F z 的系数11(1,,,...)24代表逆滤波算子f(t)有关的时间序列。
可以看出它有无限多个系数，然而它们递减的很快。
如同任何滤波过程一样，实际应用的算子都是被截断的。
1.4 最小平方反滤波当输入子波良好，其z 变换的逆可以用一收敛序列表示，则上面所描述的反滤波将得到一个很好的近似于尖脉冲的输出.将下面问题列出方程式：给定输入子波(1－0.5)，寻找一个二项滤波器(a,b)，使实际输出与期望输出(1,0,0)之间的误差在最小二乘的意义上最小。
将滤波器(a,b)与输入子波(1,-0.5)褶积以计算实际输出。
误差的累积能量L 定义为实际的与期望的输出的系数之差的平方和：222(1)()()22a b L a b =-+-+- （7） 目的是寻找系数(a,b)使L 最小，这要求变量L 随系数(a,b)而变并使之为零对上式进行简化，取L 对a 和b 的偏导数，并使其为零，得到： 522502a b a b ⎧-=⎪⎪⎨⎪-+=⎪⎩（8） 有两个方程和两个未知数即滤波器系数(a,b)，可变成下列普通矩阵形式：5/21215/20a b -⎛⎫⎛⎫⎛⎫= ⎪⎪ ⎪-⎝⎭⎝⎭⎝⎭（9） 求解滤波系数，我们得到(a,b)为(0.95,0.38)。
《地震资料数字处理》复习
《地震资料数字处理》复习地震资料数字处理围绕以下三方面工作：1、提高信噪比；2、提高分辨率；3、提高保真度。
一、提高信噪比的处理1、原理利用噪声和信号在时间、空间、频率和其他变换域中的分布差异，设计滤波因子，将噪声进行压制。
2、处理顺序提高信噪比包含消除噪声和增强信号两部分内容。
消除噪声一般在叠前的各种道集上进行，主要针对规则干扰如多次波和面波等，增强信号一般在叠后剖面上进行，主要针对随机噪声。
3、随机噪声是指没有固定的频率、时间、方向的振幅扰动和震动，其成因大致是来自环境因素、次生因素和仪器因素，其中次生干扰的强度与激发能量有关。
随机噪声在记录上表现为杂乱无章的波形或脉冲，在频率上分布宽而不定，在空间上没有确定的视速度。
随机噪声的随机性与道间距有关，如果道间距减小到一定程度，许多随机噪声表现出道间的相干性，当道距大于随机噪声的相干半径才表现出随机性。
4、一维滤波器（伪门、Gibbs现象）频率滤波器是根据信号和噪声在频率分布上的差异而设计时域或频域一维滤波算子。
它压制通放带以外的频率成分，保留通放带以内的频率成分。
Gibbs现象是由于频率域的不连续或截断误差引起的，通放带和压制带之间设置过渡带可克服此现象，设计滤波器就是控制过度带的形状和宽度。
5、二维滤波器二维滤波是根据有效信号和相干噪声在视速度分布上的差异，来压制噪声或增强信号。
通常用来压制低视速度相干噪声，在f-k平面上占据低频高波数区域。
二维滤波比较容易产生蚯蚓化现象，而且混波相现象明显，在空间采样条件不满足或陡倾角的情况下受到空间假频的影响，一般常用于压制一些规则干扰，如面波和多次波等。
6、频率-波数域二维滤波实现步骤：（1）把时间和空间窗口里的数据变换到f-k域；（2）在f-k域，通过外科切除，按径向扇形划分压制区C(乘振幅置零)、过渡区S（乘振幅置0至1变化）、通放区P (乘振幅置1) ；（3）从f-k域反变换到t-x域。
8、数字滤波有两个特殊性质：（1）数字滤波由于时域离散化会带来伪门现象，（2）由于频域截断会造成吉卜斯现象。
地震子波 数字信号实验报告
根据所学知识对实验结果进行分析；
地震子波由震源激发，在地层中传播，因为在沉积地层中，每层介质的物理性质不相同，从而使得地震波的传播速度也不相同。当地震波传播到两层介质的分界面时，会发生反射，由于每层介质的反射系数不同，所以反射波的能量也不相同，检波器接收到不同时刻的、不同能量的反射波，形成一个地震记录。
subplot(2,2,4),plot(f2,'k'),title('最小相位子波地震记录')
f1(i)=0;
f2(i)=0;
end
fori=1:n+m-1%地震记录长度
forj=1:m%地层深度
ifi-j>0&i-j<=200%满足雷克子波的时间序列长度
a1(j)=r(j)*w1(i-j);
a2(j)=r(j)*w2(i-j);
f1(i)=f1(i)+a1(j);
f2(i)=f2(i)+a2(j);
地震子波波形显示及一维地震合成记录
姓名：杨肖迪学号：050422009040专业：地球信息科学与技术2009级
实验目的
1.认识地震子波（以雷克子波为例），对子波的波形有直观的认识。
2.利用褶积公式合成一维地震记录。
实验步骤
1.雷克子波
（零相位子波）
（最小相位子波）
其中 代表子波的中心频率， 代表子波宽度,随着 的增大，子波能量后移，当 =7时，最小相位子波可视为混合相位子波，这里 = 25 Hz， =4;
end
end
end
subplot(2,2,1),plot(w1,'k'),axis([0,200,-1,1]),title('零相位子波')
子波的零相位,最小相位,混合相位
子波的零相位,最小相位,混合相位子波是地震波在地下传播时激发的复杂波形，在地震勘探中起到了重要的作用。
在地震勘探中，子波的信息可以用来了解地下结构，进行地震处理和解释。
子波有三种形式：零相位子波、最小相位子波和混合相位子波。
接下来，我们将详细介绍这三种子波的特点和应用。
首先是零相位子波。
零相位子波又称为纯实信号，它的相位谱是一个常数，即各个频率分量的相位都是零。
这意味着零相位子波是非带通的，其频谱范围是从直流到无限大。
零相位子波是地震波的极限情况，一般在理论分析中使用。
它具有宽频带、短时间的特点，可以用来表示瞬态信号。
最小相位子波是一种具有最小相位谱的子波。
最小相位谱指的是在给定幅度谱相同的情况下，具有最小相位谱的子波。
最小相位子波具有狭窄的频带、宽延迟和长时间的特点，可以用来近似地表示地震记录中的真实地震波形。
最小相位子波是地震记录的有效信号。
混合相位子波是介于零相位子波和最小相位子波之间的一种子波形式。
混合相位子波的相位谱介于零相位子波和最小相位子波的相位谱之间。
混合相位子波既有零相位子波的宽频带特性，又有最小相位子波的稳定性。
因此，混合相位子波在地震勘探中被广泛应用于地震处理和解释。
通过对地震记录进行混合相位子波叠加，可以提高地震信息的分辨率和解释能力。
在地震勘探中，选择合适的子波形式对地震数据进行处理和解释非常重要。
零相位子波可以用来分析地下结构的反射特性，最小相位子波可以用来近似地还原地震记录中的原始地震波形，混合相位子波则结合了两者的优点，可以提高地震处理和解释的能力。
为了更好地了解子波的特点和应用，我们需要对子波的频谱、振幅、相位进行分析。
在频域上，子波的频谱决定了子波的带宽，即频率范围；在时域上，子波的振幅和相位决定了子波的波形。
通过分析子波的频谱、振幅和相位，我们可以把子波应用到地震数据的处理和解释中。
子波的零相位、最小相位和混合相位是地震波形分析中的重要概念。
零相位子波是理论分析的基础，最小相位子波可以近似地还原地震记录中的地震波形，混合相位子波结合了两者的优点。
地震子波提取方法论文
地震子波提取方法论文摘要：在信噪比较高的情况下最小相位和混合子波相位均适应性良好。
在信噪比降低后，混合相位提取子波受影响较大。
通过对实际模型的处理，表明了该方法的有效性和实用性。
1 子波基本理论子波包括最小相位子波、最大相位子波、混合相位子波。
子波的Z变换是一个多项式。
若此多项式的全部零点均在单位圆外，则为最小相位子波；在单位圆内为最大相位子波，零点在单位圆的内外都有，则为混合相位子波。
2 复赛谱域提取子波法复赛谱域提取子波法研究的是一种非线性系统，或者说是一种线性滤波的推广。
对于一个线性系统，它可以满足以下迭加原理：T[ax1+bx2]=aT[x1]+bT[x2]（1）式中a、b为常数，T为表示函数。
一个信号是由多个信号相加合成的，利用线性滤波的方法，就能有效地把它们分开。
如果一个信号不是由多个信号相加，而是由多个信号中褶积合成的，那么，用这种线性滤波的方法就不能将其分开。
对于这种信号的处理，我们仍可利用（2）式的基本思想，把它推广到褶积合成的信号。
为此我们需要一个称为同态系统的变换系统D，它具有把一个褶积关系转变为相加关系的特征系统。
从（2）式可见，线性系统特别适用于分开相加合成的信号；也就是说，一个系统可以这样处理：D[（a）x1*（b）x2]=aD[x1]+bD[x2]（2）上式中（a）表示标量乘，（a）x1表示x1与自身褶积a次，（b）也是如此。
系统L是一个一般的线性系统，满足：L[ax1+bx2]=aL[x1]+bL[x2] L[x·？啄（t±？子）]=L[x]？啄（t±？子）（3）系统D-1为D的逆变换特征系统，是把经过处理的相加关系转化为褶积关系输出，即：D-1[ax1+bx2]=（a）D-1[x1]*（b）D[x2]（4）系统定义为这种D→L→D-1的标准格式，其突出的优点是，只要确定了特征系统D之后，它在整个反褶积过程中保持不变，并将一个非线性过程转化为线性运算。
地震零相位
地震零相位
地震零相位是地震子波的一种相位类型，通常指由2至3个或多个相位组成的地震脉冲，确切地说，地震子波就是地震能量由震源通过复杂的地下路径传播到接收器所记录下来的质点运动速度（陆上检波器）或压力（海上检波器）的远场时间域响应。
为了实现地震子波的零相位化，可以利用APC（Automatic Phase Correction）法，该方法运用的是最大方差模作为判断标准，通过修改地震道的相位角，当某个相位角对应的方差模最大时，那么这个相位对应的子波就接近零相位了。
零相位子波具有以下优势：第一，在相同带宽条件下，零相位子波的旁瓣比最小相位子波的小，也可以理解为能量集中在较窄的时间范围内，所以分辨率高；第二，在地震解释中零相位更有利于相位对比，更便于层位追踪解释；第三，零相位分辨率高，更有利于薄层的识别。
在地震资料处理中，零相位处理可以提高地震道的分辨率，更准确地对应地层界面，有利于薄层的识别和地震解释，是一种重要的地震道处理方法。
利用Ricker子波相移特性估计地震资料剩余相位
2021年4月名瀹錶4r f鲞找第56卷第2期•处理技术•文章编号：1000-7210(2021 )02-0249-08利用R i c k e r子波相移特性估计地震资料剩余相位刘传奇*周建科李宾王腾(中海石油（中国）有限公司天津分公司，天津300459)刘传奇，周建科，李宾，王腾.利用Ricker子波相移特性估计地震资料剩余相位.石油地球物理勘探，2021,56(2) ：249-256.摘要零相位地震资料具有分辨率高、易于解释等优点.然而最小相位假设下的反褶积处理很难保证地震资料为零相位，需通过后续技术手段估计剩余相位并做零相位校正。
通过研究Ricker子波的相移特性，发现主频不同的Ricker子波在相同相位旋转下，峰值时刻存在时差，该时差与相位旋转量具有线性关系。
在常相位假设下，对相位未知的地I数据做两次滤波处理，基于滤波后的数据求出同一反射界面的峰值时差，再通过峰值时差与相位的线性关系实现剩余相位的估计。
该方法不受判别准则的约朿，模拟合成记录和实际资料的应用结果证实了该方法的可行性和有效性。
关键词 Ricker•子波常相位相位旋转峰值时差剩余相位估计中图分类号：P631 文献标识码：A doi： 10. 13810/j. cnki. issn. 1000-7210. 2021. 02. 005〇引言众所周知，零相位地震子波具有旁瓣少、峰值时 刻与地层反射界面相对应等优点，能保证零相位地震资料不仅具有较高分辨率，还有利于地震资料的解释[19]。
确保地震资料的零相位性，无疑对构造精细解释、储层综合研究具有重要意义。
在现今实际地震数据处理中，一般基于地震子波为最小相位的假设，从地震道的相位谱中消去子波的相位谱。
然 而从现场采集到室内处理的各个环节，均可能对地震数据的子波相位特性产生影响.因此最小相位子波的假设条件通常是很难满足的，混合相位假设的子波更符合实情。
虽然混合相位地震子波假设能在某些特定条件下取得一定效果，但基于此假设的一些算法在计算量、实用性以及稳定性等方面大多还存在不足。
地震极性判断(地震剖面原理)
17/80
三、地震波频率的应用
未做滤波的正常剖面，同相轴形态清楚，
18/80
三、地震波频率的应用
高频滤波剖面，同相轴形态变胖，分辨率降低
19/80
三、地震波频率的应用
低频滤波剖面，同相轴形态变细，分辨率提高，同相轴形态不能分辨
20/80
三、地震波频率的应用
21/80
三、地震波频率的应用
频率为18HZ时，地下断裂的显示效果（最好）
正反射系数 界面为单轨
正反射系数 界面为双轨
15/80
二、地震剖面极性判断
对不同时段地震剖面提取子波，均为正极性地震子波。
500~1000ms 1000~3000ms
3）、提取地震子波判别法
3000~7000ms
16/80
汇 报 提 纲
一、地震剖面子波相位 二、地震剖面极性判断 三、地震波频率的应用 四、地震波形的应用
技术交流
地震剖面极性判断（地震剖面原理）
1/80
汇 报 提 纲
一、地震剖面子波相位 二、地震剖面极性判断 三、地震波频率的应用 四、地震波形的应用
2/80
一、地震剖面子波相位
零相位子波 最小相位子波
1）、地震子波类型
最大相位子波
混合相位子波
3/80
一、地震剖面子波相位
2）、实际地震子波
4/80
4、 对于非零相位地震剖面使用视极性的概念, 即视正极性、 视负极性； 5、 对于处理资料的极性，最终要标明
8/80
二、地震剖面极性判断
正极性
1、极性的规定
负极性
9/80
二、地震剖面极性判断
2、极性的判断方法
1）、声波人工合成地震记录法
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
·处理技术 ·
第 43 卷 增刊 2
地震数据处理中关于地震子波相位特性的探讨
苏贵仕 3 ① 沈克非 ② 丁学垠 ②
( ①东方地球物理公司研究院 ,河北涿州 072751 ; ②东方地球物理公司研究院海外业务部 ,河北涿州 072751)
苏贵仕 ,沈克非 ,丁学垠. 地震数据处理中关于地震子波相位特性的探讨. 石油地球物理勘探 ,2008 , 43( 增刊 2) : 121～124
关键词 地震子波 最小相位 预滤波 叠前反褶积 零相位化
1 引言
地震数据处理的目的是将野外采集的地震记录 用处理模块进行处理后得到成像好 、分辨率高的地 震剖面 ,地震记录可描述为地震子波与地层脉冲响 应或地下反射系数系列的褶积 。就某种意义上说 , 地震数据处理实际上是一个对地震子波不断改造的 过程 。
图 4 最小相位化因子
位的地震子波经过最小相位滤波之后才依然是最小 相位的子波 。有些公司采用零相位的滤波器进行预 滤波处理 ,严格来说是不合适的 。这是由于尽管滤 波的结果是振幅谱相乘 ,相位谱相加 ,零相位滤波器 没有改变地震数据的相位谱 ,但事实上最小相位子 波振幅谱的对数的希尔伯特变换就是最小相位子波 的相位谱 ,若最小相位子波的振幅谱被改变了 ,那它 便不再是最小相位了 。
1 24
石油地球物理勘探
位化之前的地震数据的叠加 ,为零相位的数据 。图 8 是最小相位化之后的地震数据的叠加 ,为最小相 位的数据 。图 9 是最小相位化后再进行叠前预测反 褶积之后的地震数据的叠加剖面 。从图中可以看 到 ,通过反褶积处理 ,剖面的分辨率提高了 ,这是反 褶 积的作用 。图10是未经零相位化处理 、直接经过
V=
∑y4j ∑y2j 2
→max
另一种是 Par simo ny 准则
∑ ln | y j | P -
∑| y j | Pl n | y j | P
j
∑| y j | P
→ mi n
j
式子中 y j 为滤波的输出 。经过扫描确定的相角或
相位函数就可以进行纯相位滤波 。利用纯相位滤波
进行零相位化 ,地震数据中强反射波将起主要作用 。
第 43 卷 增刊 2
苏贵仕等 :地震数据处理中关于地震子波相位特性的探讨
123
零相位的 高通 滤波 的结果 , 滤 波器 的频 带是 7 ～ 125 Hz ;第四道是经过两次零相位高通滤波的结果 ; 第五道是经过三次零相位高通滤波的结果 。通过比 较可以看出经过零相位高通滤波之后的子波强波峰 位置没有发生改变 。
2008 年
图 11 叠后地震数据子波经零相位化处理的结果
剩余静校正和精细速度拾取等处理后得到的最终叠 加剖面 。图 11 是图 10 经过零相位化之后的结果 。 零相位化因子是从叠加数据中直接求取的 。两图对 比可以看到在剖面右上方 400ms 左右尖灭更加清 晰 ,剖面的分辨率得到了提高 。
不同相位的反褶积 , 通常既是对地震子波的相 位特性而言 ,也是对反子波的相位特性而言 。从这 个意义上讲 ,反子波相位特性与地震子波相位特性 相同 。对于脉冲反褶积 :地下反射系数是白噪声和 地震子波是最小相位 。当输入的地震数据的子波是 最小相位时 ,其脉冲反褶积的结果也是最小相位的 。 对于预测反褶积 , 当预测步长是一个样点时它就是 脉冲反褶积 。当预测步长大于一个样点时 , 预测反 褶积的滤波因子不一定是最小相位 , 其结果是一次 反褶积的结果不一定是最小相位的 。对于两步法统 计子波反褶积 ,当其期望输出子波是零相位时 ,反褶 积的结果是零相位的 。当其期望输出子波是最小相 位时 ,其反褶积的结果也是最小相位的 。对于同态 反褶积和最小熵反褶积 , 它们对输入的地震数据的 地震子波相位特性没有特别要求 , 一般认为输入地
图 6 对最小相位子波应用零相位高通滤波的试验结果
2. 3 叠前反褶积 根据褶积模型 ,地震记录是地震子波和地下反
射系数的褶积 ,即 S j = w j 3 ej
其中 : w j 为地震子波 ; ej 为地下反射系数 。 反褶积处理就是要设计一个滤波因子 , 使地震
子波延续长度尽可能地减小 。如果这个滤波因子是 地震子波的反子波 , 反褶积之后的地震子波延续长 度将只有一个点即尖脉冲 , 反褶积结果即是地下反 射系数 ,这是理想化的结果 。一般情况下 ,反褶积之 后的地震子波有一定的延续长度 。
摘要 本文就实际地震数据处理中常用的地震子波最小相位化 、预滤波 、叠前反褶积和叠后地震子波零相位化 等四个处理模块对地震子波相位特性的影响进行了探讨 。当地震记录是由可控震源激发得到时 ,其地震子波 是可控震源扫描信号的自相关 ,呈零相位子波特性 ,此时需将零相位子波转化为最小相位子波 ,为后续的反褶 积处理做准备 ;当地震记录是由炸药震源激发得到时 ,其地震子波一般被认为是最小相位子波 ,在对此地震数 据进行叠前反褶积处理之前首先使用一个高通滤波器对其进行预滤波处理 (这个滤波器必须是最小相位的 ,因 为最小相位地震子波只有经过最小相位滤波器滤波后才是最小相位的) ,然后进行叠前反褶积处理 (主要包括 脉冲反褶积 、地表一致性反褶积等 ,该处理过程要求输入地震数据的地震子波是最小相位的) ,经过这些反褶积 处理后的地震数据的子波也认为是最小相位的 ;因零相位子波的地震剖面分辨率最高 ,所以最后还需再对叠后 数据的最小相位子波进行零相位化处理 。以上认识已在巴基斯坦 E 工区二维地震资料的处理中得到验证 。
4 结束语
地震数据处理的目的是提高地震数据的信噪比 和分辨率 。最终得到高信噪比 、高分辨率的地震剖 面 。而在频率域中则表现为使得振幅谱频带尽可能 宽且剩余子波是零相位的 。在处理时 ,有些模块对 地震数据的子波相位要求有假设前提条件 。一般情 况下 ,叠前反褶积输入的地震数据要求子波是最小 相位的 ,若是零相位的需要进行最小相位化 。根据 相同的振幅谱零相位的子波分辨率最高的原则 ,需 要对叠后数据的子波进行零相位化处理 。
最小相位化因子的求取需要用到可控震源激发
© 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved.
1 22
石油地球物理勘探
2008 年
的扫描信号 ,因为只有利用扫描信号进行自相关 ,才 能得到零相位子波 ,再将这个零相位子波转换为具有 相同振幅谱的最小相位子波 ,进而求取出匹配因子 , 这个匹配因子就是最小相位化的因子。图 1 为可控 震源的扫描信号 ,起始频率和终止频率分别为 7 Hz 和 72 Hz ,扫描信号长度为 12000ms。图 2 是扫描信号 的自相关 ,为零相位子波 。将这个零相位子波转换为 具有相同振幅谱的最小相位子波如 3 所示 。图 4 是 最小相位化因子 。将图 4 所示最小相位化因子应用 到地震数据就可以得到最小相位子波的地震数据 。
3 河北省涿州市东方地球物理公司研究院 ,072751 本文于 2008 年 3 月 20 日收到 。
的 。如果一个地震子波能量的大部分集中在末端 , 则我们说这个地震子波是最大相位的 。分辨率高的 地震子波在频率域上表现为振幅谱尽可能的拓宽 。 而具有同样的振幅谱的地震子波中 ,最小相位的地 震子波分辨率最高 。实际地震数据处理中对地震子 波相位特性有要求和影响的处理模块主要包括可控 震源地震子波最小相位化 、预滤波 、叠前反褶积和叠 后地震子波零相位化 。下面就这四个处理模块进行 具体的讨论 。
如果有测井合成记录 , 应以合成记录为标准进行相 位校正 。
3 实际应用
下面的地震数据选自巴基斯坦 E 工区 ,该区的 地震数据采集时利用可控震源激发 。图 7 是最小相
© 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved.
2 影响地震子波相位特性的四个处理 模块
2. 1 可控震源地震子波最小相位化 当地表条件较简单时 ,可采用可控震源激发 ,此
时因采集到的地震数据的地震子波是可控震源扫描 信号的自相关 ,所以是零相位的 。通常后续的叠前 反褶积处理要求输入的地震数据的子波是最小相位 的 ,为此需要将可控震源的地震子波进行最小相位 化处理 。
2. 2 预滤波处理 当叠前反褶积的地震数据存在比较严重的低频
干扰波时 ,通常使用一个高通滤波器对地震数据进 行滤波 ,这就是预滤波处理 。进行预滤波处理时 ,要 求这个滤波器必须是最小相位的 , 因为只有最小相
图 5 对最小相位子波应用最小相位高通滤波的试验结果
© 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved.
需要注意的是 ,地震数据经过最小相位滤波器 滤波之后 ,其强相位位置会发生变化 ,而地震数据经 过零相位滤波器滤波之后 ,其强相位位置不发生变 化 。图 5 是对最小相位的地震子波用最小相位滤波 器进行高通滤波的试验 : 图中第一道在 500ms 处 有一尖脉 冲 ; 第 二 道 是 最 小 相 位 子 波 , 其 频 带 是 7～72 Hz ,最大波峰在 500ms 处 ;第三道是经过一次 最小相位的高通滤波的结果 ,滤波器的频带是 7～125 Hz ;第四道是经过两次最小相位高通滤波的 结果 ;第五道是经过三次最小相位的高通滤波的结 果 ,子波强波峰明显发生了变化 。图 6 是对最小相 位的地震子波用零相位滤波器进行高通滤波的试 验 :图中第一道在 500ms 有一尖脉冲 ;第二道是最 小相位子波 ,其频带是7～72 Hz ;第三道是经过一次
本文得到了东方地球物理公司研究院海外业务 部的领导和专家的关心和帮助 ,在编写过程中得到 了周兴元教授的悉心指导 ,与李合群高级工程师也 进行了多次有益的讨论 ,在此表示衷心的感谢 。