(完整版)地震数据处理第一章：地震数据处理基础
地震数据处理方法(DOC)
安徽理工大学一、名词解释（20分）1、、地震资料数字处理:就是利用数字计算机对野外地震勘探所获得的原始资料进行加工、改进，以期得到高质量的、可靠的地震信息，为下一步资料解释提供可靠的依据和有关的地质信息。
2、数字滤波:用电子计算机整理地震勘探资料时，通过褶积的数学处理过程，在时间域内实现对地震信号的滤波作用，称为数字滤波。
（对离散化后的信号进行的滤波，输入输出都是离散信号）3、模拟信号：随时间连续变化的信号。
4、数字信号：模拟数据经量化后得到的离散的值。
5、尼奎斯特频率：使离散时间序列x(nΔt)能够确定时间函数x(t)所对应的两倍采样间隔的倒数，即f=1/2Δt.6、采样定理：7、吉卜斯现象：由于频率响应不连续，而时域滤波因子取有限长，造成频率特性曲线倾斜和波动的现象。
8、假频:抽样数据产生的频率上的混淆。
某一频率的输入信号每个周期的抽样数少于两个时，在系统的的输出端就会被看作是另一频率信号的抽样。
抽样频率的一半叫作褶叠频率或尼奎斯特频率fN；大于尼奎斯特频率的频率fN+Y，会被看作小于它的频率fN-Y。
这两个频率fN+Y和fN-Y相互成为假频。
9、伪门：对连续的滤波因子h(t)用时间采样间隔Δt离散采样后得到h (nΔt)。
如果再按h (nΔt)计算出与它相应的滤波器的频率特性，这时在频率特性图形上，除了有同原来的H (ω)对应的'门'外，还会周期性地重复出现许多门，这些门称为伪门。
产生伪门的原因就是由于对h(t)离散采样造成的。
10、地震子波：由于大地滤波作用，使震源发出的尖脉冲经过地层后，变成一个具有一定时间延续的波形w（t）。
11、道平衡：指在不同的地震记录道间和同一地震记录道德不同层位中建立振幅平衡，前者称为道间均衡，后者称为道内均衡。
12、几何扩散校正：球面波在传播过程中，由于波前面不断扩大，使振幅随距离呈反比衰减，即Ar=A0/r，是一种几何原因造成的某处能量的减小，与介质无关，叫几何扩散，又叫球面扩散。
地震资料处理复习总结(第1-6章)
《地震勘探资料处理》第一章～第六章复习要点总结第一章 地震数据处理基础一维谱分析数字地震记录中，每个地震道是一个按一定时间采样间隔排列的时间序列，每一个地震道都可以用一系列具有不同频率、不同振幅、相位的简谐曲线叠加而成。
应用一维傅里叶变换可以得到地震道的各个简谐成分；应用一维傅里叶反变换可以将各个简谐成分合并为原来的地震道序列。
连续函数正反变换公式：dt et x X t i ωω-∞∞-⎰=)()(~ 正变换 ωωπωd e X t x t i ⎰∞∞-=)(~21)( 反变换 通常由傅里叶变换得到的频谱为一个复函数，称为复数谱。
它可以写成指数形式 )()()(|)(~|)(~ωφωφωωωi i e A e X X ==式中)(ωA 为复数的模，称为振幅谱；)(ωϕ为复数的幅角，称为相位谱。
)()()(22ωωωi r X X A +=，)()(tan )(1ωωωφr i X X -=（弧度也可换算为角度）离散情况下和这个差不多（看PPT 和书P2-3）一维傅里叶变换频谱特征：1、一维傅里叶变换的几个基本性质（推导）线性 翻转 共轭 时移 褶积 相关（功率谱），P3-72、Z 变换（推导）3、采样定理 假频 尼奎斯特频率，tf N ∆=21二维谱分析二维傅里叶变换),(k X ω称为二维函数),(t x X 的频——波谱。
其模量|),(|k X ω称为函数),(t x X 的振幅谱。
由),(k X ω这些频率f 与波数k 的简谐成分叠加即可恢复原来的波场函数),(t x X （二维傅里叶反变换）。
如果有效波和干扰波的在f-k 平面上有差异，就可以利用二维频率一波数域滤波将它们分开，达到压制干扰波，提高性噪比的目的。
二维频谱产生空间假频的原因数字滤波在地震勘探中，用数字仪器记录地震波时，为了保持更多的波的特征，通常利用宽频带进行记录，因此在宽频带范围内记录了各种反射波的同时，也记录了各种干扰波。
地震数据处理.ppt
（2）波阻抗反演、AVO分析、
方法研究效果只有通过解释才能
谱分解处理技术完善 （3）多波多分量数据处理 （4）井下、井间数据处理流程
真正体现 只有通过解释才能发现问题，才 能买现与地质结合
发展目标是：
（5）时移地震数据处理方法
处理完了，解释也就完了
处理停止了，解释也就停止了
二、常规处理技术的精细处理
七、后记
与谱分解技术联合处理
2003年中油下属 14个油田分公司 共设 173个 老资料重新
处理解释项目
2D 91355 KM
3D 17295 KM2
共投资 17630万，取得了巨大的勘探效益
技术特色：
1、2D 和3D 连片处理
（1）数据规则化处理 （2）处理解释方法针对具体地质目标设计 （3）提高分辨率和改善深层数据品质为两个主要亮点 （4）2D连片一般叠前时间/深度偏移处理
绝大部份工作采用常规处理技术流程,常规处理技术方法成熟
动、静校正-叠加-叠后时间偏移
技术水平 = 技术应用水平 + 精细 + 处理员素质
叠前处理的目的实现同相叠加 (时间对齐、波形一致) 叠加是提高信噪比的最基本最有效手段 叠后时间偏移是近似的，精细只能从偏移策略、算法、参数、速度、 输入数据等方靣入手
技术领域均有重大进展，速度的各向异性研究开始受重视
软硬件环境
1、适应采集数据量的猛增，海上三维作业从 500km2 到 3000km2 ,甚至高 达5000 km2，效率由每日3 km2 到 25km2 , 拖缆由2根到12根。
2、利用高速卫星通信和地面ATM 网络等方式，实现采集实时交互处理与 解释。
多块3D连片处理 多次3D采集对比（时移）处理
地震第1章 地震数据处理基础PPT课件
是连续的，满足傅里叶变换存在条件，则利用上述(1一1)式可以得到其傅
里叶变换分X~ ( w )，其变量 表示圆频率，它与频率 f 之间的关 2πf
为
~
X (w)
，称为地震道x(t)的频谱。
由于傅里叶变换是可逆的，如果已知地震道的频谱
~
X
(w)
，则利用傅
里叶反变换(1-2)式可以得到原来的地震道函数x(t)。
Chapter1 地震数据处理基础
§1.1 一维傅里叶变换及其应用
傅里叶变换是地震数据处理的主要数学基础。它不仅是地震道、地 震记录分析和数字滤波的基础，同时在地震数据处理的各个方面都有着 广泛的应用。例如，在反褶积处理、叠加处理和偏移处理及地震波场的 分析中也都有其重要的应用。
一、一维傅里叶变换及频谱分析
或简写成 xn{x0,x1,x2,,xN 1}
由（1-9)式，此x(n△t)的频谱为
X ~(m f)tN1x(nt)ei2π m f nt (m,n=0,1,2,···，N-1) n0
令
Zei2eiwt
令△t=1，则 Zei，w 这时x(n△t)的频谱
~
X
(m可f )以写成以Z为变量
积。这样就可以将函数在时间域复杂的褶积运算通过傅里叶变换变成为其频谱在
频率域的乘积运算。这种转变在地震数据处理过程中是非常方便和有用的。
6、相关
如有一函数x(t)的傅里叶变换为
~
X
(w
)
，其自相关函数
rxx() x(t)x(t)dt
(1-22)
的傅里叶变换为
~
Rxx() rxx()eiwd [ x(t)x(t)d]teiwd
简谐成分 x()在t=0时的初始相位。
《地震勘探数据处理》课程教学大纲
本科生课程大纲课程属性：公共基础/通识教育/学科基础/专业知识/工作技能，课程性质：必修、选修一、课程介绍1.课程描述：该课程授课对象是勘查技术专业的本科生。
课程主要包括地震资料数据处理基础，常规处理的方法原理：地震资料数据处理基本流程、信号处理基础、反褶积、动校正、静校正、速度分析、水平叠加、偏移归位，并包括目标处理：高分辨率地震资料处理、三维地震资料处理、叠前深度偏移处理、井孔地震资料处理、多波多分量地震资料处理、时间推移地震资料处理和特殊处理：亮点及AVO分析、高精度反演、地震属性分析、相干体数据处理、可视化数据处理等内容。
2.设计思路：地震勘探是钻探前勘测石油与天然气资源的重要手段，在煤田和工程地质勘查、区域地质研究和地壳研究等方面，也得到广泛应用。
地震勘探数据处理具有鲜明的跨学科特点且发展迅猛。
因此，该课程除了介绍地震勘探数据常规处理技术外，还会将当前该领域的最新的研究成果分专题进行介绍，每种处理方法都结合大量地震勘探处理实例加以说明，力争深入浅出、通俗易懂，让学生掌握现代地震勘探数据处理的基本理论和方法。
3. 课程与其他课程的关系：地震勘探数据处理为专业知识教育层面课程，是勘查技术与工程专业最基础的专业课程群组成部分，选修学生应具备的扎实数理基础和高水平计算机应用能力及掌握现代地球物理勘探、地球物理数据处理基础等的基本理论、方法和基本技能。
先修课程为：地震勘探原理或勘探地震学。
- 1 -二、课程目标课程主要包括地震资料数据处理基础，常规处理的方法原理：地震资料数据处理基本流程、信号处理基础、反褶积、动校正、静校正、速度分析、水平叠加、偏移归位，并包括目标处理：高分辨率地震资料处理、三维地震资料处理、叠前深度偏移处理、井孔地震资料处理、多波多分量地震资料处理、时间推移地震资料处理和特殊处理：亮点及AVO分析、高精度反演、地震属性分析、相干体数据处理、可视化数据处理等内容。
使学生初步掌握地震资料数据处理的基本原理及相关处理软件等方面的知识；达到基础研究、应用基础研究科学思维和科学实践训练目标；了解该课程以及相关课程的研究现状和发展动态；具有一定的地震资料数据处理研究与开发的能力。
地震数据处理
（1）速度的用途
地震勘探的各个环节都要用到速度信息 1．野外采集：设计观测系统，确定组合检波形式 2．资料处理：动、静校正，滤波，偏移… 3．资料解释： 速度资料对地震勘探的各个环节都会产生影响，最终影响解释的精度，因 此提取、分析、利用速度是地震资料解释的重要环节。
5、时深转换
经水平叠加后，剖面已变成与地质构造特征相对应t0时间剖面。下步工作 就是将其转化为深度剖面。
对浅层畸变大的波形切除示意图
4、速度分析
速度参数在反射法数据处理至关重要。
影响地震时间剖面的质量； 最终影响到地质解释的精度 影响层速度及平均速度的计算精度。
目的：
第一：为水平叠加、偏移等提供处理的速度参数；
第二：为时深转换提供平均速度。 速度分析常采用：速度谱分析，速度扫描。
几种速度和时深转换 （1）真速度 （2）层速度 （3）平均速度 （4）均方根速度
H
或
H
t V
i 1 i
n
i
1 t 0V 2
式中：Vi－地层的层速度，ti－该地层的单程旅行时间；t0－所求深度 处的回声时间， V －平均速度。
6、偏移处理基础
偏移的概念 偏移归位、偏移成像、波场延拓成像等 射线偏移的原理 波动方程偏移原理 从原理角度： 射线偏移 波动方程偏移
预处理：原始记录数据处理之前所必须完成的工作。 目的：把原始数据进行初步加工，使之满足处理方法技术的要求。
包括：剪辑处理、切除、抽道选排。 1、剪辑处理 剪辑：挑选信噪比低的不正常记录道或炮，将其充零。 不正常道：工作不正常道、死道、极性反转道。 不正常记录：外界干扰背景严重而引起的噪声记录，应将整张记录充零。 2、切除 (1)切除强振幅的初至波，这些初至波一般是直达波和浅层折射波等干扰波； (2)切除发生相位畸变的浅层宽角反射波； (3)切除震源干扰波、相干干扰波。
地震数据处理 第一章：地震数据处理基础
3.速度分析(velocity Analysis)；
4.动校正(Normal Moveout Correction)消除由于炮检距不同引起同一
反射波达到时间的差异;
5.叠加(Stack)； 6.显示叠加剖面 (Display) (有波形、变面积、波形+变面积三种显示方式)；
从波形可看出波的振幅、周期、频率等动力学特点；从变面积的角度，它又突出了 反射层，较直观地反映地下构造形态的特点
ICTFT
f (t )
时 域 恢 复 时 域 抽 样
LT
F ( s)
S j j S
F ( j )
截 取 主 周 期 频 域 周 期 延 拓
ILT
j j n F ( e ) f ( n ) e n- DTFT : j j n f ( n) 1 F ( e ) e d 2
地震波不是简谐波，从波剖面中可得到相邻两峰或谷 间的距离称为视波长，其倒数为视波数。
地 震 波 场
地 震 波 场 时 间 切 片， 即 波 动 图
一ቤተ መጻሕፍቲ ባይዱ付里叶变换
一个正弦运动要用频率、振幅和相位才能完整 的描述。
在计算机中用快速算法实现付里叶变换（FFT）。
付里叶变换：
正变换：时域信号 分解 频域信号；
时 间 （s）
频率（Hz）
图1.1－11 几个没有相位延迟但峰值振幅相同的正弦波的总和产生一个带限对称子波， 表示在右边一道上(由星号标出)，这是一个零相位非对称子波
图1.1—12表示给在图1.l-11中的各正弦 波一个线性相位移所产生的结果。线性相 位移在频率域定义为：
时 间 （s）
模拟与数字信号 一道地震信号是一个连续的时间函数。在地震记录中，连续（模拟） 的地震信号在时间域按照固定的比例取样，叫做采样间隔。典型采样间 隔范围在1到4ms，高分辨率要求采样间隔小到0.25ms。 一般地说，给定采样间隔 ，则可恢复的最高频率为尼奎斯特（Niquist） 频率。公式如下：
地震资料数字处理
MB1 MB2 MB3
四、处理流程
（1） 观测系统定义 （2） 野外静校正 （3） 线性动校正 （4） 叠前去噪 （5） 反褶积 （6） 速度分析 （7） 动校正、切除 （8） 剩余静校正 （9） 叠加 （10）叠后去噪 （11）绘图
1、定义观测系统 运行如下模块: 3D Land Geometry Spreadsheet* 弹出如下菜单:
横线方向覆盖次数： N y P * R /(2d )
P 排列不动所需的激发点数；R 接收线数； d 束线间接收线移动距离相当的激发点数。
例 如 ： 对 于8线8炮 制 的 采 集 方 式 ， 其R和P分 别 为8； 束线间接收距为4 200m 800m，横线炮距100m, d 800/100 8， 故 N y 8 8 /(2* 8) 4
10。用MB1击,则 将数据输入表中。
● Sources (填写炮点参数表）
Source Line Station St Index x y z FFID Time Date Offset Skid Shot fold* Pat Shift Static
用户定义的震源编号 震源线号 震源站号 同一炮点识别器(1~9) 震源点的X坐标 震源点的Y坐标 震源点的高程 野外文件号 放炮时间 放炮日期 炮点垂直炮线的偏移距,+右、-左 炮点平行炮线的偏移距,+大、-小号 接收道数 排列滚动的站点值 用户定义的静校量
● Receivers (填写接收点参数表）
Station Pt index x y Elev Line Static
接收点站号(桩号) 站点识别器(1~9) 接收点x坐标 接收点y坐标 接收点高程 接收点线号 接收点静校量
本例为498~825/线 全1 为测量数据 为测量数据 为测量数据 本例线号为1~9 全0
地震勘探-地震数据处理
第三节 水平叠加(stacking)
在地震资料数字处理中,水平叠加是常规处理 方法中最基本、最必要的一环。
叠前：（静校正、动校正）
水平叠加
叠加：（常规、自适应）
叠后：（道内均衡、道间均衡）
第一部分 叠前
（1）动校正处理 （消除炮检距引起的时差）
（2）静校正处理 （消除地表因素引起的时差）
2、为什么要进行地震资料处理
野外地震资料中包含着有关地下构造和岩性的信息，但这些信 息是叠加在干扰背景上且被一些外界因素所扭曲，信息之间往 往是互相交织的，不宜直接用于地质解释。因此，需要对野外 采集的地震资料进行室内处理。
地震资料处理
野外地震记录
处理后地震记录
3、地震处理三个基本阶段
（1）预处理：将野外采集数据转换成适合计算机 处理的格式，并对数据作相应编辑和校正。
第二节 预处理
预处理概念：
指在对数据作实质性处理之前为满足计算机和软件系统以及 处理方法的要求，对输入的原始数据所必须完成的一些准备工作。
预处理的主要内容：
（1） 数据解编 （2） 道编辑 （3） 增益恢复 （4） 抽道集 （5） 初至切除
1．数据解编（重排）:将按时序分道排列的野外 磁带记录数据转换成按道分时排列的形式（实质 上就是矩阵的转置）。野外磁带数据是按时序排 列的，即依次记下每一道的第一个采样值，各道 记完后，再依次记下各道的第二个采样值，由此 类推，直至结束。矩阵表示如下：
A11 A12
A21 A22
A1N AN2
A1M A2M ANM
A11 A21
A12 A22
地震数据处理【精选文档】
地震数据整体流程不同软件的地震数据处理方式不同，但是所有软件的处理流程基本是固定不变的,最多也是在处理过程中处理顺序的不同。
整体流程如下：1 数据输入（又称为数据IO）数据输入是将野外磁带数据转换成处理系统格式，加载到磁盘上，主要指解编或格式转换。
解编：将多路编排方式记录的数据（时序)变为道序记录方式，并对数据进行增益恢复等处理的过程.如果野外采集数据是道序数据，则只需进行格式转换，即转成处理系统可接受的格式.注:早期的时序数据格式为记录时先记录第一道第一个采样点、第二道第一个采样点、……、第一道第二个采样点、第二道第二个采样点、……直至结束。
现在的道序记录格式为记录时直接记录第一道所有数据、第二道所有数据、……直至结束,只是在每一道数据前加上道头数据.将时序数据变为道序数据只需要对矩阵进行转置即可。
2 置道头2.1 观测系统定义目的为模拟野外，定义一个相对坐标系，将野外的激发点、接收点的实际位置放到这个相对的坐标系中。
即将SPS文件转换为GE-Lib文件，包括1)物理点间距2）总共有多少个物理点3)炮点位置4）每炮第一道位置5）排列图形.2.2 置道头观测系统定义完成后,处理软件中置道头模块，可以根据定义的观测系统，计算出各个需要的道头字的值并放入地震数据的道头中。
当道头置入了内容后，我们任取一道都可以从道头中了解到这一道属于哪一炮、哪一道？CMP号是多少？炮间距是多少?炮点静校正量、检波点静校正量是多少？等等.后续处理的各个模块都是从道头中获取信息,进行相应的处理,如抽CMP道集，只要将数据道头中CMP号相同的道排在一起就可以了。
因此道头如果有错误，后续工作也是错误的。
GOEAST软件有128个道头，1个道头占4个字节，关键的为2(炮号）、4(CMP号)、17（道号)、18（物理点号）、19（线号)、20（炮检距）等。
2。
3 观测系统检查利用置完道头的数据,绘制炮、检波点位置图、线性动校正图.3 静校正（野外静校正）静校正为利用测得的表层参数或利用地震数据计算静校正量，对地震道进行时间校正，以消除地形、风化层等表层因素变化时对地震波旅行时的影响.静校正是实现共中心点叠加的一项最主要的基础工作。
地震数据处理讲课PPT课件
地
30年代
由折射地震法改进为反射法
震
50年代
出现多次覆盖技术
勘
60年代
出现数字地震仪及数字处理技术
探 发
70年代初
偏移归位成像技术
展
70年代中期 三维地震勘探出现
90年代
高分辨率与三维地震结合
一、引言
数据处理的主要机型及软件：
微机 — 简单常规处理，主要用于野外质量监控。 工作站、服务器 — 具有单个或多个CPU，处理能力较强， 能满足常规处理及特殊处理。 微机集群 — 含有多个节点，具有较强处理能力和速度， 主要用于大数据量三维处理、三维连片处理及迭前时间域或深 度域偏移迭加。 兰德马克公司的PROMAX处理系统 帕拉代姆公司的FOCUS处理系统 西方地球物理公司的OMEGA处理系统 法国CGG公司的CGG处理系统 中石油球公司的GRISYS处理系统
瑞利（Rayleigh)准则：一个反射波的分辨率的极限是1/4波长。
横向分辨率不常用，而且在迭偏剖面上很难讨论清楚。
一、引言 处理技术涉及的基本概念
3、“三高”处理 三高是指高保真度，高信噪比，高分辨率。这是我们在数
据处理中应当严格执行的一条准则。 高保真度是指处理的成果剖面主要目的层构造特征清晰、断
三、常规处理流程
二维地震资料处理流程
三维地震资料处理流程
三、常规处理流程
地震数据输入— 解编
由于野外数字地震仪记录的地震数据是按时序排列的，也 就是说在磁带上记录的地震数据是按采样间隔以时间序列排列， 既依此记录每道的第一样值，各道记完后，再依此记录以下各 道的第二个样值，依此类推。目前在地震数据处理中多半是采 用单道或多道处理，因此要求地震数据是按道序排列，这种将 野外磁带上按时序排列的数字样值，转换为适合计算机处理的 按道序排列样值转换过程，称为数据解编或数据重排。
地震数据处理
1：我们要新建一个Project，然后点左上角的加号加载数据进来。
选择加载数据的数据类型为“segy”，对加载的数据进行重命名。
如下图：3：选择数据名前边的曲线图标，展现出原始的曲线,该曲线为由模型得到的,模型由七层构成。
第一层的倾角为0，以后每层的倾角逐渐增大，深度也逐渐增大，第四层的深度和第五层的深度相差很小。
4：进行观测系统的建立，选择左上角的“interactive”下的“geometry widow display”，自己重新建立一个观测系统，注意不要load观测系统，得到如下图：5：进行观测系统的参数的设置，增量参数设置好了以后，点击软件上的加号，添加第一个炮点的参数信息，如下图所示6：加好第一行的信息以后，按住加号不放，选择添加多行的炮点信息，参数如下表所示：7:最后得到了，炮点的观测系统的信息，如下图：8：点击左上角的小旗子进行接收点的观测系统的设置，和炮点的设置方法类似也是先设置增量，然后设置第一行的参数，参数设置情况见下面两图：9：然后点击计算CMP和检波器叠加次数和炮检距离，以便写入道头和相应的位置，完成后续的处理，可以点击查看检波器的叠加次数，并可根据这个判断建立的观测系统是否正确。
\可以看出，CMP的最大覆盖次数为六次，根据已知的模型参数信息，由公式可以计算出理论上的最大覆盖次数为六次，因此可以说明上图建立的观测系统是正确的，可以按住保存钮不放，在出现的按钮组合中选择，将相关的信息写入道头，点击这个按钮还可以选择将炮检关系及他们的坐标写入文本文件以便以后查看，。
加入观测系统后(主要是将相关的信息写入道头)，便可对数据的输入按所要的处理模块选择不同排列方式。
10：至此，观测系统已经建立完成，下面进行速度分析，选择菜单栏的“Job Flow”下的“Open flow file+command window”，把模块拖拽进来如下图所示：设置为5，即抽取以选择的CMP为中心，左右各5道共十一个cmp道集如下图所示进行连线进行设置，先建立可选的CMP道集，设置参数如下图5000，50，速度谱法的速度扫描范围为500，5000，50，如下面图所示：对输出的文件名进行设置，文件名可以任意设置点击左上角菜单栏的最后一个图标运行上面的工作流，执行完后得到三个输出的数据集分别是速度谱，最佳炮检距道集，常速度叠加数据集。
地震数据处理过程及格式说明
向相反取负值) （分米） 。
2
PDF 文件使用 "pdfFactory Pro" 试用版本创建
11 21—22
12 23—24 13 25—26 14 27—28 15 29—30 16 31—32 17 33—34 l8—l 35
18—2 36
19 37—38 20 39—40 21 41—42 22 43—44
45—48 炮点的地面高程（cm）。 49—52 炮井深度(正数，cm)。 53—56 接收点基准面高程（cm）。 57—60 炮点基准面高程（cm）。 61—64 炮点的水深（cm）。 65—68 接收点的水深（cm）。 69—70 对 41 一 68 字节中的所有高程和深度应用此因
子给出真值。比例因子=l，土 10，土 100，土 1000 或者 土 10000。如果为正，乘以因子；如果为负， 则除以因子。(此约定中= -100) 71—72 对 73—88 字节中的所有坐标应用此因子给出真 值。比例因子=1，土 10，土[00，土 1000 或者 土 10000。如果为正，乘以因子；如果为负，则 除以因子。(此约定中= -10) 73—76 炮点坐标—X（分米）。（如果坐标单位是弧度·秒， 77—80 炮点坐标—Y（分米）。X 值代表径度，Y 值代表 81—84 接收点坐标—X（分米）。纬度；正值代表格林威 85—88 接收点坐标—Y（分米）。治子午线东或者赤道北的
4＝其他 34—2 68 135—136 扫描道起始斜坡长度，以 ms 表示。 35—1 69 137—138 扫描道终了斜坡长度，以 ms 表示。 35—2 70 139—140 斜坡类型：1＝线性；2＝COS²；3＝其他 36—1 71 141—142 滤假频的频率(如果使用) 36—2 72 143—144 滤假频的陡度 37—1 73 145—146 陷波频率(如果使用) 37—2 74 147—148 陷波陡度 38—1 75 149—150 低截频率(如果使用) 38—2 76 151—152 高截频率(如果使用) 39—1 77 153—154 低截频率陡度 39—2 78 155—156 高截频率陡度 40—1 79 157—158 数据记录的年 40—2 80 159—160 日 41—1 81 161—162 小时(24 小时制) 41—2 82 163—164 分 42—1 83 165—166 秒 42—2 84 167—168 时间代码：1＝当地时间；2＝格林威治时间；
(完整版)地震资料数字处理复习题
地震资料数字处理复习题一、名词解释〔20分〕1、速度谱把地震波的能量相对于波速的变化关系的曲线称为速度谱。
在地震勘探中，速度谱通常指屡次覆盖技术中的叠加速度谱。
2、反滤波又称反褶积，是指为提高纵向分辨率，去掉大地滤波器的作用，把延续几十至100ms的地震子波b〔t〕压缩成原来的震源脉冲形式，地震记录变成反映反射系数序列的窄脉冲组合。
3、地震资料数字处理就是利用数字计算机对野外地震勘探所获得的原始资料进行加工、改良，以期得到高质量的、可靠的地震信息，为下一步资料解释提供可靠的依据和有关的地质信息。
4、数字滤波数字滤波就是指用数学运算的方式用数字电子计算机来实现滤波。
对离散化后的信号进行滤波，输入、输出都是离散数据。
5、水平叠加将不同接收点受到的来自地下同一反射点的不同激发点的信号，经动校正叠加起来。
6、叠加速度在一般情况下，都可将共中心点反射波时距曲线看作双曲线，用一个同样的式子来表示：t2=t2+x2/V2，其中，V 就是叠加速度。
0 αα7、静校正把由于激发和接收时地表条件变化所引起的时差找出来，再对其进行校正，使畸变了的时距曲线恢复成双曲线，以便能够正确地解释地下的构造情况，这个过程叫做静校正。
8、动校正消除由于接受点偏离炮点所引起的时差的过程，又叫正常时差校正。
9、假频一个连续信号用过大的采样得到的离散序列实际上含有连续信号中高频成分的奉献。
这些高频成分折叠到离散时间序列中较低的频率。
这种现象是由连续信号采样缺乏引起的，称作假频。
10、亮点技术所谓“亮点〞狭义地说是指地震反射剖面上由于地下油气藏存在所引起的地震反射波振幅相对增强的“点〞。
利用地震反射波的振幅异常，同时也利用反射波的极性反转、水平反射的出现、速度的降低及吸收系数的增大等一系列亮点标识综合指示地下油、气藏的存在，进而直接寻找油、气藏的技术。
11、相关定量地表示两个函数之间相似程度的一种数学方法。
12、自相关表示波形本身在不同相对时移值时的相关程度。
地震数据处理-知识点
第一章概述1.1 地震数据处理的目的是对地震采集数据做各种处理提高反射波数据的信噪比、分辨率和保真度以便于解释。
地震数据处理主要包括地震反褶积、叠加和偏移成像三大技术。
地震反褶积是通过压缩地震子波提高地震时间分辨率；叠加的目的是压制随机噪声提高地震信噪比；偏移成像包括射线偏移和波动方程偏移两大类，主要目的是实现反射界面的空间归位和恢复反射界面空间的波场特征、振幅变化和反射系数，提高地震空间分辨率和地震保真度。
1.2地震数据处理包括预处理、常规处理和特殊处理三个阶段。
常规处理包括反褶积、叠加和偏移三大技术。
预处理是把野外数据格式转换成适合计算机处理的格式并对数据做相应编辑和校正。
它包括数据解编、格式转换、编辑、几何扩散校正、建立野外观测系统和野外静校正等。
数据解编：把按时分道的数据记录方式变换成按道分时的数据记录方式。
道编辑：噪音道、带有瞬变噪音的道或单频信号道都要删除；极性反转的道要改正。
几何扩散校正：通过给数据加一增益恢复函数，以校正波前（球面）扩散对振幅的影响。
野外静校正：对路上资料，把所有炮点和接收点位置均校正到一个公共基准面上，以消除高程、低降速带和井深对旅行时的影响。
反褶积的基础是最佳维纳滤波。
特殊处理主要包括T-P变换、小波变换、三维叠前深度偏移、子波处理、属性分析和反演等。
T-P变换：将偏移距-时间域变换到射线参数-截距时间域，可用来压制面波和多次波。
小波变换：小波变换与多尺度分析可用于去噪、数据压缩、提高分辨率处理、信号增强和解波动方程等。
第二章数字滤波2.1 滤波器可以分为模拟滤波器和数字滤波器采样定理时域实参数的滤波器，其频率振幅谱是偶对称的，而相位谱是奇对称的。
一个滤波器如果是稳定的，这是指当输入信号为有限信号时，其输出也是有限信号。
最小相位，在时间域中也称最小能量延迟，在频率域则常称为最小相位滞后。
纯振幅滤波器也称零相位滤波器。
又称为理想滤波器。
2.2 理想滤波器常设计成四种类型：低通滤波器、带通滤波器、带陷滤波器和高通滤波器。
地震波处理
精选ppt
2
• 二、编辑
• 在地震数据采集中，由于施工 现场复杂，外界干扰大，难免出 现一些不正常道和共炮点记录， 这些记录信噪比低，如果参与叠 加处理会严重影响处理效果。
精选ppt
3
▪在正式处理之前，需要对这些不正常 的记录进行编辑处理，例如对信噪比很 低的不正常道进行充零处理，发现极性 反转的工作道对它们进行改正等。
另外，还要显示有代表性的记录并 观察初至同相轴，以便进行初至切除。
精选ppt
4
切除是为了消除包括噪声的记录 开始部分所存在的高振幅，这样 做对避免以后处理时出现的叠加 噪声有好处。切除的方法就是用 零乘需要切除的记录段。
•
精选ppt
5
三, 抽道集
抽道集也叫共深度点 选排，是把具有相同 炮检距点的记录道排 成一组，按共深度点 号次序排在一起。抽 道集处理后，磁带上 记录的次序是以共深 度点号（CDP)为次序 的记录，以后所有的 处理都将方便地以共 深度点格式进行。
16
• 由于地下介质的复杂性，故参数提取 是一个十分复杂而艰巨的任务，只能 用一些简化的方法和近似的假设条件 来求取。
• 所用方法不同，可得到不同定义的速 度参数。地震速度分析中普遍采用速 度谱分析和速度扫描技术，得到平均 速度、均方根速度、层速度等速度参 数。
•
精选ppt
17
设共反射点道集内有N个记录道， 其炮检距分别为x1，x2，…xn ，各个 记录道对应的正常时差分别为t1，t2 ，…tn ，炮检距为xi的第i个记录道的 正常时差为
精选ppt
48
地 震 道 的 记 录 样 式
精选ppt
49
结果第一个样值就是原记录道上的第 五个样值，然后把尾部的四个单元冲 零。
基于小波分析的地震数据处理技术研究
基于小波分析的地震数据处理技术研究第一章引言地震是地球内部能量释放所引起的自然灾害之一，也是一种让人们感受到地球活力和自然力量的重要现象。
随着地震测量技术的发展，人们获得了大量的地震数据，可是如何对这些数据进行分析处理一直是一个难题。
小波分析作为一种新兴的信号处理方法，已被广泛应用于地震数据处理领域，成为一种有力的手段。
本文将就小波分析在地震数据处理中的应用进行探讨，以期为地震数据处理提供一些新思路和新方法。
第二章小波分析原理及方法2.1 小波分析概述小波分析（Wavelet Analysis）是一种非平稳信号分析方法，与傅里叶分析相比，小波分析具有更好的时域局部性和频域紧凑性。
小波分析实质是一种时频变换方法，能够将信号分解成不同频率的成分，以及在不同时间尺度上的成分。
2.2 小波分解小波分解是指将原始信号分解成一组特定的基函数（小波函数），然后根据这组小波函数的系数来重构信号。
小波分解有两种方法，分别是连续小波变换（CWT）和离散小波变换（DWT）。
CWT是在连续时间域上对信号进行小波变换，而DWT是在离散时间域上进行小波变换，DWT由于运算速度快，被广泛应用于实际工程中。
2.3 小波重构小波重构是指利用小波分解的系数来重构原始信号。
重构时主要采用级联重构和逐层逆向重构等方法。
第三章小波分析在地震数据处理中的应用3.1 小波包分析小波包分析是指在DWT的基础上对小波系数再进行分解，从而得到能够更好地描述信号局部时序特性的子带。
小波包分析可用于地震数据短时序特征提取、识别和分类等方面。
3.2 小波分析在地震数据压缩中的应用地震数据一般具有高维、高精度、大容量等特点，这要求将地震数据进行压缩处理。
小波分析作为一种有效的压缩方法，已被广泛应用于地震数据处理中。
3.3 小波分析在地震数据去噪中的应用地震数据中常常存在着各种干扰噪声，这会干扰地震数据的分析处理。
小波分析可以通过对地震数据进行小波分解，然后对不同频率小波系数进行滤波，从而实现地震数据的去噪。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
学时：64（理论48，实习16） 学分：4（必修课）
绪论
地震资料野外采集 地震勘探三个环节： 地震资料数据处理
地震资料地质解释
光点技术 地震勘探数字化进展： 模拟记录
数字化技术
地震资料处理
提高信噪比 提高分辨率 提高保真度
为目的。
Why is Do Data processing?
n 0
1
N 1
F (n)e j(2 / N
n )k
n
N n0
图1.1-2 标有星号的时域信号（道)，可用一组具有 不同频率、振幅 和相位延迟 的正弦运动来表示。
也可以用振幅谱和相位谱来表示：
1.1-3 图 标有星号的时域信号的振幅谱（下）和相位谱（上） 振幅谱上的每一点相当于该频率正弦曲线的峰值振幅。 相位谱上的每一点相当于该频率正弦曲线的波峰或波谷相对于t＝0的时间延迟。
1
j2
j
F
(
s
)e
st
ds
j
CTFT
:
F () f (t)
1 2
f (t )e j tdt F ()e j
t d
f (t)
时
时
域
域
恢
抽
复
样
f (n)
CTFT
ICTFT
LT
S j
F (s)
ILT
j S
z e ST
将不容易 解释的原始 资料变成容 易解释的时
间剖面；
1.预处理(Preprocessing) (解编,不正常道、炮的处理,抽道集) 2.静校正(Static Correction)消除表层因素(低降速带厚度、速度变化、地
表起伏不平)造成的时差影响;对同一道而言,从浅到深,有相同的校正量,故称静校正。
3.速度分析(velocity Analysis)； 4.动校正(Normal Moveout Correction)消除由于炮检距不同引起同一
反射波达到时间的差异;
5.叠加(Stack)； 6.显示叠加剖面 (Display) (有波形、变面积、波形+变面积三种显示方式)；
从波形可看出波的振幅、周期、频率等动力学特点；从变面积的角度，它又突出了 反射层，较直观地反映地下构造形态的特点
7.数据滤波和反滤波(Filtering and AntiFiltering)； 8.偏移归位处理(Migration Processing)
偏移：
是通过数值计算把地面记录延拓为地下波场的过程，在此 过程中，绕射波得到收敛，倾斜界面反射波得到归位，波 场干涉得到分解，波前回转现象得到消除，界面折射得以 校正(深度偏移)，从而使地层构造、断层分布、断点、尖灭 点、边缘、异常体和岩性变化得到清晰成像和准确归位。
DFT
:
F(k) f (n)
N 1
f (n)e j(2 / N )k n
n 0
1
N 1
F (n)e j(2 / N )k
n
N n0
振幅谱与相位谱也可以写成离散的形式。
F(s)
f (t )est dt
LT :
f (t)
（1）振动图——固定空间位置，观察r处质点位移随时间变 化规律的图形。（一道地震记录=一个振动图）
周期T： 一次全振动所需要的时间；
频率f=1/T: 单位时间内全振动的次数。
地震波不是简谐波，从振动图中 可得到相邻两峰或谷间的时间称 为视周期T*，其倒数为视频率 f*。
（2）波剖面——固定某时刻，观察质点位移随距离变化规律 的图形。
时间延迟也可表示为相位延迟（时间延迟/时间周期）
时 间 (s)
频率(Hz)
图1.1-4 图1.1-2的部分放大图，以便更好地从一个频率到另一个频率勾划相位曲线的趋势。用正 峰P所显示的趋势与图1.1－3中的相位谱比较
相位定义为负的相位 延迟。这样，一个负 时延相当于一个正相 位值。
图I-12 经六家处理公司处理的同一条地震测线。（数据由British Petroleum Development. Ltd: Carles Exploration Ltd; Clyde Petroleum Plc;Goal Petroleum Plc; Premier Consolidated Oilfields Plc; and Tricentrol Oil Corporation Ltd提供）
T
频域抽样,即 (2 / N )k
频域恢复
F (k )
DTFT
ZT
:
F(z) f (n)
n-
1
f
(n)z n F ( z)z
n1dz
j2 c
DFT
:
F(k) f (n)
N 1
f (n)e j(2 / N )k
S (ln z) / T
ZT
z e j
F(z)
IZT
e j z
IDTFT
F ( j)
截
取 主
频 域 周
周
期
期
延 拓
F (e j )
DTFT
:
F
(e
j
)
f
( n)
f (n)e j n
n-
1 F (e j )e j nd 2
在计算机中用快速算法实现付里叶变换（FFT）。 付里叶变换：
正变换：时域信号 分解 频域信号； 逆变换：频域信号 合成 时域信号。
Hale Waihona Puke 一维付里叶变换一维付里叶变换
f(t)的振幅谱表示f（t）频率为时简谐成份的振幅值，相位谱表示频率为时 简谐成份在t=0时的初始相位。
一维付里叶变换
一维付里叶变换
离散的付氏变换：
从简谐波的波剖面中可以得到:
波长： 传播一个波的距离
波数: k
—单位距离内传播的波的个数。
地震波不是简谐波，从波剖面中可得到相邻两峰或谷 间的距离称为视波长，其倒数为视波数。
地 震 波 场
地 震 波 场 时 间 切 片， 即 波 动 图
一维付里叶变换
一个正弦运动要用频率、振幅和相位才能完整 的描述。
偏移方法：
具有算法多、类别多、变种多、杂交多和应用多的特点。 一般而言，有二维/三维、叠前/叠后、时间/深度偏移。
成像中的偏移
第一章 地震数据处理基础
本章内容：
第一节 一维付里叶变换及其应用 第二节 二维付里叶变换及其应用 第三节 基本（地质-地球物理）模型及地震数据处理特点
波是质点的振动在介质中的传播。故波动是时间t和空间位 置r的函数，即 u(t, r) 。