子波提取方法

合集下载

地震子波提取方法综述

地震子波提取方法综述

地震子波提取方法综述
地震子波是指地震波在地下介质中传播时所经历的各种反射、折射和散射等作用后形成的波形。

地震子波是地震勘探中重要的信息来源,能够提供地下介质的物理特征,如密度、速度、厚度等信息。

因此,地震子波的提取是地震勘探数据处理的重要步骤。

目前,针对地震子波的提取方法主要可以分为时域方法和频域方法。

时域方法主要包括叠加法、全波形反演法、小波变换法等。

叠加法是一种经典的地震子波提取方法,它通过多次叠加同一接收器上的不同地震记录得到地震子波。

全波形反演法则是一种利用弹性波方程直接求解地震子波的方法,它能够提高地震子波的分辨率和准确度。

小波变换法则是一种将地震记录分解成不同尺度和频率的方法,它能够提取地震信号中不同频率的成分,从而得到更加详细的地下介质信息。

频域方法主要包括卷积模拟法、稀疏表示法、奇异值分解法等。

卷积模拟法通过将地震记录与已知的地下介质模型进行卷积,得到地震子波。

稀疏表示法则是一种通过求解稀疏线性方程组来提取地震子波的方法,它能够提高地震子波的信噪比和分辨率。

奇异值分解法则是一种将地震记录矩阵分解成奇异值矩阵和特征向量矩阵的方法,从而得到地震信号中的主要成分。

综上所述,地震子波的提取是地震勘探数据处理中的一个重要环节,不同的提取方法各有优劣,应根据实际情况进行选择和应用。

未来,随着地震勘探技术的不断发展,地震子波提取方法也将不断地进
行改进和创新。

地震子波提取方法综述

地震子波提取方法综述

地震子波提取方法综述
地震子波提取方法综述
地震勘探是一种通过使用反射波测量地下岩层和地质结构的方法,以
便了解地下情况的技术。

能够提取地震信号中的子波,是地震勘探中
非常重要的技术。

下面将综述几种常用地震子波提取方法:
1.匹配滤波
匹配滤波是一种常用的地震子波提取方法,其基本思想是用一个已知
的波形去匹配地震记录中的波形。

匹配滤波的主要作用就是对地震信
号进行滤波增强,提高信噪比。

该方法在提取精细地震子波方面的效
果比较好。

2.小波变换
小波变换是一种将时间和频率相互联结的数学工具。

对于地震子波提
取来说,小波变换能够使原始信号中的各个频率分量得到充分的展开,并且可以将高频噪声和低频信号有效分离,从而提高地震信号的信噪比。

3.奇异值分解
奇异值分解是一种用于分解矩阵的数学技术。

在地震子波提取中,通
过将地震记录矩阵分解成多个低能量层和高能量层,可以得到最佳的
地震子波提取结果。

该方法对于提取高频率的子波有着很好的效果。

4.模拟退火
模拟退火是一种常用的优化算法,用于解决函数优化问题。

在地震子波提取中,使用模拟退火算法可以搜索地震信号的最优解,并提取出更加精细和准确的地震子波。

该方法在提取特定类型的井测距数据中效果比较好。

以上是几种常用地震子波提取方法的综述。

不同的提取方法各有优缺点,需要根据具体情况选择合适的方法进行使用。

在实际应用中,也可以将不同的提取方法进行组合,以达到更好的效果。

子波提取方法

