子波与反子波的最佳提取

格式：pdf
大小：304.35 KB
文档页数：5

下载文档原格式

地震数据数字处理总结

中国石油大学（北京）《地震数据处理方法》勘查2011级复习重点总结第一章地震数据处理基础1、地震信号的特点：1）实信号2）离散3）有限长4）能量有限5）非周期2、采样定律内容：一个连续信号，如果其最高频率小于尼奎斯特折叠频率，即信号的采样频率大于信号最高频率的两倍，则利用离散采样后的信号可以恢复原始信号。

3、采样定律的应用条件：信号的采样频率大于信号最高频率的两倍，即：最高频率至少要在一个周期内采到两个样点4、采样频率、折叠（尼奎斯特）频率、信号最高频率定义：5、假频的定义：高于尼奎斯特频率的高频成分以尼奎斯特频率为中心向低频方向折叠，形成假的频率成分，称为假频。

6、假频的判断和计算：7、地震信号的频谱特点：1）有限带宽（带限）2）有一定主频（主频越高，分辨能力越强）8、判别相位性质的三种办法：1）相位延迟（不常用）2）能量延迟3）Z变换的多项式求根（根都在单位圆外，为最小相位（延迟）信号）9、一维数字滤波实现方法、具体步骤：1）频率域：实现方法：（以零相位为例，翻译略）具体步骤：a、地震频谱分析：确定分析有效频率范围b、设计滤波器：压制噪声保留有效信号c、地震记录FFT变换：标准化变换长度d、进行滤波运算：振幅谱相乘相位谱相加e、滤波结果IFFT2）时间域：（也叫褶积滤波）实现方法：（以零相位为例，翻译略）具体步骤：a、地震记录频谱分析：确定中心频率、带宽b、设计滤波器：确定滤波算子长度（频带越宽，长度越短）c、确定滤波因子离散值：双边对乘实参数d、进行滤波运算：地震记录与滤波因子褶积10、伪门的定义：对连续的滤波因子用时间采样间隔离散采样后，得到离散的滤波因子，若再按离散的滤波因子计算出与它相应的滤波器的频率特性，这时在频率特性的图形上，除了有同原来连续的滤波因子的频率特性对应的“门”外，还会周期性地重复出现很多“门”，这些门称为“伪门”。

产生“伪门”的原因：由于对滤波因子离散采样。

11、吉布斯现象：当对滤波因子用有限项代替无限项时，在原始信号突变点（间断点）处，通过信号出现的明显的振荡现象。

常用地震子波提取方法简介

科学论坛
Ｉ叠
ＣｈｉｎａｓｃｉｅｎｃｅａｎｄＴｅｃｈｎｏｌｏｇｙＲｅｖｉｅｗ
基于学生个性化能力培养行动导向在数控技术实训
教学中的应用研究
高生祥
（浙江工业职业技术学院３１２０００）
２．行动导向简介
育对象的气质、情绪、认知、兴趣、能力、性格、价值观和信念等进行人格整合和
个性优化。通过应用行动导向法教学能够发挥学生Байду номын сангаас主体的主动性，调动学生的求知欲，促进学生个性化发展。将六步法行动导向应用在数控技术实训教学中，
具体教学实施过程如下：
１）、资讯
主要是下达典型工作任务，明确任务的教学目的、重点和难点，通过任务的下达，让学生去自主查阅完成任务的相关资料，培养学生的个性化查阅资料的
能力。
２）、计划
在２０世纪８哞代德国的双元制职业教育，就是一种以“ 实践为导向” 的教育，它注重实践性教学环节，突出职业实践能力的综合培养，其实质是：教学组织以学生为中心，教师是学习过程的组织者与协调人，遵循“ 资讯、计划、决策、实施、检查、评估这一完整的“ 行动 ” 过程序列。行动导向的学习主要探讨认知结构与个体活动之间的关系；与之不同的是，行动导向理论以人为本，强调人在实现目标过程中进行反思的重要性，注重

地震数据处理第三章：反褶积

b(t ) o(t ) * g (t ) * (t ) * d (t ) * i(t ) o(t ) * fg (t ) * f d (t ) (3 - 4)
式中 o(t ) — 震源子波； g (t ) — 地层响应； (t ) — 透射响应； d (t ) — 地面接收响应；
i (t ) — 仪器响应；
（3-34）
将上式两端乘以
A( z ) R( z ) Z
M
zM
，则有：

