数字处理ch32子波提取与子波整形反褶积
地震子波处理的二步法反褶积方法研究
地震子波处理的二步法反褶积方法研究
地震子波处理的二步法反褶积方法研究
针对玛湖斜坡区三块三维地震资料和赛汉塔拉凹陷二块三维地震资料连片处理中的特点,结合地质任务和处理目标要求,提出了地震数据连片处理中的地震子波处理的方法.该方法主要体现了两次反褶积,一次是采用地表一致性反褶积,将不同震源的频带拓宽到一个标准上;再一次采用相位校正反褶积,将不同震源的数据校正到相同相位上.为了保证提取的相位校正反褶积算子稳定,采用叠后地震道提取(主要考虑到叠后地震道信噪比高,算子稳定性强),然后将该算子应用到叠前地震道,进行相位校正.
作者：王西文胡自多田彦灿王红旗王述江肖明图马建华WANG Xi-wen HU Zi-duo TIAN Yan-can WANG Hong-qi WANG Shu-jiang Xiao Ming-tu Ma Jian-hua 作者单位：中国石油勘探开发研究院西北分院,兰州,730020 刊名：地球物理学进展ISTIC PKU英文刊名：PROGRESS IN GEOPHYSICS 年，卷(期)：2006 21(4) 分类号：P631 关键词：连片处理子波处理地表一致性反褶积相位校正。
地震数据处理第三章反褶积PPT
B(ej)|X(ej)|
零相位地震子波：
b(n)21 |X(ej)|ejnd （3-32）
3、多项式法
选择一段质量较高的地震记录
x(n),(n0,1 ,2,.M ..),
褶积模型： x(n)b(n)*r(n)
设反射系数为白噪声序列，则记录自相关与子波 自相关等价，即
xx()N0bb (), N0为系数
其中：
s(t, x) 偏移距为x的地震记录； (t) 垂直入射时的反射系数序列； T1(t, ) 多次波效应(为其周期)； M (t, x) 炮检距对时差的影响； SL(t) 排列损失或球面扩散效应； T2(t, ) 时变吸收或非弹性衰减效应； RI(t, x) 与炮检距有关的浅层混响效应和记录系统的影响； b(t) 震源子波； N(t) 全部附加噪声（包括规则的和随机的）。
压缩地震子波，提高分辨率。 可以压 制多次波和短周期鸣震等干扰，提高地震 资料信噪比。
震源爆炸使地下介质形成三个区域：
震源爆炸产生尖脉 冲传播到弹性区起始 边界时，已经变成了 有一定延续时间的稳 定波形——地震子波。
地层对震源脉冲的改造作用，相当于对 它进行了一次低通滤波，此滤波器常称为 大地滤波器。
地震记录自相关rxx( )的z变换为
Rxx(z) X(z)X(z1)
B(z)B(z1)(z)(z1)
B(z)B(z1)
将 zej 代入，有：
X (e j )X (e j ) B (e j )B (e j )
由于 b：
X(ej) X*(ej) B(ej) B*(ej)
2、自相关法
选择一段质量较高的地震记录，时窗长度为T：
x(n)(,n0,1,2,..M .),M , T1 t
利用VSP子波替换提高地面地震资料的分辨率
利用VSP子波替换提高地面地震资料的分辨率于茜【摘要】子波反褶积是地面地震处理中提高成像分辨率的常用方法,而做好子波反褶积的关键是计算反褶积算子.地面地震资料受本身分辨率和信噪比的影响,处理中难以得到精确的反褶积算子.本文利用VSP资料的高分辨率、高信噪比的优势,从中提取子波,然后用地面子波向VSP子波进行匹配,计算出较为精确的反褶积算子,用于叠后地面地震资料子波替换反褶积处理.实例处理结果表明,此方法能使地面地震剖面的分辨率明显提高,同相轴的连续性得到改善,波阻特征有所加强,整体质量得到较大改进.【期刊名称】《物探与化探》【年(卷),期】2010(034)004【总页数】3页(P546-548)【关键词】地震勘探;VSP;子波;反褶积;分辨率【作者】于茜【作者单位】中海石油(中国)有限公司天津分公司,天津,300452【正文语种】中文【中图分类】P631.4子波反褶积是地震资料在提高分辨率处理过程中常用的一种方法。
该方法以某个子波作为地震波传播的已知子波，以子波的逆(即反子波)作为实际地震记录的反褶积算子［1］。
但由于地面地震接收干扰严重，高频成分在中深层衰减和吸收作用较大，地面地震得到的叠后剖面中深层频率偏低，处理中难以提取真实的地震子波及反褶积算子。
所以纯粹利用地面地震资料的信息来改善地震资料质量的能力是有限的。
VSP在地面潜水面以下激发，在井中接收，避开了低速带的影响，射线传播距离短，信号强，噪声相对地面弱，因而VSP资料具有较高的信噪比和分辨率，特别零偏VSP资料往往具有高频子波信息，所以可以利用VSP资料提高地面地震资料的分辨率［2］，国内外已有这方面的介绍［3－6］。
笔者提出了一种子波替换法，利用VSP资料的特有优势，从VSP资料中提取较为准确的地震子波，然后对过井的地面地震资料进行子波替换反褶积处理，实际数据处理结果效果明显。
从零偏VSP记录上截取时窗长度为m的数据: w(0)，w(1)，…，w(m－1)，w(m)。
地震资料数字处理技术
X (ω ) = X (ω ) H (ω )
∧
………….. (2-1-9)
四学时
傅氏变换 连续信号傅氏变换
X (ω ) = ∫ x(t ) e
−∞ ∧ +∞ − jωt
dt
1 x(t ) = 2π
∫
+∞ ∧
−∞
X (ω ) e dω
jωt
离散信号傅氏变换：
X Δ ( f ) = Δ ∑ x( n)e
x(n) ≤ M 。
n = −∞
∑ h(n) < ∞
+∞
………………
(2-1-14)
滤波器能量有限输出
这 是 指 如 果 输 入 x(n) 的 能 量
§1.2 地震处理流程
几何扩散校正：通过给数据加一增益恢复函数以 校正波前（球面）扩散对振幅的影响。 建立野外观测系统 ：把所有道的炮点和接收点 位置坐标等测量信息都储存于道头中以保证各道 的正确叠加 。 野外静校正 ：对陆上资料，把所有炮点和接收 点位置均校正到一个公共基准面上以消除高程、 低降速带和井深对旅行时的影响。
剩余静校正
地层断开引起的原因：由于障碍物 的影响没有布校波器
叠后处理
绕射波，滤波后可 直接作解释
偏移处理
顶
前积反射 断点比较清楚
第二章 数字滤波
本章主要回顾和介绍数字滤波器的有关 知识，以及利用干扰波与有效波在频率、 传播方向、速度以及能量等方面的差异进 行干扰波压制或消除，从而突出有效波， 提高地震资料的质量和精度的方法原理。 §２.1 概述（４） §2.2 一维滤波 （６） §2.3 二维滤波 （４）
– 特殊处理（目标处理）：针对不同目的采用
的特殊处理手段。
子波提取——精选推荐
子波提取褶积模型是所有反演的基础:地震道=地震子波*反射系数+噪声频率域内, 褶积则为乘积的关系.反演相当于地震道除以地震子波, 得到反射系数:反射系数=地震道/地震子波频域内窄频段的子波限制了信息的获取范围.The narrow band wavelet restricts the available range of information in the frequency domain.地震子波完全由它的振幅谱和相位谱来定义:The Wavelet is defined completely by its amplitudespectrum and its phase spectrum:在有限频率范围内, 相位谱通常可近似为一条直线. 直线的截距是子波的常数相位旋转, 它是子波的最佳表征. 直线的斜率标示着子波的时移.The intercept of the line is the constant phase rotation which best characterizes this wavelet.The slope of the line measures the time-shift of the wavelet.极性的约定:极性约定是一个特殊的子波相位问题. 默认的约定便是: 声阻抗的增加在零相位的地震数据上代表一个波峰.A special wavelet phase issue is the Polarity convention.The default convention is that an increase in acoustic impedance is represented as a peak on zero-phase seismic data:另一个默认的约定便是: 声阻抗的增加在零相位的地震数据上代表一个波谷.The alternate convention is that an increase in acoustic impedance is represented as a trough on zero-phase seismic data:使用ì极性约定菜单î可以设置极性约定:The polarity convention is set using the SyntheticPolarity Convention menu:地震子波在时间和空间上都存在着变化, 即具有时变性和空变性, 这是基于以下几个原因:Wavelets in the earth vary both laterally (spatially) and temporally for a variety of reasons:近地表效应(空变)Near surface effects (space variant)频率吸收(时变和空变)Frequency-dependent absorption (space and time variant)层间多次波(时变和空变)Inter-bed multiples (space and time variant)NMO 拉伸处理过程中的人为因素Processing artifactsSTRATA 假定子波是常数, 不随时间和空间变化: 时间不变性: 这意味着反演就是在有限的时窗内求最优化的波阻抗Time invariant: This means that the inversion is optimized for a limited time window.空间不变性: 这意味着去除子波的空变后被最优化处理. 通常, 许多方法有可以用来提取子波. STRATA中用了以下几种:In general, a variety of methods can be used for wavelet extraction. Some are available in STRATA. (1) 仅用地震数据估计地震子波的振幅谱. 假设相位谱已经从别的渠道得知.子相关autocorrelation最大熵谱分析maximum entropy spectral analysis交互谱分析cross spectral analysisSTRATA中统计子波的提取用自相关: Statistical wavelet extraction uses the autocorrelation(2) 单独使用地震数据估算振幅谱和相位谱Estimate both amplitude and phase spectra from the seismic data alone.最小熵子波估计高阶力矩法higher order momentsSTRATA 不用这种方法, 因为STRATA认为该方法不可靠.(3) 使用给定的测量数据估计振幅谱和相位谱Estimate both amplitude and phase spectra using deterministic measurements.海洋信号marine signaturesVSP 分析STRATA中, 以ASCII文件形式读入外部子波(4) 用地震和测井资料估算振幅谱和相位谱Estimate both amplitude and phase spectra using both seismic and well log measurements.STRATA中用测井资料提取全子波.(5) 用地震资料和测井资料估算振幅谱和常数相位谱STRATA中用测井资料提取常数相位子波.STRATA中提取子波的方法:第一步, 是否用测井资料来估算子波的相位. 关键是看测井资料与地震资料的相关性是否好. 通常情况下, 必须首先进行手动校正测井曲线. The critical issue for this decision is how well the logs tie the seismic data. Usually, manual correlation must be done to align the logs first.1 提取统计子波(不用井资料):这个过程只是通对地震道进行自相关计算子波的振幅谱, 并假设已知子波的相位.主要参数:ï道范围(通常设置为较大值以增加统计所用的道数) Trace range (usually set this large to increase statistics)ï时窗(至少应该为子波长度的两倍)ï子波长度(取决于层厚和分辨率, 层厚一般取200ms, 薄层取50~100ms).2 用测井资料提取子波:用测井资料提取子波:此方法用测井资料估算子波的振幅谱和相位谱. 效果取决于测井曲线和地震道的相关程度.主要参数:选择要用的井(只用标定效果好的井)道范围(距井的距离)时窗子波长度3 用测井资料计算单一常数相位值该方法使用地震道的自相关计算子波的振幅谱, 与统计子波提取方法中一样, 用测井资料计算子波的相位谱, 并且相位谱被近似为一个单一的常数谱.This procedure calculates the amplitude spectrum of the wavelet using the autocorrelation of the seismic traces,exactly as in the statistical procedure.The phase spectrum is approximated as a single constant value, using the well logs.这种方法比较稳定, 特别是测井资料与地震数据的相关性较差时.This procedure is more robust than the full phase spectrum calculation, especially when the tie between logs and seismic is poor.计算相位的步骤:(1) 用统计子波提取方法计算子波(不用井资料).(2) 对所提取的子波进行一系列的常相位旋转(3) 用每一次旋转后的子波计算合成道, 并且与地震道进行相关.(4) 选出与地震道产生最大相关值的相位旋转子波提取中的问题:用井提取子波时, 必须首先求出测井曲线之间的最优化相关To extract a wavelet using logs, an optimum correlation must be done first.正确地相关必须以子波已知为前提To perform correlation properly, the wavelet must already be known.实际子波提取的流程:(1) 用统计子波提取来确定一个初步的子波, 假设子波的近似相位已知.(2) 拉伸或压缩测井曲线来标定地震道.(3) 使用新的测井曲线来提取新的子波.(4) 重复第(2)、（3）步，直到提取的子波达到要求为止.。
常用地震子波提取方法简介
常用地震子波提取方法简介[摘要]子波在地震处理和解释中都是一个极其重要的概念，提取制作一个适合地震工区的子波，在作合成地震记录标定及反演工作时都具有极其重要的意义。
[关键词]子波振幅谱相位谱统计子波中图分类号：p315 文献标识码：a 文章编号：1009-914x（2013）09-0161-02子波在地震处理和解释中都是一个极其重要的概念，地震子波是地震记录褶积模型的一个分量，通常指由2至3个或多个相位组成的地震脉冲，确切地说，地震子波就是地震能量由震源通过复杂的地下路径传播到接收器所记录下来的质点运动速度（陆上检波器）或压力（海上检波器）的远场时间域响应。
一个子波可以由它的振幅谱和相位谱来定义，相位谱的类型可以是零相位、常数相位、最小相位、混合相位等；对零相位和常数相位子波而言，可简单将其看作是一系列不同振幅和频率的正弦波的集合，所有的正弦波都是零相位或常数相位的；在频率域中，子波提取问题由两部分组成：确定振幅谱和相位谱，确定相位谱更加困难，并且是反演中误差的主要来源。
通常子波提取方法分为三个主要类型：①、完全确定法：直接应用地表检波器或其它仪器直接测量子波。
②、纯统计法：只根据地震数据确定子波，这种方法很难测定可靠性的相位谱。
通常是为了回避确定子波的困难，使用自相关统计原理，从地震数据中提取地震道的振幅谱来作为子波的振幅谱。
对子波的相位谱则做最小相位或零相位的假设。
③、使用测井曲线法：使用测井曲线提取常相位子波，子波的振幅谱由地震数据自相关获得。
其相位谱假定为一常数，由解一个任意度的最小平方整形滤波器得到平均相位。
使用测井曲线提取全子波，由测井声阻抗计算的反射系数与地震道用最小平方拟合的方法求取子波的振幅谱和相位谱。
此方法的优点在于能够计算一个精确的子波。
但它对井和地震道之间的匹配关系非常敏感。
即使用测井曲线与地震数据结合，理论上这种方法能够提取井点位置精确的相位信息，但问题是该方法要求测井和地震间必须要有良好的对应关系，而将深度域样点转换为双程旅行时的深时转换可能产生不恰当的对应关系，而这种不恰当的对应关系必将影响子波提取的结果。
地震子波的提取及其在最小平方反褶积中的应用
地震子波的提取及其在最小平方反褶积中的应用作者：王超来源：《科教导刊·电子版》2015年第09期摘要本文通过希尔伯特变换法从实际地震资料中提取最小相位子波，将最小相位子波与通过求解托普利兹方程组得到的反滤波因子做褶积，实现了对地震资料的最小平方反滤波，提高了地震剖面的分辨率。
关键词希尔伯特变换地震记录最小相位子波托普利兹方程组最小平方反滤波地震剖面分辨率中图分类号：TN911.7 文献标识码：A为了精细地解释地下地质情况，在实际地震资料处理过程中，常常会对偏移（或叠后）数据体进行修饰性处理，从而提高地震剖面的分辨率、改善其波组特征。
对地震资料做最小平方反滤波就是提高地震剖面分辨率的一种常用方法。
该方法主要分为两个步骤：第一步为地震子波的提取，第二步为最小平方反滤波。
准确提取地震子波是进行反褶积的前提，而反褶积是提高地震资料分辨率的重要手段。
地震子波是一段具有确定的起始时间和有限能量，有限连续长度的信号，它是地震记录中地震的基本单元；一般认为，地震震源激发时所产生的地震波仅是一个连续时间很短的尖脉冲，随着尖脉冲在介质中的传播，尖脉冲的高频成分很快衰减，波形随之增长，变成了一个具有有限频带宽度和一定延续时间的地震子波。
实际的地震子波是一个很复杂的问题，子波是震源波形经过地层滤波和检波器耦合与仪器效应等综合因素共同作用的结果，具体因素的作用会有所不同，对具体勘探目标的影响也不相同。
在地震子波的选取过程中，由于其能量和相位特征较难确定，实际地震资料通常很难满足子波提取的一些假设，使得地震子波提取技术成为制约地震资料处理质量进一步提高的关键因素。
由于本文的最终目的是对地震资料做最小平方反滤波，因此本文选择提取的子波是最小相位子波，采用的是希尔伯特变换法。
反褶积作为地震数据处理中的一个基本的处理环节和提高地震分辨率的一种重要手段，主要是通过压缩地震记录中的基本地震子波以及交混回响和周期多次波的方式，提取反射系数序列；最小平方反滤波，其又称为脉冲反滤波，是一种比较简单但是很实用的滤波方法；最小平方反滤波就是将期望输出设计为尖脉冲，构造托普利兹矩阵方程组，通过解方程最终求解出反滤波因子，然后用求出的反滤波因子与已求出的地震子波做褶积，从而达到对地震记录做最小平方反滤波的目的。
数字处理ch子波提取与子波整形反褶积ppt文档
§ 3.5 子波提取与子波整形反褶积
一、 子波提取：
一、 子波提取： 13．．i.自直相接海关观上法测勘和法探最:小由相于位盐子度波和的温提度取不同，海希 K水o尔l通m伯o常g特o分r变o成v换公两法式层或，W由ol震d—源出 2． 直接被观该测界法面反射海返上回及海陆面地下勘的探检有波差器异，记录下来可以近似作为地震 3．ii.Z 变陆换上法震源子波单一位般圆很外少无记根录，（尽管有井口检波器）； 4i．ii.测井VS资P 料中求井取下子记波录的初至波排除套管波要 度等测求干井有扰资高后料质可及量作井的为声旁子地波波震、；记密 5．iv.同态对法于可控震源，用震源扫描信号与录接收道的相关结果作子波。
设子子设理理波波子论论为波为上上有为可零可限有得得相得长限到到到位度长无无无。，限度限限b长，长长(t的)b的的(反t反(反)b子(子子0(波b)波,波(b0，(，)，1,实)b,实理(实.1际..)论际.际,.处...,b上处.处.(理.n可.理,b理)中)(得中n中只)到)只只取无取取主限主要主长要的要的的有的反有限有子限项限波项，项，，如，实如m如际项m处m项：项理：中：只取主 理 论在理aa上((数论tt))可上字得((可滤aaa(到((得波m无mm到、000限)无)),,反,a长aa限(((褶的m长mm积反0的00和子1反11)反波,)).子,.,...演.，...波..,...,中实,a，a(a(际经(实mam处常m(际0t0)0理用处m中m(零理m)a))只())中相) 取m只位0主)取子, a要主(波的要m。有0的限有1)项,限....，.项., 如a，(如mm项0m：项m：)) a(t利)a利(t用用()a(ZZ(域a域域m(0求求求)m,根a根根0()的,的的am方(方方0 法m法法10可可)可,.以.1.以以.).,.推.,推.推.a.导..(导,导出am出出(有0有利有m限0用限m限长)长)Z长m子域子)子)波求波波与根与与所的所所求方求反求法反子反可子波子以波（波推（只（导只取只出取有取有有限有限限项限长项）项子）波）与所求 利用利的的的Z用关关关域系系系Z求域：：：根求的根方的法方可法以可推以导推出导有出限有的长限关子长系波子：与波所与求所反求子反波子（波只（取只有取限有项限）项） 的123）））关的123123123123）））））））））最最负混系）））关：小大轴合系最最负混最最负混最最负混最最负混相相取相：小大轴合小大轴合小大轴合小大轴合位值位位相相取相相相取相相相取相相相取相子；子子位值位位位值位位位值位位值位位位波波波子；子子子；子子子；子子子；子子：：：波波波波波波波波波波波波反mm：：：：：：：：：：：：00子反反mm==mm反反mmmm0n波00子00子，+00子00==子====m正==0n波0n波反0n，波0n，+波，++，负，+m正m正子反m正反，m正反，反，半反负，负波子负子反负子反子反轴半子反半在波波半子123波半子波子均轴）））波子轴时时在波轴在轴波在波有均在波均间时时间最最负混均时时均时时值有时时有轴间轴间小大合有间间有间值间，间间值轴的轴相相取相的值轴轴值轴，轴决轴，轴的位值位的正位负，的的，决的定的的决的正负半子；子子半决正负决定正于负定正半负半波波波轴轴定半半于定半半Z于轴半轴半：：：；，于轴轴于轴域轴Z；，轴反轴时mmZ；，Z域；，根Z00域时子；，间域==时根域的时0n间波根时-根，+间的n根分间m正-的间的反，n-分处布的-负n分-n分子反处布开n。分半处布处布波子开。始处布轴开。开在。波始向开。均始时时始向始有向间间向值向轴轴，的的决正负定
地震资料处理中子波的有效提取和反褶积计算
地震资料处理中子波的有效提取和反褶积计算孔维新【摘要】提取地震子波是数据分析和处理的关键,采用不同的子波提取计算方法,比较不同算法求取子波同原始子波间的差异,分析其原因和问题所在；另外通过反褶积获取新的地震记录,对比剖面分辨率的变化情况.其主要目的是提供相关处理参数的算法,比较之间的差异和优劣,寻求最佳匹配算法,计算精度是本文的重点.【期刊名称】《工程地球物理学报》【年(卷),期】2011(008)005【总页数】9页(P579-587)【关键词】反褶积;子波;算法【作者】孔维新【作者单位】中国石化西南石油局第五物探大队,湖南湘潭411104【正文语种】中文【中图分类】P631.3地震波在传播中，受到地层介质的滤波作用，使得垂向分辨率极大降低；而反褶积（deconvolution）又称反滤波（inverse filter），是能消除这种滤波作用，尽量恢复激发信号波的处理方法，它是一种通过压缩基本子波来提高地震数据垂向分辨率的处理过程。
在理想情况下，反褶积能压缩子波长度并衰减多次波，最后在地震道上仅仅保留地下反射系数（Yilmaz，1987）［1］。
而反褶积的实现需要知道反滤波因子，而这得先知道其地震子波，但地震子波的准确求取是个难点。
本文就提取子波的不同算法进行介绍，并对这些算法提取的子波进行对比分析，得出Claerbout提取子波法较合适的结论。
同时本文的论证过程也可为实际地震资料处理提供一些计算思路和办法。
目前提取子波的方法主要有直接观测法（只适用于海上勘探）、自相关法（Wold －Kolmogorov法）［2］、Claerbout算法、多项式求根法（无法处理单位圆上的根）、同态滤波法［3］（抗噪声能力弱，但没有苛刻的假设条件，只要求子波的对数谱和反射系数的对数谱序列在时间上可分，本试验采用同态滤波的递归实现方式）以及利用测井资料提取子波（只适应于有井区域，可通过谱除法和最小平方法求解）等方法。
井控子波提取与反褶积及Q值估算
井控子波提取与反褶积及Q值估算1井控子波提取与反褶积子波反褶积是提高地面地震资料分辨率的一种常见方法，但受表层条件、传播路径、上覆地层吸收衰减等因素的影响，往往很难估算出一个准确的地震子波。
VSP资料中包含了丰富的上、下行波信息，由于受噪声影响相对较小，从下行波场中能观测到比较单纯的地震子波。
因此利用VSP资料提取反褶积算子，能够提高地面地震资料处理的质量。
实际地震资料处理前后频谱对比表明，井控子波反褶积可以使地震资料主频提高10Hz左右。
2井控Q值估算与Q值补偿Q值反映了岩石对地震波传播能量的吸收衰减作用，它对地震信号的振幅和相位都有直接影响。
与地面地震资料相比，由于VSP记录只受到一次表层低降速带的衰减作用，下行波能量较强，是用于Q因子提取的理想数据。
利用VSP资料估算Q值的方法很多，最常用的是谱比法，但谱比法要求地震资料具有较高的信噪比，否则计算的Q值误差较大，甚至出现没有物理意义的数值。
子波模拟法是由Jannsen等人在1985年提出的，它通过Q因子扫描使相邻两道VSP初至波形的相关程度达到最大来获取Q值。
由于相关过程中消除了随机噪声的影响，因此该方法具有更好的抗噪性。
另外，在进行最佳Q因子扫描时增加了一个约束条件，将扫描范围限定在实际岩石品质因子的分布范围内，避免了无意义的负值。
大庆喇嘛甸工区Q值计算结果与利用李氏公式进行拟合的结果，二者趋势非常一致。
为方便资料处理，最终采用时变Q因子对地面地震资料进行Q值补偿：地震反射时间小于750ms 时Q值取40，反射时间为750～900ms时Q值取75，反射时间为900～1700ms时Q值取100。
经过Q值补偿，地面地震资料主频提高了15Hz，频带也扩宽了10Hz左右。
尽管高含水油田开发后期地震油藏描述的目标尺度小、精度要求高，但是由于井网密度大，地质、测井和油藏动态信息丰富，通过井控处理、井控解释、井震联合反演、地震约束油藏建模和数值模拟将地震资料与油藏静态和动态资料融合，可以明显改善油藏描述的效果，如井控小断层解释可以识别断距为2m以上的小断层，随机地震反演可以大幅提高厚度在2m以下的储集层预测精度，地震约束油藏建模和数值模拟可以提高油藏建模和历史拟合的精度等。
子波提取方法
1、提取的子波存入Openworks数据库（左图）；2、提取的子波存成Landmark ASCII文 本格式文件（右图）；3、中下图为提取的子波与其振幅谱和相位谱，可以同时显示五个 子波；4、提取子波编辑器（中上图）：可以对提取的子波在相位上与时间延迟上进行处 理
编辑处理列表－－滤波处理 ：利用提取的子波与其他各种子波，对合成记录进行分 时间段滤波处理，以使合成记录更好地匹配地震剖面
相 关 测 试
SeisWell模块提 供了质量控制工具， 三个统计显示工具：1、 常态测试；2、稳态测 试；3、相关测试。这 些统计工具仅仅在子 波被拾取和应用之后 才有效，他们能够帮 助我们分析计算结果 的有效性。
附注
MIRA 初始化文件（.mirarc）包 括以下七个部分： 1、MiraSystem; 2、SynTool; 3、Prints; 4、PrintSetup; 5、CGM; 6、OpenWorks; 7、General Transforms .mirarc文件中的各个部分包含了” 关键字＝值”设置或变量。你可以手 工定制这些设置或变量。
图中最右侧面板为Checkshot—ITT诊断面板，绿色曲线为原始时差曲线；粉色曲线为 校正过的时差曲线；绿色叉号为参考Checkshot标志；红色叉号为激活Checkshot标志。 通过观察选项，可以查看每一个Checkshot诊断面板中各种曲线与标志的含义。
对于每一个Checkshot诊断面板，其中激活的Checkshot标志是可以编辑的，并可以 将编辑结果存储成新的checkshot到数据库中。 Checkshot图形编辑应该注意：1、更新置换速度；2、预置拐点；3、向下广播速度； 4、x方向Snap;5、z方向Snap。
绘图步骤： 绘图步骤： SynTool—File— Print—OK 注意： 1、起始深度； 2、是否比例（垂 直与水平）； 3、重新计算绘图 尺寸。
子波整形提高合成地震记录质量
子波整形提高合成地震记录质量单刚义,韩立国,张丽华,潘保芝吉林大学地球探测科学与技术学院,长春 130026 摘要:阐述合成地震记录的制作原理,利用测井数据计算出反射系数。
根据地震波传播机制,通过三维地震数据统计得到子波,由此计算出原始合成记录。
再根据过井剖面建立褶积模型,对子波进行整形,得到一个随时间变化的地震子波,最终得到与地震数据的波组特征相吻合的合成地震记录。
为精细储层描述打下了坚实的基础。
关键词:合成地震记录;褶积模型;子波整形;反射系数;深时转换;子波提取中图分类号:P 631.4 文献标识码:A收稿日期:2006207211作者简介:单刚义(19702),男,山东莱阳人,硕士研究生,主要从事地震资料的处理和解释研究工作,E 2m ail :ganyi 2shan@yahoo .com .cn .Usi ng W avelet Shap i ng to I m prove Syn thetic Se is m ogram Qua l itySHAN Gang 2yi ,HAN L i 2guo ,ZHAN G L i 2hua ,PAN B ao 2zh iColleg e of GeoE xp lora tion S cience and T echnology ,J ilin U n iversity ,Chang chun 130026,Ch inaAbstract :Syn thetic seis m ogram is the b ridge that connects seis m o logic data w ith w ell logging data ,and is also the basis of the in terp retati on of structu re and reservo ir’s seis m ic litho logy .T he accu racy of syn thetic seis m ogram directly affects the accu rate calib rati on of seis m ic geo logic ho rizon and the accu racy of reservo ir’s litho logy in terp retati on .T h rough th is p aper illu strates the p rinci p le of syn thesizingseis m ogram .T he reflecti on cofficien ts are ob tained w ith the w ell logging data .B y the p rinci p le of seis m ic w aves ,the w avelet is ex tracted from the seis m ic data of th ree di m en si on ,and the raw syn thetic seis m ogram is fin ished .T h rough bu ilding convo lu ti on m odel ,th is paper shapes w avelet ex tracted from seis m ic data ,th is w avelet varies w ith the ti m e .A t last ,the seis m ic litho logy of it is as the sam e of seis m ic data .T hu s it lays a strong foundati on fo r reservo ir fine descri p ti on .Key words :syn thetic seis m ogram ;convo lu ti on m odel ;w avelet shap ing ;reflecti on coefficien t ;dep th -ti m e conversi on ;ex tracti on of w avelet0 引 言测井和地震等资料的有机结合与综合运用是油藏描述中研究地质构造、岩性、岩相、油藏类型及预测有利含油气区等问题的基础。
感应测井自适应反褶积处理与地层参数的提取方法
（ f ）= h b1 （ f ）= h b2 （ f ）= h b!
T（ 1 f） （ f） g b1 T（ 2 f） （ f） g b2 （8）
T（ - T（ 1 f） 2 f） （ f） g b1
"
-
（ z） h（ g（ c z' = " b z - z' ） b z' ）
（2）
由式 （8） 的 Fourier 逆变换可求得相应的空间域滤波器 （ z） 、 （ z） 和 h b! （ z） ， 进而由 h b1 h b2 （ z ）= g b， T
地层中 HRI 仪器测得的视电阻率是随地层深度变化的， 对测井资料进行反褶积处理首先遇到的问题就是如何 选择背景电导率， 自适应反褶积处理是按测量点逐个选 取它的背景电导率。针对 0 . 05 m 的测量间距， 以该
f $ fc 式中， 的第一盲频率的某一频率， 称 f c 是稍小于 g（ b f） 为 g（ 的截止频率； f 0 是渐变段起点处的频率。由于 b f） 深感应的截止频率比中感应的小得多， 为达到中深感应 纵向分辨率之间的相互匹配可以进一步通过设计差值 滤波器的方法， 用中感应的高频成份补充深感应由于盲 频率低而失掉的高频成份， 同时通过调整窗函数中的 f 0 和 f c， 使中深感应响应函数能在一个较大的背景电导率 范围内达到相互匹配， 即通过设计空间频域的中、 深感 应的以及差值滤波器
第 26 卷 ・ 第1期
李曙光等： 感应测井自适应反褶积处理与地层参数的提取方法
・1・
感应测井自适应反褶积处理与地层参数的提取方法
李曙光
摘
杨善德
要
（吉林大学物理系） 李曙光， 杨善德 . 感应测井自适应反褶积处理与地层参数的提取方法 . 测井技术， （1） ： 2002， 26 40 ~ 45 首先对 3 个典型地层模型的高分辨率感应测井 （HRI） 模拟资料进行了自适应反褶积处理， 提高了纵向分辨能力和减 少围岩的影响。为了提取地层的纵向边界和每层的电阻率以便于测井解释， 提出了一种提取地层参数的特征识别方法， 即利用测井曲线和经过自适应反褶积处理得到的测井曲线对地层深度的相对一次导数的绝对值曲线的极大值位置来确 定地层的纵向边界， 取经自适应反褶积处理得到的测井曲线对地层深度的相对一次导数绝对值曲线极小值处的响应值 作为该层的电阻率。对 3 个典型地层模型处理结果表明， 这种特征识别方法是有效的。 关键词： 感应测井 自适应反褶积 地层参数 提取 特征识别
用VSP资料提取子波对地面地震记录作整形反褶积
用VSP资料提取子波对地面地震记录作整形反褶积
王秀槐
【期刊名称】《石油地球物理勘探》
【年(卷),期】1990(000)006
【摘要】选择与地面地震记录接近的 VSP 资料提取子波 b(t),经快速傅氏变换求得B(f),将|B(f)|平滑滤波成|(f)|后作为滤波因子对 B(f)进行滤波,再经反傅氏变换,即可得到与地面地震记录子波最为接近的子波 d(t)。
然后求取零相位整形因子,再对地震记录作整形反褶积及一些保持高信噪比的处理,便可输出高分辨率的零相位地震剖面。
【总页数】8页(P744-750,761)
【作者】王秀槐
【作者单位】大港石油管理局物探公司计算中心
【正文语种】中文
【中图分类】TE13
【相关文献】
1.地震资料处理中子波的有效提取和反褶积计算 [J], 孔维新
2.VSP优化预测反褶积与VSP子波替换法反褶积 [J], 孙哲;刘洋;王静;田洪;苏华;赵前华
3.震源子波反褶积在海洋地震资料处理中的应用 [J], 王红丽;丁在宇;桂德军;黄金连
4.利用VSP子波替换提高地面地震资料的分辨率 [J], 于茜
5.子波零相位反,反褶积与地震记录分辨率的关系 [J], 刘金俊;王修田
因版权原因，仅展示原文概要，查看原文内容请购买。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
§ 3.5 子波提取与子波整形反 褶积
希尔伯特变换法或Wold— Kolmogorov公式
海上及陆地勘探有差异
单位圆外无根
要求有高质量的声波、密 度测井资料及井旁地震记 录
自学并查相关文献完成作业
§ 3.5 子波提取与子波整形反
褶积
据你所知都有哪些相位的信号？？
同态反褶积前 同态反褶积后
第三章 反褶积
§ 3.1 § 3.2 § 3.3 § 3.4 § 3.5 § 3.6 § 3.7 § 3.8
反褶积的概念 最佳维纳滤波 最小平方反褶积 预测反褶积 子波提取与子波整形反褶积 同态反褶积 地表一致性反褶积 反Q滤波及谱白化
§3.7 地表一致性反褶积
Surface Consistent Deconvolution
§3.7 地表一致性反褶积
常假定为最小相位
§3.7 地表一致性反褶积
§3.7 地表一致性反褶积
（分辨率提高、振幅均匀）
地表一致性反褶积效果比较
(a)普通反褶积
(b)地表一致性反褶积
第三章 反褶积
§ 3.1 § 3.2 § 3.3 § 3.4 § 3.5 § 3.6 § 3.7 § 3.8
反褶积的概念 最佳维纳滤波 最小平方反褶积 预测反褶积 子波提取与子波整形反褶积 同态反褶积 地表一致性反褶积 反Q滤波及谱白化
§ 3.6 同态反褶积
一．复赛谱Cepstrum和D特征系统
两边取对数 记 为
傅氏反 变换
该词来自英文Cepstrum，Spectrum的反写 由褶积变为加法！
§ 3.6 同态反褶积
பைடு நூலகம்
§ 3.6 同态反褶积
§ 3.6 同态反褶积
两边取对数
泰勒展 开
傅氏反 变换
复赛谱
§ 3.6 同态反褶积
对复赛谱作高通滤波，通过逆D系统，可得反射系数，实现同态反褶积。
例放大原来各振幅谱的值，使外包线展平为“白色”的宽频谱。 在不改原先的相位谱的情况下，做傅氏反变换，即得谱白化结果。
二、谱白化
在时域中进行时： 1. 先将频谱区间划分成3-4个滤波频段，每个滤波频段由梯形组
成，它的相邻两个边是互补的(加起来为1)。 2. 分完滤波频段之后，用分频档滤波的方法将记录分为四个频档
另外，Q因子本身也是一种很好的吸收属性,它还可在粘油In弹藏el性a描st介ic
述中得到应用。
A质tte衰nu减ation
Q —— 品质因子，quality factor
反Q滤波（Q补偿）——Q Compensation
1、Q因子的物理含义
物理意义
Q因子反映的就是能量被吸收损耗的比率，一个波长λ内，原？地震波
非因果信号
零相位子波比最小相位子波有更高的分辨率！
零相位子波
最小相位子波
§ 3.5 子波提取与子波整形反 褶积
§ 3.5 子波提取与子波整形反 褶积
n
n
n
子波
子波与反子波 时域分布关系
n
非因果的
反子波
§ 3.5 子波提取与子波整形反 褶积
？？？
如何获得？？
找同学上来求！
§ 3.5 子波提取与子波整形反 褶积
消除了漂移，匹配良好
存在漂移的原因？？
第三章 反褶积
§ 3.1 反褶积的概念 § 3.2 最佳维纳滤波 § 3.3 最小平方反褶积 § 3.4 预测反褶积 § 3.5 子波提取与子波整形反褶积 § 3.6 同态反褶积 §一3．.复7 赛地谱C表ep一st致rum性和反D特褶征积系统 § 3.8 反Q滤波及谱白化
因为褶积因子具有逐点变化的时变性，无法求得精确解，需要做如 下近似：
假设在同一地层中，也就是在一段时间 t 中，Q 有确定值，滤波因子h不变。这样，就可得到：
显然，这种假设在 Q ∞时，即没有衰减的情况下是成立的。 注意：对于每一个时间t都要在频率域计算一个积分，计算量是巨大的。
因此有必要采取某种措施，以避免积分运算！
李氏经验公式：
§3.8 反Q滤波及谱白化
3、反Q滤波的实现方法
大地滤波滤波器的振幅谱为
假定大地滤波过程为最小相位，利用Hilbert变换根据振幅谱求取相位谱：
改造后 频谱 令
则反Q滤波器的频谱为
时间域的滤波因子如下：
具有时变性，逐点 变化
Q 补偿后，即未衰 Q 补偿前，即吸收
减前
衰减后
§3.8 反Q滤波及谱白化
？？
当期望输出与子波具有相同的相位（最大）延迟时 ，误差能量最小！！
§ 3.5 子波提取与子波整形反 褶积
子波 期望输出 整形反褶积算子 褶积输出
最佳结果
§ 3.5 子波提取与子波整形反 褶积
子波整形反褶积前后资料的比较
扫描信号的相 关得到接近零 相位的子波
消除了漂移，匹配良好
最小相位
零相位
的能量E与传播中吸收所损耗的能量ΔE之比即为Q。
§3.8 反Q滤波及谱白化
§3.8 反Q滤波及谱白化
一、反Q滤波
1、Q因子的物理含义
2、Q值的估算
教材中有介 绍，自学
实际生产中 常用
具有速度资 料时可使用
常有振幅包络法，频谱比法，常Q扫描法和李氏公式法等。
常Q扫描法：
选一段记录，用不同的Q值从小到大各扫一次做反Q滤波，组成许多 拼起来的图，然后处理人员在这些图上从浅至深选择一个个Q值。
Dave Hale（SEP，1982）将Q(f,t)用Taylor级数展开，设G(f)对应的时域信 号为g(t) ，则有：
泰勒 展开
俞寿朋 《高分 辨率地震勘探》 中对于Q提取及 影响因素有详细 的讨论。
常Q扫描
实例一
实例二
实例二（放大）
3.8反Q滤波及谱白化
二、谱白化(Spectrum Whiltening)
一、反Q滤波
二、谱白化
§3.8 反Q滤波及谱白化
一、反Q滤波 地震资料处理的常规处理步骤之一
弹性能转换成
振幅衰减主要是几何（球面）扩散与吸收。 热能的过程 由为于恢地复层地的震吸波收原作来用的，能地量震，波提经高地分层辨传率播，后处，理能时量必被须衰做减吸损收耗，
频补率偿变，低即。Q补特偿别。是从深滤层波，的分角辨度率讲大，大就下是降反。Q滤波。
• 谱白化是一种展宽频谱的方法。不过它不改变子波的相位谱，是 一种“纯振幅”的滤波过程。
• 谱白化处理可以在频率域中完成，也可以在时间域里进行。
频率域实现方法
时间域实现方法
二、谱白化
• 在频率域中： 1. 可将地震道振幅谱的各尖峰连成没有极小值的外包线， 2. 再将其实行平滑滤波。 3. 再在纵坐标上加适量“白噪成分”，然后求一个倒数，依此倒数比