子波提取方法
.mirarc文件中的各个部分包含了” 关键字=值”设置或变量。你可以手 工定制这些设置或变量。
SeisWell模块:新的合成子波提取程序,它能完成如下任务:1、在用户定义参数基础上,比
较多地震道和合成记录道;2、计算通过地震数据的各个时间范围内的各个地震道的品质因素;3、 以一种或多种颜色码显示信噪比,助你快速识别最佳匹配子波位置;4、从你选择的位置提取子波, 并显示子波谱;5、在当前提取子波的基础上,重新计算合成地震记录;6、产生比例绘图和其他统 计框图,使你能够使用它来评价计算结果;7、使用多种技术提取和应用多种子波;8、提取子波并 自动显示它,或将其存储到数据库中。
在测井曲线采集过程中,由于各种因素的
影响,如井壁垮塌、基线漂移、电缆拉伸等, 需要对曲线进行编辑处理,制作出高精度的合 成记录 曲线编辑方法:1、表格编辑;2、块编辑;3、 厚度编辑;4、鼠标编辑 块编辑:在一深度段内对测井曲线作如下调整: 1、新的平均值;2、乘以一个因子;3、上下
曲线编辑--厚度编辑:
合成地震记录的存储:首先存储时深表至数据库,然后存储合成地震记录至数据库。
注意:存储时深表和合成地震记录时,可以存储成激活的,激活的时深表与合成记录可以 直接在SeisWorks中应用.
在一体化解释过程中,SeisWorks2D/3D软件可以直接调用存入数据库中的时深表和 合成记录,但需要将其激活,用来进行层位标定与钻井地质分层的时深转换,并且在 SeisWorks中,可以直接编辑合成记录,再存入数据库中。
Checkshot测量:可以选择一个激活的Checkshot,并可以编辑此Checkshot;还可以选择一个参考 的Checkshot。编辑Checkshot时,有三个选项:1、编辑双程旅行时;2、冻结层间速度3、向下传 播时深对编辑

地震子波提取方法综述

地震子波提取方法综述
通过最大后验估计(maximum posteriori, MAP)或马尔可夫链-蒙特卡罗(Markov chain Monte Carlo,MCMC)模拟方法得到地震子波。
1.5 循环迭代法[7]
方法考虑到地震子波在空间变化的特点, 首先由多道相关法提取初始的地震子波的振 幅谱,然后结合测井资料确定初始地震子波的 相位谱,然后根据离散反演理论迭代求取精细 的井旁地震子波。
统计性子波提取方法是通过地震道自身 来估计子波,又可分为基于二阶统计量和基于 高阶统计量两种方法。二阶统计量方法首先由
1
Robinson(1975)提出,它基于这样的假设,即 地震子波是时不变的,地下的反射时具有白噪 谱的随机序列,则观测到的地震道的自相关就 给定了地震子波的自相关的一个估计,也就是 已知了地震子波的振幅谱,对于子波的相位 谱,则必须给出一定的假设,如假设地震子波 是零相位、最小相位、最大相位,而实际上地 震子波是一种混合相位的,因此,基于二阶统 计量的自相关统计的方法提取的子波也是不 准确的。80年代后期以来,许多学者开始使用 高阶统计方法来解决地震子波估计的问题,这 些方法大多源于60年代发展起来的累积量和 多谱理论[2],T. Matsuoka和T.J. Ulrych(1984) 最早将它们用于混合相位地震子波的估计,高 阶统计量地震子波提取的新思想是由 zear(1993)和D.R.Velis(1996)提出的,他 们将非高斯信号处理中的四阶累积量用于子 波的估计,为解决混合相位子波估计问题提供 了一条全新的思路。
的情况下,仅仅根据系统的输出D来辨识系统
函数W。
输入 R
输出 D
系统函数 W
图1 输入和系统响应未知的盲系统辨识示意图
统计性地震子波提取与盲系统辨识问题 类似,如果将地层反射系数看作输入 R,将系 地震子波看作是系统函数 W,统计性地震子波 提取就是在地层反射系数和地震子波都未知 的情况下,仅仅根据观测到的地震记录 D 来估 计地震子波。

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