M
M ( ) Z
M
(M ) Z 2 M (M 1) Z 2 M 1 (0) Z M (1) Z M 1 ( M ) Z 0
（3-18）
B(e ) | X (e ) | e
j
j
j ( e j )
（3-19）
( e j ) 未知，现在来确定它
•假如地震子波是最小相位的物理可实现序列，则其z变换为：
B( z) b0 b1z 1 b2 z 2
B( z ) 0 , 对下式由物理可实现性知：当| z | 1 时，
根据“最小相位序列z域零点在单位圆内”这一特点，选出模小于1的根，便可组成最小相位子波，其z变换为：
B1 ( z ) b0 (1 z1 z 1 )(1 z 2 z 1 ) (1 z M z 1） b0 b1 z -1 bM z -M
由于 ( ) ( )
A( z ) 应有2M个根。鉴于系数均为实数，所以显然， 2M个根是M对互为倒数的，即若
z01 e j , (| | 1)
则另一根为：
1 1 j z02 e z01
根据这M对根在单位圆内、外的位臵，可以组成2M个不同相位的地震子波，其中必有一个是最小相位，一是最大相位的。

地震反演方法综述

地震反演技术简介在上世纪70~80年代，地震反演作为地球物理学的一个重要进展得到了广泛的赞扬，获得广泛应用；地震反演技术能够帮助解释人员确定地层单元而不仅仅是通过反射波确定地层单元的边界，而且能直接进行深度域成图。

在一个竞争的市场环境中，开发出了很多不同的反演算法，在基本递归反演方法的基础上不断取得进进展，一下简要介绍几种基本的地震反演方法。

主要分三大类：1、基于地震数据的声波阻抗反演：其结果有两种：相对阻抗反演（常说的道积分）与绝对阻抗反演。

主要算法有：递归反演（早期的地震反演算法）与约束稀疏脉冲反演（优化的地震反演算法）。

这种反演受初始模型的影响小，忠实于地震数据，反映储层的横向变化可靠；但分辨率有限，无法识别10米以下的薄砂层。

2、基于模型的测井属性反演：此种反演可以得到多种测井属性的反演结果，分辨率较高（可识别2-6米的薄层砂岩）；但受初始模型的影响严重，存在多解性，只有井数多（工区内至少有10口以上的井，分布合理，且要求反演的属性与阻抗相关），才能得到较好的结果。

3、基于地质统计的随机模拟与随机反演：此种算法可以进行各种测井属性的模拟与岩性模拟，分辨率高（可识别2-6米的薄层砂岩），能较好的反映储层的非均质性，受初始模型的影响小，在井点处忠实于井数据，在井间忠实于地震数据的横向变化，最终得到多个等概率的随机模拟结果；但要求工区内至少有6-7口井，且分布较合理，才能得到好的模拟结果。

道积分道积分技术出现，为广大少井无井地区岩性及油气预测提供了新的途径，它能得到类似于虚速度测井的新方法，其结果对应于地层的波阻抗，它最大优点是不像虚速度测井那样依赖于井的资料和地球物理学家的经验。

尽管道积分剖面不能像GLOG波阻抗剖面那样反映地层绝对速度，而只能反映其相对速度大小，但是它反映出的层位与GLOG剖面是一样的，甚至在反映的细节上还比它多，对薄层识别也非常有利，因此道积分剖面能用于岩性和油气层解释。

谱比法求取地震子波

振幅
0 5 0 0.02 0.04 0.06 0.08 0.1 0.12 时间 (s) 0.14 0.16 0.18 0.2
0.2
0
0
50
100
150
200
250
雷克子波，反射系数序列及合成的地震记录
原始子波（红色）、合成记录（蓝色）、提取的子波（绿色）振幅谱
绿色的是雷克子波的振幅谱，蓝色的是合成记录的振幅谱，红色的是用二次谱模拟方法提取的子波的振幅谱图，提取的子波振幅谱和原始的雷克子波振幅谱几乎相同
0
100
200
300
400
500
600
700
选取250m处与0m处的子波振幅谱通过谱比法拟合得到的直线
7
2.3 谱比法存在的问题不同时刻的子波的振幅谱难以确定
两种解决办法： 1.地震记录的振幅谱直接作为子波振幅谱进行计算（王秋成 2011） 2. 从地震记录的振幅谱中近似模拟出子波振幅谱（Rosa 1985,1991)
谱模拟法模拟子波（Rosa 1985,1991）
8
2.4 研究思路
基于谱比法，利用二次谱模拟法来提取不同时刻的地震子波振幅谱, 进而可以估算Q值
1.模拟不同类型子波在地层中传播情况 2.进行二次谱模拟方法的分析与应用
3.Q值的实际应用：反Q滤波
9
2.5 衰减子波波传播模拟
B( f , t ) A(t ) B0 ( f ) exp(
15
2.8 二次谱模拟方法的分析与应用（1）实际地震资料的不同时刻子波提取
-2 0
-1.5
பைடு நூலகம்
-0.5
0.5
1.5
-1

探地雷达信号子波提取及互相关滤波

t =1
式中 rxy (τ) 是雷达子波和雷达信号之间的互相关
系数。 T 是互相关滤波的时窗。
3 数值模拟
311 水平层状介质模型剖面图 1 是水平层状介质合成雷达记录。在该剖
面中 ,反射波的起跳点在深度 150m 处。图 2 是加入随机噪声后的雷达记录。在该雷达记录中 ,反射波的起跳点完全淹没在随机噪声中。图 3 是从图 2 的剖面的记录道中提取的零相位子波。图 4 是用提取出的子波做适当截取后和对应的加入随机噪声后的记录道进行互相关滤波处理后的雷达记录。由图 2 和图 4 可见 ,经过互相关滤波处理 , 记录道的波形明显变得较平缓 ,雷达信号中的噪声被有效压制 ,而且在深度 150m 处有一个明显的极大值表示反射波的到达。
3 Rogerb L Roberts and Jeffrey J Daniels. Modeling near - field GPR in three dimensions using the FDTD method. Geophys2 ics ,1997 ,6 (24) :1114～1126
4 O Lazaro - Mancilla , E Gomez - Trevino. Ground penetrating radar inversion in 1 - d : an approach for the estimation of electrical conductivity , dielectric permittivity and magnetic permeability. Journal of applied Geophysics , 2000 , ( 43) : 199～213
由于反射系数ξ( t) 为白噪声 ,所以

基于动态子波估计的吸收系数提取方法

摘要：在基于褶积模型的子波提取中，多假设地震子波为最小相位性质。虽然在地震资料大
处理中进行了相位转换和Ｑ补偿处理，从实际地震记录中提取的地震子波，多也表现为混但大合相位性质。提取吸收系数大多采用地震记录谱模拟方法，即用地震记录振幅谱平滑代替子波
基金项目：国家自然科学基金委员会与中国石油化工股份有限公司联合基金资助项目（０３９３４７９０）
收稿日期：２１００—１１—１６改回日期：２１ —０２０１３— ３
３期
曹健等：于动态子波估计的吸收系数提取方法基
２５６
关键词：动态子波；子波提取；吸收系数；波估计子中图分类号：Ｐ６１４４３．３文献标识码：Ａ
剖面表现为分辨率低，于进行储层的精细解难
０前言
吸收系数是油气预测与各向异性研究的一个
释。
辨率，但这种转换只是一种近似。反Ｑ滤波方法具有补偿大地的滤波作用，能够调整地震记录中子
不是白噪音，震记录为子波与反射系数的褶积。地
在剖面中任意提取两个反射层的子波，根据前面的
波的振幅谱和相位谱。在理论上，种反Ｑ滤波这可将大地非弹性吸收引起的相位移校正为震源子
气。
谱对数差的斜率计算吸收系数，其剖面分辨率高，容易对储层进行高精度解释研究。

4.反褶积

iφ ( eiω )
假设地震子波是最小相位的，则地震子波满足因果关系，具体讨论见教材。
3.2反滤波
再假设地震子波是零相位的，地震信号满足
自相关法
(ω ) |= 1 |R
(ω ) | | W (ω ) || R (ω ) | | W (ω ) | |X = =
(ω ) |2 =| X (ω ) |2 |W
得到反子波时间序列并与地震记录进行褶积滤波
w '(t ) = {w '0 , w '1 , w '2 , , w 'm }
= r (t )
w '(τ )x(t − τ ) ∑ τ
最小相位-稳定其他相位-不稳定
3.3最佳维纳滤波及最小平方反褶积
最小二乘拟合/优化思想已知样点（x1,y1）,(x2,y2),…,(xn,yn)
= r (t ) w = 't )* x(t ) w '(t )* w(t )* r (t )
子波与反子波满足
w '(t ) * w(t ) = δ (t )
已知地震子波求出反子波，将反子波与地震记录褶积，即可求出反射系数，这个过程叫作反褶积。
3.2反滤波
地震子波的求取
在进行反褶积处理时，通常必须知道地震子波的形状。地震子波求取得是否准确对反褶积结果的影响很大。求取地震子波的方法较多,常用包括：（1）直接观测法（2）自相关法（3）多项式求根法（4）测井资料求子波（5）对数分解法
基本原理
最佳维纳滤波是数字滤波中的一大类滤波方法。它是在滤波器实际输出与期望输出的误差平方和为最小的情况下，确定滤波器的滤波因子的，因而称为最小平方滤波。已知输入信号 b(t ) = {b(0), b(1), b(2), , b( n)} 现在要求设计一个滤波器，其滤波因子为使得滤波后的实际输出为

地震反演原理及其应用

地震波阻抗反演是地震数字处理的终结结果之一，是油气藏勘探、开发中储层横向预测的重要手段。而且随着勘探、开发程度的提高，对反演资料的精度要求也越来越高。但是，由于现有的地震资料受地震波激发和采集系统本身条件的限制，就决定了所记录的地震波仅有相似的振幅信息，分辨率有限且包含各种干扰。因此，波阻抗反演无疑是一个多解的非线性问题。
目前，以物探局范祯祥先生开发的《非线性波动方程地震波反演技术软件—ANGEL2001》为代表的地震反演软件包已投入生产应用。本技术的基本技术思路是：采用有限元波动方程对地震波剖面数值模拟，借助于最优选择理论对地震波反映的物性参数与几何形态进行逐次逼近，以实现物性参数的反演。在此基础上借助于神经网络分析对所反演的物性参数进行非线性标定。小波技术，模糊识别技术，分形技术穿插应用其中。
即而反演过程则是估算个子波的地震反演概述正演合成记录tntwtrtsiiiiiiiir?1111?而反演过程则是估算一个子波的逆反子波用反子波和地震道进行褶积运算得到反射系数rt然后由上式导出的递推公式逐层递推计算出每一层的波阻抗即用gardner公式从波阻抗中分离出速度和密度
地震反演原理及其应用
引言：
s(i) r ( j ) w(i j 1) n(i)
jHale Waihona Puke （ 1）（ 2）这里＊意思是褶积
关于褶积模型的假设条件: •叠后数据 – 地震道是零炮间距的； •没有多次波；
•没有 AVO 效应；
•噪声是随机的, 即是白噪的, 与地震不相关；没有相干噪声。 •子波是恒定的 – 不随时间变化。
这里
( j ) = 第j层的密度 ( j ) = 第j层的纵波速度
•
反演的目的是为了从地震道本身估算地层速度。显然，首先需要从褶积模型中提取反射系数的一个估算值，因此，我们可进行反褶积的相关处理。在讨论反褶积和反演之前，让我们更详细地讨论一下褶积模型的两个主要组成部分：反射系数和地震子波。

地震子波提取方法综述

地震子波提取方法综述
杨培杰;印兴耀
【期刊名称】《石油地球物理勘探》
【年(卷),期】2008(043)001
【摘要】地震子波提取方法总体上可分为确定性子波提取方法,统计性子波提取方法两大类.确定性子波提取方法不需要对反射系数序列的分布作任何假设,能得到较为准确的子波,但很容易受各种测井误差的影响,尤其是声波测井资料不准而引起的速度误差会导致子波振幅畸变和相位谱扭曲;统计性子波提取不需要测井信息,也可以得到子波的估计,但需要对地震资料和地下反射系数序列的分布进行某种假设,所得到子渡精度与假设条件的满足程度有关.文中对确定性子波提取方法和统计性子波提取方法进行了详尽的评述,并认为可将通信工程中的子空间分解法的思想用于混合相位子波提取.子空间法混合相位子波提取基于信号和噪声子空间正交的原理,建立二次型的目标函数,不需要对地层反射系数做任何假设,因此子波提取准确度高且稳定性好.
【总页数】6页(P123-128)
【作者】杨培杰;印兴耀
【作者单位】山东省东营市中国石油大学地球资源与信息学院,257061;中国石油大学(华东)地球资源与信息学院
【正文语种】中文
【中图分类】P61
【相关文献】
1.基于高阶统计量的地震子波提取方法评价 [J], 陈健;戴永寿;张亚南;魏玉琴;丁进杰
2.高效ARMA模型高分辨率地震子波提取方法 [J], 张亚南;戴永寿;王少水;彭星;牛慧
3.地震子波及其提取方法研究 [J], 陈建;王科文
4.基于矢量预测的SIMO系统混合相位地震子波提取方法研究 [J], 王德营;李振春;张丰麒;董烈乾;丁成震;黄建平
5.深度域地震子波提取方法综述 [J], 张杰;陈学华;蒋伟
因版权原因，仅展示原文概要，查看原文内容请购买。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

年华木石油学院学报第期

子波与反子波

的

最佳提取

陈伟钟

物探教研室

提要

本丈分折了有升声阻杭技术中对井时差与有限数据最小平方法中,期望猪出

位笠对子波与反于波提取的影响在此墓袖上,提出了利用实际猪出与期望输本之

间的相似来数作为别别准瓤村记录段时窗进行滑动扣描,来选择最佳时窗位主,从而解决了这个呵题理论和实际资朴的试井结果表明该方法是切实可行的而且

是很有必要的

一目前的方法及存在的问题

原理

子波与反子波的提取是当前声阻抗技术的关键有井声阻抗技术中,提取子波与反子

波使用的资料是时域测井反射系数与井旁记录考虑到一般相位的子波其反子波是无限长的双边序列所以用下列双边反褶积方程来提取子波,幻

刀一、一厂一了…一卜一一

…

……

二一

式中

—井旁记录。的自相关

。—井旁记录,与测井反射系数序列氛的

互相关

一反子波

实际计算岭和。时所用的记录段。与反射系数段夸扭是根据用户提

供

的井段反射系数序列的起始时尚来确定反射系数序列的时间坐标然后在同一

时窗内截取

的么存在向题

测井资料往往不是从地面起始的井段反射系数序列的起始时间是用户按统计资料或地震测井资料提供的这往往会存在甲定的误差从而导致同‘时窗内截取的反射系数段言

多界与记录段。之间存在系统时差这种系统时差通常称为对井时差由于它的存在即使是

理想的褶积模型时窗一端附近的。的也不能为同一时窗内某个部分的反射系数与子波的褶积,而另一端附近的反射系数也与时窗内记录无关利用这样的资料必然会使提取的反子波存在误差有限区间最小平方法,虽然是解决时窗边界影响的但是当这种系统时差超过了一个子波长度时此法也就失去了意义此外,一般相位的子波其反子波是无限长的双边序列并且能量的分布也不一定是对称的而求解方程的阶数总是有限的所以把反子波安排在解时窗中的不同位置等价于期望输出的位置不同将会使提取的反子波有不同的精度

二子波与反子波的最佳提取方法

基本原理

由以上分析可知对井时差对反子波提取的影响实质上是由记录与反射系数序列之间的相对误差引起的因此为讨论方便起见设测井反射系数序列的时间坐标是准确的而

认为地震记录的时标存在系统时差

在反射系数序列时标确定的情况下子波零点的变化等价于记录零点的变化或记录时

标的变化,相反也成立因此子波零点定义的不同将会影响到对井时差的定义确切地

说对井时差是在子波定义为物理可实现的前提下地震记录与测井反射系数序列在时标上存在的系统时差同理子波零点的变化也会引起其反子波零点的变化但变化方向相反,反之亦

然

因此调整记录时标等价于调整子波零点同时也等价于调整反子波的零点双边反褶积方程决定了反子波零点在解时窗的中点所以调整记录时标也等价于调整反子波在解时窗中的位置由此可见解快前面提出的两个问题都归结为调整记录时标为此,将仅有的不太长的反射系数段固定让记录段时窗滑动扫描从中选取一个合适的记录时窗位置为确定这

个合适的时窗位置选择实际输出与期望输出之间的相似系数作为判别准则因为在声阻抗技术中只利用实际输出与期望输出之间的相对幅度关系且相似系数是归一化的这样

就能获得这一准则下的最佳反子波理论记录试验试验采用了某井的测井反射系数并用其井旁记录上提取的子波制作的合成记录试验分

别采用了无限区间最小平方法与有限区间最小平方法并对两种方法进行了比校

无限区间最小平方法此法即为方程式此方程是以所用的记录。满足褚的完整褶积为

基础的且误差能量是在一十上计算

的

试验采用方程的阶数为时窗长度为采样率为试验结果如图所示图中给出了相似系数与对井时差的关系曲线可以看到最大相似系数并不出现在对

井

时差为零的地方这正是前面提到的反子波的能量不一定关于原点对称分布的原因将最佳位置处求得的子波与合成记录所用的理论子波比较示于图由图可见所求的子

波

是

相当好的此时获到的近似反射系数与实际反射系数相比较两者在波形上不存在什么差

别

多革相似系效对井叶差一》

浮才子碑子刃口劝可卜日日

图最佳反子波的选择

相似系数一对井时差曲线理论子波最佳

反

子波位置上提取的子波最佳位笠上未得的反子波

有限区间最小平方法及其改进效果

上面这种无限区间最小平方法是建立在完整褶积基础上的这在实际应用上是不可能满足的。因为实际工作中只能从几秒长的记录上截取一段用于提取反子波所以即使

在

最

佳位置,仍然存在时窗边界的影响尤其在时窗较短时窗端点附近有较强反射层时,影响就更大文献〔〕提出的有限区间最小平方法就是专门解决这一问题

的

有限区间最小年方法求解反子波的方程为

·“,,耳。“‘,一‘,“‘,

一‘,·暴卢‘,’“一‘’‘

一”

”

式中。反射系数时窗起点反射系数时窗的终点

︸一︸卜

丁—反子

波

。—记录

占—反射系数序列

将式写成矩阵的形式即为一个实对称矩阵方程理论记录试验所用的资料与前面的相同采用方程的阶数为,时窗长度为试

验分别用有限区间最小平方法及无限区间最小平方法计算了相似系数曲线试验结果如图所示可见,与无限区间最小平方法相比有限区间最小平方法对短时窗有明显的效果同时可以看到对井时差较大时有限区间最小平方法所采取的改进措施也将失去意义应用实例

对某测线上井和井提供的测井资料及井旁记录进行了实际计算井的测井资料在

时域长度为井为,,井旁记录与侧井反射系数序列都是以采样由于井反射系数段较长采用了无限区间最小平方法而井较短采用了有限区间最小平方

法用户提供井的井段初始时间为。井为计算结果如图

所示曲线的横坐标为记录段滑动时窗的起点纵坐标为相似系数

。

名口区问最小平方法

勺毕

户理论子波

上对少抑寸差

一一一或丽

一

了丽一一一一

间

有限区间尿小乎方法且州山

反子波

相似系数曲线比较

图有限区间最小平方法与无限区间最小平方法的效果比较

对井在用户提供的初值附近计算出来的相似系数很小甚至有的出现了负值进一步扫描才得到了图的曲线通过这一计算就有了一个量的标志表明对这口井若按用户提供的时间起点转换声阻抗剖面将一定是失败的从曲线可以看出最佳时窗起点

偏离用户提供的时间达。最佳位置处求得的近似反射系数与实际反射系数比较于图中,可见两者有较好的相似性因此换算出来的阻抗在井段范围内与测井得到的阻抗

也能较好地相似井的计算结果表明最佳时窗起点也偏离用户提供的起始点达。在最佳位置处求得的近似反射系数与实际反射系数比较于图中两者也有较好的相似性

相似系

睡

二,几卫卫

吕冈犯

图“最佳位筑求出的近似反封系吸,己录时窗起始时间而两矿气丽厂面石一汤顽厂一丽苏一丽井最佳位五的选择及效果

最佳位置选出以后,分别用机上模块求出子波如图所示其频谱如图,所示可见在最佳位置处求出的子波不仅在波形视周期

而且在频谱上都是较为合理的儿卫、土,

,,,

图井最佳位呈的选择

及效果

扩子

波的振幅诺

二空井最佳位置求得的子波

卜井最佳位置求得的子波

子波卜的振幅讲

图最佳位主选择对实际资料上子波提取的效果

三结论

实际资料中存在的对井时差对子波与反子波提取的影响是不可忽视的否则提高提取反子波精度的两种措施—有限区间最小平方法及双边反褶积就失去了意义因此有必要采取文中提出的扫描选取法

在此项工作中得到了物探局研究院刘家瑜老师的指导及周兴元工程师的帮助,在此表

示感谢

年月收到

参考文献

〔〕刘雯林等,合成声波测井,石油地球物理勘探,年月

〔傅才芳层序反褶积与合成声波测井石油地球物理勘探年月

〔〕查中沂、提取子波的一种方法,石油地球物理勘探年月