常用地震子波提取方法简介
科 学论 坛
I 叠
C h i n a s c i e n c e a n d T e c h n o l o g y R e v i e w
基 于 学 生个 性 化 能 力培 养行 动 导 向在数 控 技 术 实训
教 学 中的 应 用 研 究
高生 祥
( 浙 江工 业职 业 技术学 院 3 1 2 0 0 0 )
2 . 行 动导 向简 介
育对象的气质、 情绪、 认知、 兴趣、 能力 、 性格、 价值观和信念等进行人格整合和
个 性优 化。 通 过应用 行动 导 向法 教学 能够发 挥学生Байду номын сангаас主体 的主动 性 , 调 动学 生的 求知 欲 , 促进 学生个 性化 发展 。 将 六步法行 动导 向应 用在数控 技术 实训 教学 中 ,
具 体教学 实施 过程 如下 :
1) 、 资 讯
主要是 下达典型 工作任务 , 明确任务 的教学 目的 、 重点和难 点, 通过 任务 的 下达 , 让学 生去 自主查 阅完成 任务 的相关 资料 , 培养学 生 的个性化 查 阅资料 的
能力 。
2) 、 计 划
在2 0 世纪 8 哞 代德 国 的双 元制 职业教 育 , 就 是一 种 以“ 实 践为 导 向” 的教 育, 它注 重实践性 教学 环节 , 突 出职 业实践 能力 的综 合培 养 , 其实质 是 : 教学 组 织 以学生为 中心 , 教 师是学 习过程 的组 织者 与协调 人 , 遵 循“ 资讯 、 计划、 决策 、 实施 、 检查 、 评 估 这一 完整 的“ 行动 ” 过 程序列 。 行动导 向的学 习主要探 讨认 知结 构 与个体 活动 之间 的 关系 ; 与之 不 同的 是, 行动导 向理论 以人为本 , 强调 人在实 现 目标过程 中进行反 思的重要 性 , 注重

子波基本理论与提取方法

子波基本理论与提取方法

---------------------------------------------------------------最新资料推荐------------------------------------------------------子波基本理论与提取方法子波基本理论与提取方法 1 地震子波基本原理由震源激发、经地下传播并被人们在地面或井中接收到的地震波通常是一个短的脉冲振动,称该振动为振动子波。

它可以理解为有确定起始时间和有限能量,在很短时间内衰减的一个信号。

地震子波其振动的一个根本属性是振动的非周期性。

因此,它的动力学参数应有别于描述周期振动的振幅、频率、相位等参数,而用振幅谱、相位谱等概念来描述。

子波一般是物理可实现的,特别是地震子波,作为一个物理滤波器的响应函数,自然是物理可实现的,所有必定为非零相子波,但不同子波相位延迟不同。

子波包括最小相位子波、最大相位子波、混合相位子波。

子波的 Z 变换是一个多项式:nnzbzbzbbzB...)(221 若此多项式的全部零点均在单位圆外,则为最小相位子波;在单位圆内,为最大相位子波;零点在单位圆的内外都有,则为混合相位子波。

最小相位最大相位混合相位 2 地震子波的数学模型实际中的地震子波是一个很复杂的问题,因为地震子波与地层岩石性质有关,地层岩石性质本身就是一个复杂体。

为了研究方便,仍需要对地震子波进行模拟,目前普遍认为雷1/ 4克提出的地震子波数学模型具有广泛的代表性,即称雷克子波。

最小相位的地震子波的数学模型为:ftetbat2sin)(2式中:f为子波的主频;)ln(22Mf为子波衰减系数;|/|21mmM 为最大波峰值1m和最大波谷值2m 之对比。

其波形大致如图所示:3 地震子波提取的基本方法地震子波的提取方法有两大类:第一类是确定性子波提取方法;第二类是统计性子波提取方法。

确定性子波提取方法指的是利用测井资料首先计算出反射系数序列,然后结合井旁地震道由褶积理论求出地震子波,它的优点是不需要对反射系数序列的分布作任何假设,能得到较为准确的子波,而统计性方法的优点是不需要测井信息也可以得到子波的估计,但缺点是需对所用的地震资料和地下的反射系数序列的分布进行某种假设,所得子波理论上的精度不是高很。

数字处理-ch3-2子波提取与子波整形反褶积

数字处理-ch3-2子波提取与子波整形反褶积
第三章
反褶积
§ 3.1 反褶积的概念 § 3.2 最佳维纳滤波 § 3.3 最小平方反褶积 § 3.4 预测反褶积 § 3.5 子波提取与子波整形反褶积 § 3.6 同态反褶积 一、 子波提取: § 3.7 地表一致性反褶积 二.3. 子波整形反褶积 直接观测法 § 3.8 反Q滤波及谱白化
i. 海上勘探:由于盐度和温度不同,海水通常分成两 1. 子波的相位与分辨率
1. 子波的相位与分辨率
在振幅谱相同的子波中, 零相位子波的分辨率最高, 最小相位子波次之。 在振幅谱相同的子波中, 零相位子波的分辨率最高, 最小相位子波次之。 在振幅谱相同的子波中, 零相位子波的分辨率最高, 最小相位子波次之。 在振幅谱相同的子波中, 零相位子波的分辨率最高, 最小相位子波次之。
海上勘探和陆上爆炸震源产生的地震子波最接近最小相位,可控震源的
在数字滤波、反褶积和反演中经常用零相位子波。 在数字滤波、反褶积和反演中经常用零相位子波。 在数字滤波、反褶积和反演中经常用零相位子波。 在数字滤波、反褶积和反演中经常用零相位子波。 子波为零相位。
在振幅谱相同的子波中, 零相位子波的分辨率最高, 最小相位子波次之。 可控震源的 海上勘探和陆上爆炸震源产生的地震子波最接近最小相位, 海上勘探和陆上爆炸震源产生的地震子波最接近最小相位,可控震源的 海上勘探和陆上爆炸震源产生的地震子波最接近最小相位, 可控震源的
2110误差能量1621期望输出为010求解方程为21162117期望输出为001求解方程为21202120这与期望输出的误差能量为213535子波提取与子波整形反褶积子波提取与子波整形反褶积2110误差能量1621期望输出为010求解方程为21162117期望输出为001求解方程为21202120这与期望输出的误差能量为212110误差能量1621期望输出为010求解方程为21162117期望输出为001求解方程为21202120这与期望输出的误差能量为212110误差能量1621期望输出为010求解方程为21162117期望输出为001求解方程为212021

谱比法求取地震子波

谱比法求取地震子波

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

子波提取: 子波提取 从井旁地震道提取子波,用提取的子波与反射系数褶积,求取合成地震记录,以使合成
记录与井旁地震道匹配的更好。 子波提取方法: 子波提取方法 1、自相关算法(比较常用);2、维纳-莱文森算法 子波输出类型: 子波输出类型 1、Openworks数据库;2、Landmark ASCII文本格式;3、老式Mira ASCII文本格 式
在一体化解释过程中,SeisWorks2D/3D软件可以直接调用存入数据库中的时深表和 合成记录,但需要将其激活,用来进行层位标定与钻井地质分层的时深转换,并且在 SeisWorks中,可以直接编辑合成记录,再存入数据库中。
SynTool绘图:SynTool软件可以直接生成CGM绘图文件或PS文件,用于绘制 SynTool面板图形,若机器上安装了ZEH或SDI绘图软件,且配置了绘图仪,如HP或 VERSATEC绘图仪,就可以直接绘图了
Vint(Cor-Raw):显示原始的和校正过的时差曲线之间的差作为层速度;3、Vint(Blocked):显示各 Checkshot深度范围内原始的和校正过的时差曲线作为平均层速度;4、ITT:显示原始的和校正过的 时差曲线作为层间传输时间;5、Time(2 Way Travel):显示从原始的和校正过的时差曲线计算出的 总旅行时,总旅行时是通过综合的原始时差曲线和校正过的时差曲线计算出来的;6、Time(1 Way Drift):显示校正过的和原始的综合旅行时之间的差;7、Vave:显示从综合校正过的和原始的时差曲 线计算的平均速度。
SynTool—Extract—SeisWell…启动子波提取程序--初始化地震工区--选择三维 数据体。编辑三维子波参数输入表选项:1、欲扫描的地震道中心线号;2、欲扫描的地震 道中心道号;3、扫描线两边的线数;4、扫描道两边的道数;5、反射系数相关时窗的开 始时间;6、第一个地震相关时窗的延迟时间;7、各地震道相关的数目;8、相关时窗长 度;9、平滑时窗长度。


相 关 开 始 时 间
左图为沿 某一主测线的 各个CDP点与 各个相关开始 时间的信噪比 观察图,图中 黄色叉号为扫 描框图内的品 质因数的最佳 统计匹配位置, 色彩值对应于 品质因数值, 即信噪比值。 利用此图,可 以快速识别最 佳匹配子波位 置。
点击此处,可以前后移动各线对应的信噪比观察 图
Checkshot诊断 诊断:可以将不同类型的Checkshot诊断面板加入到SynTool面板中,然后对 诊断 Checkshot诊断数据进行编辑修改,来校正Checkshot时深对。
Checkshot诊断面板的种类:1、Vint:显示原始的和校正过的时差曲线作为层速度;2、 诊断面板的种类: 诊断面板的种类
左图 为SeisWell 对3D数据 的扫描结果, 图中红圈为 井眼位置, 红叉位置为 全部的品质 因数(或称 吻合度)值 最高的位置。 上图为各个 地震道对应 的最大的品 质因数图; 下图为最大 的品质因数 出现的相关 延迟时间图
道 号
道 号
线 号
线 号
左图为上一页中的上图,即信噪比观察图,图中色彩对应于品质因数值;右图为上一页中 的下图,即延迟时间观察图,图中色彩值对应于延迟时间值。
从你选择的位置提取子波,并显示子波谱;在当前提取子波的基础上,重新计算合成地震记录; 还可以将其存储到数据库中或ASCII文本文件中。
提取子波编辑器:可以对提取的子波在相位上与时间延迟上进行处理。
常态 测试
对 正 态 分 布 的 预 期 的 误 差 值
理论值 实际 值
残余量的归一化值
稳态 测试
统 计 上图为(下图中的上中下四分位数的)79点平滑 平均(理想状态:三条红线为三条平行直线) 统 计 残 差 时 间
Checkshot测量:可以选择一个激活的Checkshot,并可以编辑此Checkshot;还可以选择一个参考 测量: 测量 的Checkshot。编辑Checkshot时,有三个选项:1、编辑双程旅行时;2、冻结层间速度3、向下传 播时深对编辑 Checkshot方法:1、不应用Checkshot;2、层间传输时间;3、速度+层间传输时间;4、速度;5、 方法: 方法 多项式;6、多项式+速度+层间传输时间。层间传输时间是最常用的Checkshot方法,它强制综合 时差曲线精确匹配时深对,与各Checkshot时深对对应的时差样点都要被减去或加上一个常数值。并 可以切除在应用Checkshot后超过某些限度的时差值。
曲线编辑--厚度编辑: 曲线编辑--厚度编辑 --厚度编辑
在一深度段内对测井曲线作如 下调整:1、删除剖面;2、拉 拉 压缩(比较常用);3、替 伸/压缩 压缩 换上部;4、替换下部。 注意:曲线编辑是在深度域进 行的
Checkshot
在应用 Checkshot之前,必 须施加TVD校正 (上一页的时深关 系未加Checkshot, 本页的时深关系加 入了Checkshot)。 Checkshot的目的也 是为了合成记录更 加匹配井旁地震剖 面
图中最右侧面板为Checkshot—ITT诊断面板,绿色曲线为原始时差曲线;粉色曲线为 校正过的时差曲线;绿色叉号为参考Checkshot标志;红色叉号为激活Checkshot标志。 通过观察选项,可以查看每一个Checkshot诊断面板中各种曲线与标志的含义。
ቤተ መጻሕፍቲ ባይዱ
对于每一个Checkshot诊断面板,其中激活的Checkshot标志是可以编辑的,并可以 将编辑结果存储成新的checkshot到数据库中。 Checkshot图形编辑应该注意:1、更新置换速度;2、预置拐点;3、向下广播速度; 4、x方向Snap;5、z方向Snap。
在SynTool面板中加入其他测井曲线,如伽玛曲线、中子曲线、自然电位、视电阻率 曲线等,用于分析层序特征。
在合成记录面板中加入子波和发射系数分量贡献,从而使解释人员能够更准确地分析 各反射界面与地震剖面反射同相轴的对应关系,准确地标定出储层。
在测井曲线采集过程中,由于各种因素的 影响,如井壁垮塌、基线漂移、电缆拉伸等, 需要对曲线进行编辑处理,制作出高精度的合 成记录 曲线编辑方法:1、表格编辑;2、块编辑;3、 曲线编辑方法 厚度编辑;4、鼠标编辑 块编辑:在一深度段内对测井曲线作如下调整: 块编辑 1、新的平均值;2、乘以一个因子;3、上下 漂移;4、赋予一个常数;5、插值
1、提取的子波存入Openworks数据库(左图);2、提取的子波存成Landmark ASCII文 本格式文件(右图);3、中下图为提取的子波与其振幅谱和相位谱,可以同时显示五个 子波;4、提取子波编辑器(中上图):可以对提取的子波在相位上与时间延迟上进行处 理
编辑处理列表--滤波处理 :利用提取的子波与其他各种子波,对合成记录进行分 时间段滤波处理,以使合成记录更好地匹配地震剖面
绘图步骤: 绘图步骤: SynTool—File— Print—OK 注意: 1、起始深度; 2、是否比例(垂 直与水平); 3、重新计算绘图 尺寸。
SeisWell模块:新的合成子波提取程序,它能完成如下任务:1、在用户定义参数基础上,比 模块: 模块
较多地震道和合成记录道;2、计算通过地震数据的各个时间范围内的各个地震道的品质因素;3、 以一种或多种颜色码显示信噪比,助你快速识别最佳匹配子波位置;4、从你选择的位置提取子波, 并显示子波谱;5、在当前提取子波的基础上,重新计算合成地震记录;6、产生比例绘图和其他统 计框图,使你能够使用它来评价计算结果;7、使用多种技术提取和应用多种子波;8、提取子波并 自动显示它,或将其存储到数据库中。 图中的提取选项定义了三个内容: 图中的提取选项定义了三个内容:1、提取子波的存储方式与观察方式;2、地震面板来源选 项;3、合成记录覆盖方式。
合成地震记录的存储 对于制作好的合成记录可以五种方式存储:1、以时间域存入数据库;2、替换数据 库中的合成记录;3、以深度域存入数据库;4、存成ASCII文本文件;5、存成磁盘SEGY 文件
合成地震记录的存储: 合成地震记录的存储:首先存储时深表至数据库,然后存储合成地震记录至数据库。 注意: 注意:存储时深表和合成地震记录时,可以存储成激活的,激活的时深表与合成记录可以 直接在SeisWorks中应用.
相 关 测 试
SeisWell模块提 供了质量控制工具, 三个统计显示工具:1、 常态测试;2、稳态测 试;3、相关测试。这 些统计工具仅仅在子 波被拾取和应用之后 才有效,他们能够帮 助我们分析计算结果 的有效性。
附注
MIRA 初始化文件(.mirarc)包 括以下七个部分: 1、MiraSystem; 2、SynTool; 3、Prints; 4、PrintSetup; 5、CGM; 6、OpenWorks; 7、General Transforms .mirarc文件中的各个部分包含了” 关键字=值”设置或变量。你可以手 工定制这些设置或变量。
相关文档
最新文档