地震处理成像系统新技术

Seismic Processing and Imaging
地震处理成像系统新技术

一体化叠前处理成像系统
IMAGING
GeoDepth
Velocity Model Building & Depth Imaging
Focus
Seismic Data Processing
3D Brick
Interpretation With Modeling

Focus 最新技术进展
一、实现了完整的叠前时间偏移处理流程二、完善了VSP处理功能，实现零偏和非零偏VSP处理三、实现了交互式工区管理，使Focus成为完全交互处理系统四、完善了多次波处理的SRME反演方法和PRT方法五、增加了时频谱分解，提供了薄层检测、断层研究以及能量衰减研究的新方法；六、增加了叠前/后去噪和二维滤波模块，去噪手段得到加强；七、增加了LIFT智能去噪技术，增强了Focus处理低信噪比等复杂资料的能力八、采用Epos3.0 体系结构，支持PNS管理方式加强了库函数和系统的稳定性九、采用新的并行处理机制，增强了中断处理和负载均衡功能

Focus 最新技术进展
? ? ? TFDCOMP TFCLEAN VSPARC
新增地球物理模块
Attributes with time-frequency decomposition Noise Suppression in t-f domain
First Arrival Offset Correction Factor for near offset ? VSPREF Refracted Ray Path Correction based on check shot or near VSP data ? VSPMAP VSP to CDP Transform using iterative refining ray trace scheme ? VSPCDPS Elasticized Trace CDP Stacking

Focus 最新技术进展
? LIFT
LFAF AVO
新增地球物理模块
Leading intelligence Filter Technology (Low Frequency Array Filters ) (reflectivities from band-limited seismic
data )
AMPSCAL (Noise Attenuation) HDRMAN (Signal Reconstruction )
? Suppres band-limited noise suppression

Focus 最新技术进展
? ? ? ? ? ? ? ? ? ? ? EXTRACT STATICR UNIFORM REGLO QCSTATS FILT2D CORINT SMACTRM VELDEF CVELDEF ETADEF
修改地球物理模块
Database Geometry Definition from Data 3D Residual Statics Solution Offset or signed offset regularization Prestack trace interpolation Seismic trace statistics 2D user-defined filter application 3 Component vertical orientation correction Surface related multiple attenuation Velocity Analysis Vp/Vs ratio analysis Eta analysis

一体化处理系统的一些新技术
1、时间域处理流程的优化和模块功能增强 LIFT叠前去噪带限噪音压制 SRME 多次波压制高效并行
2、精确的速度建模方法 - 丰富的速度修正手段 - 高精度层析成像技术 - 完善的实体模型建立手段 - 各项异性参数分析 3、适应条件更广，更精确、更快的偏移算法：考虑VTI、TTI条件下的 PSTM 、 PSDM 新的PSPCPI –Phase 算法比以前速度快3-5倍 3D 炮域波动方程叠前时间偏移转换波叠前偏移 CRAM共反射角叠前偏移

New Geophysical Applications
Noise Suppression – LIFT?
Signal Model
Mostly signal
Input
AVO Reconstruction FX Prediction or Projection Low Frequency Array Forming
Output
Suppressed Noise
T-F Decomp
Higher N to S ratio
Noise

Signal Reconstruction and Adaptive Subtraction (LIFT) Data Examples

3D Inline Without LIFT
Refraccion Noise
Refraccion Multiple
AirBlast

3D Inline With LIFT

Band-Limited Noise Suppression
Band-Limited Noise Suppression SUPPRES performs noise suppression for band-limited noise types by band-pass filtering and envelope equalization.
Application to Ground Roll (Above) And Swell Noise (Below)

地震资料解释

地震资料解释期末复习（王松版） 1地震资料解释——以地质理论和规律为指导，运用地震波传播理论和地震勘探方法原理，综合地质、测井、钻井和其它物探资料，对地震数据进行深入研究、综合分析的过程。 2地震子波（wavelet）：地震勘探过程中，爆炸产生的尖脉冲传播到一定距离时波形逐渐稳定。 3褶积模型的应用：已知r(t)和w(t)，求s(t)：正演问题已知w(t) 和s(t) ，求r(t) ：反演问题已知s(t) 和r(t)，求w(t)：子波处理 4同相轴：指地震时间剖面上相同相位的连接线 5极性判断 6有效波的识别标志 1）强振幅：叠后资料往往经提高信噪处理，反射波能量大于干扰波能量 2）波形相似性：子波相同、同一界面反射波传播路径相近，传播过程影响因素相近，相邻地震道上的波形特征（主周期、相位数、振幅包络形状等）是相似的。 3）同相性：同一个反射波的相同相位，在相邻地震道上的到达时间也是相近的，每道记录下来的振动图是相似的，形成一条平滑的、有一定长度的同相轴，也称相干性。 4）时差变化规律：在共炮点道集上，直达波、折射波是直线，反射波、绕射波、多次波等为曲线。在动校正后的剖面上，原来直线的同相轴被校正成曲线，一次反射波成为直线，多次波、绕射波为曲线。 1、2用于识别波的出现； 3、4用于识别波的类型、特征及地层界面特征的判断。 7水平叠加剖面的特点（1）在测线上同一点,根据钻井资料得到的地质剖面上的地层分界面,与时间剖面上的反射波同相轴在数量上、出现位置上,常常不是一一对应的。（2）时间剖面的纵坐标是双程旅行时t0 ，而地质剖面或测井资料是以铅垂深度表示的,两者需经时深转换,其媒介就是地震波的传播速度,它通常随深度或空间而变化。（3）反射波振幅、同相轴及波形本身包含了地下地层的构造和岩性信息,如振幅的强弱与地层结构、介质参数密切相关。但是反射波同相轴是与地下的分界面相对应,同相轴与界面两侧的地层、岩性有关。必须经过一些特殊处理（如声阻抗反演技术等）才能把反射波所包含的“界面”的信息转换成为与“层”有关的信息后,才能与地质和钻井资料进行直接地对比。（4）地震剖面上的反射波是由多个地层分界面上振幅有大有小、极性有正有负、

地震勘探复习资料

绪论 1、地球物理勘探的概念（1）简称“物探”，是通过观察存在地球及其周围的地球物理场的特征和岩石的各种物理特性来研究地质规律和勘查各种矿产的各种方法的总称。（2）是以物理学原理为基础，利用电子学、计算机的数字处理、信息论等科学技术中的新技术所建立起来的一整套勘探地下矿产的方法。（3）是借助于各种物探仪器在地面观测地下岩石的各种物理参数，从而解释和推断地下岩石的构造特点、岩石性质等，从而到达勘查地下矿产（金属非金属矿产、煤、油气等等)的目的。 2、地球物理勘探的分类，不同勘探方法的优缺点。重力勘探：利用岩石的密度差异磁法勘探：利用岩石的磁性差异电法勘探：利用岩石的电性差异地震勘探：利用岩石的弹性差异放射性勘探：利用岩石的放射性差异地震勘探的优点：精度高，分辨率高，穿透深度大，能较详细地了解由浅至深一整套地层的地质规律。缺点：成本高 3、地震勘探的概念、分类，目前地震勘探以何种方法为主。概念：利用岩石的弹性差异来进行矿产勘察。是通过人工激发地震波，研究地震波在弹性不同的地下地层中传播的规律，以查明地下的地质构造，达到油气或其他勘探目的的一种物探方法。分类：地质法（优：在找油初期，可以起到一个指向作用，避免了盲目性，成本低。缺：野外地质方法很难准确了解地下地质情况！）；钻探法（优点：精度最高，缺点：一孔之见，而采用大量的钻井，不仅成本高，而且效率低）；物探方法（优点：精度高于地质法，成本低于钻探法；不足：精度低于钻探法，成本高于地质法）。应用最多的方法：物探方法 4、地震勘探的三个阶段地震资料野外采集、地震资料室内处理、地震资料解释。第一章各种介质的概念重点：①物体是否为弹性、塑性介质与受力大小、时间及温度有关。②均匀介质与各向同性介质的关系。（1）理想弹性介质：当介质受外力后立即发生形变，而外力消失后能立即完全恢复为原状的介质；（2）粘弹性介质：当外力消失后不是立即恢复原状，而是过一段时间后才恢复原状的介质称为粘弹性介质。（3）塑性介质：当外力消失后不能完全恢复原状,保留了一部分形变的介质称为塑性介质。（4）各向同性介质：凡介质的弹性性质与空间方向无关的介质称为各向同性介质（5）各向异性介质：凡介质的弹性性质与空间方向有关的介质称为各向同性介质（6）均匀介质：弹性性质（波速）不随空间坐标的变化而变化，是常数。（7）非均匀介质：弹性性质（波速）随着空间坐标的变化而变化，不是定值。（8）层状介质：如果非均匀介质的物理性质呈层状分布，则称这种介质为层状介质。层状介质中各层的弹性系数是不变的。层状介质模型已经成为地震勘探中常用的物理模型。（9）连续介质：层状介质的层数无限增加，每层的厚度无限减小时，层状介质就可以视

地震数据处理方法(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，是一种几何原因造成的某处能量的减小，与介质无关，叫几何扩散，又叫球面扩散。为了消除球面扩散的影响，只需A0=Ar*r即可，此即为几何扩散校正， 13、反滤波(又称反褶积):为了从与干扰混杂的地震讯息中把有效波提取出来，则必须设法消除由于水层、地层等所形成的滤波作用，按照这种思路所提出的消除干扰的办法称为反滤波，即把有效波在传播过程中所经受的种种我们不希望的滤波作用消除掉。 14、校正不足或欠校正：如果动校正采用的速度高于正确速度，计算得到的动校正量偏小，动校正后的同相轴下拉。反之称为校正过量或过校正。 15、动校正:消除由于接受点偏离炮点所引起的时差的过程，又叫正常时差校正。 16、剩余时差：当采用一次波的正常时差公式进行动校正之后，除了一次反射波之外，其他类型的波仍存在一定量的时差，我们将这种进过动校正后残留的时差叫做剩余时差。

地震解释技术

随着锦州油田油气勘探开发的不断深入，先进的三维地震解释技术及相关的属性分析技术的使用凸显重要。利用最新采集处理的三维地震资料,采油厂加大了相关地震配套软件的使用，2011年锦州采油厂计划引进SeisWare地震解释系统及landmark地震解释工作站，使得利用各种地震属性研究储层的技术得到了加强。利用高精度三维地震叠前时间偏移数据体，可以在精细地层小层对比、整体解剖精细评价的基础上针对目标层段内的砂泥岩薄互层砂组进行多种地震属性的处理，引进landmark解释工作站的多体多属性地层追踪及快速高效的储层描述方法，能从整体上描述储层的空间展布及小断块内储层的分布特征，计算机技术的飞速发展及相应的层位自动追踪技术、三维可视化技术等解释手段的发展极大地提高了解释工作的效率及准确度,同时最大限度地发挥了三维数据体的优势。利用最新采集处理的三维地震资料,经过地震资料品质分析后,优选具有较高的信噪比,偏移归位合理,目的层波组特征明显的资料,在合成记录标定的基础上,搭建格架剖面并进行人工解释,然后采用人机联合波形对比层位自动追踪技术进行全区层位解释,采用相干、倾角扫描以及层面光滑度分析技术进行断层平面组合分析,能精细落实研究区的构造特征和断层展布特征。

LandMark 一体化系统通过强有力的可视化技术提供给用户一个真三维的解释平台，可对海量的三维地震数据进行快速准确地构造解释，能快速搜索地质目标，精确雕刻；并提供了一个多学科协同和决策环境，可以实现构造解释、储层预测、叠前AVO分析、可视化处理以及井轨迹设计和钻井实时监控。其三维可视化手段可应用于地震资料处理、构造解释、全区目标搜索、精细目标解释、储层预测等三维连片解释的所有阶段。 LandMark 一体化系统特点：储层自动追踪ezTracker 基于波形的层位自动追踪，可同时拾取多个种子点，可以保存种子点信息，灵活定义追踪的波形时窗，对追踪结果可进行多种灵活编辑，如遗传删除、门槛值调整和多边形删除点集自动追踪Autopick 可根据种子点值的大小，或人工定义数据体值的范围，快速追踪地质体。也可利用多种属性（如在波阻抗体和相位体上）共同约束追踪地质体三维形态，如河道、扇体等，直接形成地质体顶底t0面。点集可自由转换为层位。三维体雕刻Geobody 可用三维体追踪点集,层位,断面作为约束条件雕刻三维地质体，利用透明度和颜色来彰显地质异常体，突出空间展布。异常体快速搜索GeoAnomaly 依据多数据体振幅值和数据连通性，快速搜索满足定义条件的异常体。 SeisWare软件的地震地质解释功能灵活方便，适于在勘探/开发阶段进行综合地震解释、随钻跟踪分析、油气层识别、储量计算以及新区预探、老区扩边、部署调整等研究工作。其特点包括：多工区，不同类型地震资料的连片解释；断层追踪识别功能可以直观方便的显示地震剖面上断层的平面要素，实时地观察断层面的空间走向及展布趋势。欢西油田是一个地质条件和油藏来信十分复杂的断块油田，断距从十几米至几百米不等的不同级次断层纵横交错，断块分隔凌乱，油层埋藏差异大，储层沉积特征不一，发育不稳定，诸多因素都给地质研究带来困难。面对复杂断块，Seisware地震解释系统的技术优势是，可以直观方便地显示地震剖面上断层的平面要素，实时地观察断层面空间走向及展布趋势，并使三维数据断层解释过程自动化。地震解释人员可以能够在较短时间内进行高精度的断层解释，即使在构造情况复杂地区或资料品质较差地区也能实现，其直观的编辑功

地震资料处理解释大作业(处理部分)

地震资料处理/解释大作业（处理部分）专业：勘查技术与工程班级：12-4 姓名：封辉、孙运庆、何瑞川学号：2012011236、2012011249、2012011239 2016年 1 月 15 日评分标准：第三章和第四章各20分，其余各章10分

目录第一章数据加载和观测系统定义 (2) 第二章道编辑和真振幅恢复 (4) 第三章反褶积 (6) 第四章速度分析 (7) 第五章动校正和水平叠加 (8) 第六章静校正 (10) 第七章偏移 (12) 第八章总结和体会 (13)

第一章数据加载和观测系统定义地震资料处理流程第一步为数据输入和预处理。预处理是地震数据处理前的准备工作，将地震数据正确加载到地震资料处理系统，进行观测系统定义，并对数据进行编辑和校正。原数据是SGY格式的地震记录文件，用Promax对其进行处理需要格式转换，将其格式转换成软件定义的格式。图1.1是原始数据炮集。格式转换后可对数据进行加载与处理，但是处理需要的各种测网信息需要进行定义，所以我们做观测系统定义，用FFID（野外文件号）和CHAN（记录道号）为索引将测网的各检波器与炮点坐标、高程、CDP 号等信息与数据的各道联系起来。观测系统定义分为炮点定义，检波点定义与炮检关系定义。图1.3是CDP覆盖次数。图1.1 原始数据炮集

图1.2a 炮点与检波点信息图1.2b 炮点与检波点信息

图1.3 多次覆盖次数第二章道编辑和真振幅恢复通常的地震采集中，由于检波器数量很多、野外干扰因素复杂等原因，不是每一道都能很好的反应地下反射界面带回来的信息，最基础的我们需要挑出其中坏检波器采集的道与极性不正常的道，称为道编辑（如图2.1）。在记录图中使用picking进行编辑。点击picking，有编辑错道和编辑极性翻转道。拾取所有的错道和翻转道集后，分别放在两个文件里面。由震源引发的地震波，会随着波前面变大，底层吸收衰减等因素而能量减小，而我们需要的通常是深部的地层信息，所以我们需要对地震波进行振幅恢复（如图 2.2），经过真振幅恢复以后，深层反射波能量相对增强了,反射界面变得清晰，但面波等干扰波也增强了。

石油地震勘探资料处理

石油地震勘探资料处理 1.地震资料数字处理是怎么回事？既然野外地震已经采集到了反映地下地质情况的地震记录，为什么还要进行地震资料数字处理呢？这是因为野外采集的地震记录仅仅是把来自地下地层的各种信息以数码形式记录在磁带上或光盘上，还不能直接反映出地下地层的埋藏深度及起伏变化情况，还需要将地震记录拿到室内输入到运算速度非常快、存贮量非常大、专业功能非常强的计算机系统中，在专家的指令下进行反复计算和分析，才能获得直接反映地下地层真实情况的数据和图像，专业上把这一过程叫做地震资料数字处理。这个过程有点像我们生活中使用的数码照相机（或数码摄像机）的显像过程，将数码照相机拍摄到的图像输入到室内的电脑上，根据需要，对显示在屏幕上的影像进行修改、调整、增加、删减，满意后可通过屏幕拷贝、彩色打印输出图片来，也可以录制到光盘上存贮以供调用，这个过程叫做编辑，也叫处理。不过地震资料的数字处理所用的硬、软件则要复杂得多。因为数码相机拍摄到的图像仅是几米到几十米远的景物，而地震资料数字处理要对从地面开始到地下五六千米甚至上万米深范围内的地震数据进行处理，不仅将上面第一套地层，还要将下面很多套地层逐层搞清楚。这些地层在不同地区形态都不一样，有的很平，有的像喜马拉雅山似的高山，有的像雅鲁藏布江似的河谷。可见地震数字处理要把地下数千米深的看不见、摸不着，又极其复杂的地层情况搞清楚，这是多么难的一门学科。不过，近些年来由于将迅速发展起来的计算机技术、信息技术等许多高新科学技术引用到地震资料数字处理中，为搞清地下地层情况，寻找深埋地下的油气田提供了条件，提供了可能，而且提高了油气勘探的成功率。经过数字处理后的成果有好几十种。专业上把反映地层的埋藏深度、厚度以及形态的图件叫做水平叠加剖面（简称叠加剖面）、偏移剖面。把反映地层岩石（砂岩、泥岩等）组成及其物理性质（速度高低、孔隙大小等）等的成果叫地震属性资料。将经过数字处理的这些剖面和属性资料录制到数字磁带或光盘上，可提供给下道工序（解释）使用。

地震数据处理vista软件使用手册

Vista 5.5的基本使用方法数据输入地震分析窗口一维频谱二维频波谱观测系统工作流一、数据输入 1.1 把数据文件加入Project 首先选择File/New Project，新建一个Project，按住不放，出现按钮组合，可以选择不同类型的数据集，选择，向Project中增加一个新的2-D数据集，按住不放，出现按钮组合，可以选择加入不同类型的地震数据，选择，选择一个SEG-Y数据，即可将该数据文件加入新建的数据集。 1.2 命令流中数据的输入双击进入如下界面 1.2.1 Input Data List 数据输入列表，选择已加入到Project的数据集，下面的文本框中会显示选择的数据的基本信息。 1.2.2 Data Order 选择输入数据的排列方式，对不同的处理步骤可以选择不同的数据排列方式 Sort Order a. NO SORT ORDER 输入数据原始排列方式 b. SHOT_POINT_NO 输入数据按炮点排列方式 c. FIELD_STATION_NUMBER d. CMP_NO 输入数据按共中心点排列方式 e. FIELD_STATION_NUMBER 1.2.3 Data Input Control 数据输入控制右键-->Data Input Control a. Data Input 进入Flow Input Command（见上） b. Data Sort List 查看数据排列方式的种类 c. Data/header Selection 输入数据的选择，可以控制输入数据的道数和CMP道集查看所有已经选择的数据如果没有定义任何可选的数据信息，则如下图所示：可以选择一种选择方式，单击并设置选择信息。定义有可选的数据信息后，在查看，则如下图所示，会显示选择的信息。选择共炮点集单击后，会弹出如下界面：

地震数据处理重点

1.一维傅里叶变换及其应用：傅里叶变换是地震数据处理的主要数学基础。它不仅是地震道、地震记录分析和数据滤波的基础，同时在地震数据处理的各个方面都有着广泛的应用。 2.采样定理：设x(t)是连续的时间函数，x(t)的最高截止频率为fn，则可用采样间隔为Δt=1/2fn的离散序列X(nΔt)唯一的确定。采样过程：从模拟地震信号到数字地震信号的过程。采样间隔/采样率：采样所用的时间间隔。 3.数字滤波：利用频谱特征的不同来压制干扰波，以突出有效波的方法。 4.频率域滤波的步骤： ①对已知地震道进行频谱分析；②设计合适的滤波器：为了滤去干扰波的频谱成分，应当设计一个带通滤波器，保留有效波频率，把干扰波频率成分滤掉； ③进行滤波运算；④对输出信号谱X(w)进行傅里叶反变换，便得到滤波后的输出X(t). 5.相位性质：最小相位也叫相位滞后或最小能量延迟，实际上最小相位滞后是指频率域，而最小能量延迟则是指时间域而言。最小能量延迟子波：能量聚集在首部；最大能量延迟子波：能量集中在尾部；混合延迟子波：能量聚集在中部。 6.褶积滤波的物理意义：单位脉冲响应：在时间域的表示方法中，令一个单位脉冲通过一个滤波器，然后观测滤波器的输出，这个滤波器输出的自然过程曲线称为滤波器的脉冲响应。也称滤波器的时间特性。褶积滤波的物理意义：它相当于把地震信息x(t)分解为起始时间、极性、幅度各不相同的脉冲序列，令这些脉冲按时间书序依次通过滤波器，这样在滤波器的输出端就得到对输入脉冲序列的脉冲响应，这些脉冲响应有不同的的起始时间、不同的极性和不同的幅度（这个幅度是与引起它的输入脉冲幅度成正比的），将它们叠加起来就得到滤波后的输出x(t). 7.数字滤波的特殊性质：离散性：数字滤波是对离散的信号进行运算，这是所谓的离散性；有限性：在数字计算机上进行计算时，滤波因子不可能无穷项，而是取有限项，这就是所谓的有限性。 8.产生“伪门”原因：由于对A离散采样造成的，可以证明“伪门”在频率域出现的周期为A，为了避免“伪门”造成的影响，可以适当的选择采样间隔A，使第一个“伪门”出现在干扰波的频谱范围之外。9.波谱：以任何一种形式展示电磁辐射强度与波长之间的关系，叫波谱。波数：波长的倒数。K0=1/λ 二维频率-波数域中的二维频率-波数谱（简称二维频-波谱）分析是对地震波场进行分析的重要手段，它是建立在二维傅里叶变换的基础上。 10.空间假频：频率不变，倾角越大或者倾角不变，频率越高越容易产生空间假频。产生条件：地震信号的频率f一定时，地震信号倾斜时差δt越大，其频-波振幅谱中的波数k0也越大，而当地震信号频率f 增大时，具有相同倾斜时差δt的地震信号的频-波振幅谱中的波数k0随之增大，当频率f增大到某一个门槛频率fmax时，便开始产生空间假频。 11.二维滤波器的设计：一般二维滤波是指对于波动函数X(t,x)所进行的频率-波数域滤波。这时设计的滤波因子是时间-空间的函数h(t,x)，滤波过程类似一维滤波在时间-空间域，可用二维褶积公式表示A. 12.共中心点CMP叠加及叠后处理流程图：野外采集地震数据-解编-预处理-反褶积-抽CMP道集-速度分析-动校正-CMP水平叠加-叠后时间深度偏移。13.共中心点叠加优点：①压制多次波；②压制规则干扰波；③压制随机噪声。综上，共中心点叠加可以有效地压制各种干扰波，增强有效波，使地震剖面的信噪比明显提高，掀桌改善地震剖面的质量。 14共中心点水平叠加存在的问题：当反射界面为弯曲界面时，其反射旅行时存在如图1所示的畸变；当反射界面为，其射旅行时发生如图2所示的畸变；当覆盖介质速度横向变化时，其反射旅行时存在如图3所示的畸变；当覆盖介质速度各向异性时，其反射旅行时存在如图4所示的畸变. 15.块状介质模型地震数据处理的特点：①介质呈块状分布，它不仅有顶部和底部界面，而且其侧面也由断层面或岩层界面所封闭；②由于剧烈的构造运动作用，界面往往呈弯曲界面，界面陡、倾角较大；③介质速度往往沿水平方向变化较快。 16.共反射点CRP叠前处理基本流程图：野外采集地震数据-解编-预处理-反褶积-抽CRP道集-层速度场-速度深度模型-叠前深度偏移 ①②③④⑤⑥⑦ 1.预处理：指地震数据处理前的准备工作，是地震数据处理中的重要基础工作，一般定义为将野外采集的地震数据正确加载到地震资料处理系统，进行观测系统定义并对地震数据进行编辑和校正的过程。预处理包括：数据解编、格式转换、道编辑、观测系统定义等工作。 2.解编：就是按照野外采集的记录格式将地震数据检测出来，并将时序的野外数据转换为道序数据，然后按照道和炮的顺序将地震记录存放起来。 3.野外观测系统定义：观测系统就是以野外文件号和

地震勘探资料处理

本科生实验报告实验课程基于 Vista 系统的地震资料处理学院名称地球物理学院专业名称勘查技术与工程（石油物探）学生姓名学生学号指导教师唐湘蓉实验地点5417 实验成绩 2015年3月- 2015年5月

基于 Vista 系统的地震资料处理一、实验目的及要求 1）认知熟悉地震资料处理软件系统--vista软件的基本功能，了解其并熟练掌握vista软件运行的基本操作； 2）了解并掌握地震数据处理的基本流程，掌握地震数据处理的流程和基本方法，选择合适的处理参数以提高地震数据处理的精度； 3）对比地震资料处理与解释的理论与实际资料处理的结果，深入理解理论，并在理论指导下提高处理解释的水平、提高资料处理的质量； 4）提高综合分析问题的能力与编写实验报告或生产报告的能力。二、实验内容总流程图1 总流程图 1）加载数据打开Vista软件后选择加入2D的SEG-Y格式的原始地震数据，本实验

所用数据为给定的SHOT-20。加载后的原始地震数据如图2：图2 原始地震数据显示 2）道均衡各个道由于炮检距的不同，导致的反射波的振幅的变化，因为在共反射点叠加中，要求每一个叠加道的振幅都应该相等，每一道对叠加所做的贡献是等价的，无特殊情况，一般就以记录图中间的振幅为基准，使近激发点的地震道振幅减少，增加远离激发点的地震道记录的振幅。道均衡流程模块如图3，道均衡结果如图4：图3 道均衡流程模块

3）建立观测系统图5 观测系统显示4）初至拾取初至拾取结果显示如图6：

图6 初至拾取结果显示 5）初至切除地震记录上的初至波包括直达波和浅层折射波，它们能量强且有一定延续时间，对紧接而来的浅层反射波有干涉和破坏作用。另外，动校正后会引起波形畸变，浅层尤其厉害。对这些强能量初至波和动校正畸变引起的处理办法是“切除”，即将这些波的采样值全部变为零值（充零）。初至切除流程模块如图7，初至切除结果如图8：图7 初至切除流程模块

地震勘探资料数字处理

中国地质大学（北京）课程名称：应用地震学教师：段云卿第25册第四章：地震勘探资料数字处理野外采集到的原始资料是以二进制的数字形式记录在磁带上，必须经过计算机的各种运算，才能输出供地震地质解释的各种资料，或直接输出某些解释成果，本章介绍如何进行数据处理。 §4.1校正和叠加处理一、动校正 1．动校正的含义：（§3.5）（1）对于一次覆盖共炮点资料来说，把双曲线型或近似双曲线型反射波同相轴拉直，也就是消去炮检距不为0对反射波旅行时的影响，使同相轴能直观地反映地下界面的构造形态。（2）对于共反射点道集来说，把各道均校正成共中心点M 处的自激自收道，再叠加起来作为共中心点M 处的叠加道，使一次波同相叠加而加强，多次波等干扰波非同相叠加而减弱。 2．动校正公式（§3.5） 2 022V t x t = ? (6.2-26) 3．计算动校正量（使用共反射点道集） (1)公式为了对共反射道集的每一道的整个道进行计算，将（6.2—26）改写为： 2 002 ) (2i i j ij t V t x t = ? （j=1,2,……,n ； i=1,2,……,m ） (6.4-1) j —— 道序号。 i —— 采样点序号。 x j —— 第j 道的炮检距。 n —— 覆盖次数。 M ——道长 t 0i ——为第i 个界面共中心点处自激自收时间。 (2)问题不知什么地方有反射界面，就不知什么地方有反射波。不知反射波的t 0时间。

中国地质大学（北京）课程名称：应用地震学教师：段云卿第25册（3）解决方法地震道上有一个采样值就有一个反射波。地震道上每一个采样点的时间i △，都看成一个t 0时间，记为t oi 。（4）例子 ①设采样间隔△=4ms ②长为0.5S -4.5S 的记录,就有1001个t 0值： )(5.00,0s t = )(004.05.01,0s t += )(004.025.02,0s t ?+= )(004.05.0,0s i t i += )(004.010005.01000,0s t ?+= ③对任意一道就有1001个动校正量。例如炮检距为1000m 的第j 道，动校正量为： )(207.0) 5.0(5.021000 2 2,0s V t j =??= ? ) (205.0) 504.0()504.0(21000 2 2,1s V t j =??= ? ) (204.0) 508.0()508.0(21000 2 2,2s V t j =??= ? ) ?() 004.05.0()004.05.0(21000 2 2,s i V i t j i =+?+?= ? )(000.0) 5.4()5.4(21000 2 2,1000s V t j =??= ?

地震资料解释报告材料

地震资料解释报告序言勘查技术与工程卓越班的实践性很强，加强实践教学可以提高学生的动手能力和处理实际问题、分析解决实际问题的能力、使之能更好的适应毕业后实际工作，是一个非常重要的教学环节，也是进一步提高教学质量的重要途径之一。我们的地震资料解释实践共分两步完成，第一是在学校手工地震资料构造解释课程设计，第二是在东营对news软件的学习。此次实习是在完成了《地震勘探原理》和《地震资料解释》的基础上完成的实习，通过此次实习的机会我们得以理论联系实际并用实践以检验所学理论，各项安排有条不紊的展开。在每一步的实习过程中都有老师的带领，手工地震资料构造解释课程设计由杨国权老师负责，news软件的学习由张繁昌老师负责。实习过程中注意理论和实际的结合，在老师的带领及同学的相互帮助下，我们顺利的完成了实践所要求的所有内容。

目录一、实习目的及意义 (4) 二、实习内容 (4) 三、地震资料构造解释 (5) 四、News学习 (7) 五、结论与建议 (26)

一、实习目的及意义通过课程的学习，对解释软件系统、数据的地质地球物理解释过程等有基本的认识和掌握，通过实习熟悉了勘探方法的整个工作原理和处理解释流程以及实习报告编写等过程。了解到了反射波的追踪对比、地震资料的地质解释、构造图的绘制、以及研究成果的提交等过程。培养实际技能及对分析和解决实习问题的能力；掌握仪器的工作原理，并学会操作和使用；掌握各方法的基本数据分析和处理技能。对本专业所从事工作的性质、手段、方法以及新技术、新方法有一个全面的了解，培养学生的实际操作和计算技能以及综合分析问题的独立工作能力，巩固已学过的专业知识，为下一步进入专业课程和毕业论文阶段以及今后走上本专业的工作岗位打下基础。二、实习内容地震自资料的构造解释内容主要有工区的地质情况总结、地震资料解释流程、对地震构造解释的分析、体会和建议等。News 的实习内容主要在理论学习好的基础上，学会利用软件完成地震资料解释的整个过程，并得出理论成果。三、地震资料的构造解释构造解释是以水平叠加时间剖面为主要资料，利用由地震资料提供的反射波旅行时间、速度等信息，查明地下地层的构造形态、埋藏深度、解除关系等，通过构造解释成果，即使提供钻井井位。构造解释的三大环节：

成都理工大学地震勘探资料处理及解释复习资料及答案

1----断层在时间剖面的特征标志 1）标准层反射同相轴发生错断，是断层在地震剖面上表现的基本形式。2）标准层反射波同相轴数目突然增减或消失，波组间隔发生突变，断层下降盘地层加厚，上升盘地层变薄。3）反射同相轴形状和产状发生突变，这往往是断层作用所致。4）标准层反射波同相轴发生分叉、合并、扭曲及强相位转换等。5）断面波、绕射波等异常波的出现，是识别断层的主要标准。 2----伪门条件及消除方法滤波处理的是离散信号，由付氏变换的特性可知：离散函数的频谱是一个周期函数，其周期为1/△，即有：DFT(h(n))=H(k)=H(k+1/Δ)则通频带以1/△为周期重复出现，若称第一个门为“正门”，则其它的门为“伪门”。②克服的方法：a)选择适当的采样间隔△使伪门出现在干扰波频率范围之外，一般采样间隔△取得越小，伪门处于频率越高的地方，离正门越远，b) 在离散采样之前让信号通过“去假频”滤波器，滤掉高频成分。 3--反滤波原理及影响因素地震记录是地层反射系数序列r(t)与地震子波b(t)的褶积，x(t)=r(t)*b(t)，b(t)就相当地层滤波因子。为提高分辨率，可设计一个反滤波器，设反滤波因子为a(t)，并要求a(t)与b(t)满足a(t)* b(t)=(t)，用a(t)对地震记录x(t)反滤波x(t)* a(t)= r(t)*b(t) * a(t)= r(t)* (t)= r(t),其结果为反射系数序列，即为反射波的基本原理。影响因素：1）各种反滤波方法都必须有若干假设条件；2）反射地震记录的褶积模型问题；3）噪声干扰的影响；4）原始地震资料的质量问题。4----．爆炸反射界面成像原理（叠后偏移成像原理）①把地下地质界面看成具有爆炸性的爆炸源。②爆炸源的形状、位置与地质界面一致。③爆炸源产生的波的能量、极性与地质界面反射系的大小、正负对应。④并假定当t=0时，所有爆炸源同时起爆，沿界面法线方向发射上行波到达地面观测点。（5）用波动方程式将地表接收的波场（地震记录）作反时间方向传播（向下延拓），当波场延拓到（t=0)时的波场的值就正确地描述了地下反射界面位置，即自动实现偏移成像。说明：爆炸反射界面成像原理适用于叠后的地震资料。即自激自收剖面，自炮点发出的下行波到达反射点的路径与自该点反射返回地面的上行波的路径完全一样。只考虑上行波，若将时间剖面中时间减半，或将传播速度减一半，就可将自激自收剖面看作在反射界上同时激发的地震波沿界面法线传播到地表所接收的记录。偏移时，只需把速度减半，用单程波动方程延拓法，把波场从地面延拓到反射界面，令t=0，即可实现偏移。 5.有限差分法波动方程偏移有什么特点 ①是求解近似波动方程的一种近似数值解法，是否收敛于真解，取决于差分网格的划分和延拓步长的选择。②能适应速度的纵、横向变化，偏移噪音小，在剖面信噪比低的情况下也能做的优点；③受反射界面倾角的限制，当倾角较大时，产生频散现象，使波形畸变。法波动方程偏移有什么特点偏移结果好，精度高，稳定性好，噪音低，运算速度快，无倾角限制，无频散现象。优点：偏移结果好，精度高，稳定性好，噪音低，运算速度快，无倾角限制，无频散现象。缺点：假定传播速度为常速，速度横向变化时，会使反射界面畸变，对偏移速度误差较敏感。7克希霍夫积分偏移有什么特点与绕射扫描叠加的区别是什么不受倾角限制，能适应任意倾角地层，做三维偏移较容易实现，对网格要求较灵活。优点：不受倾角限制，能适应任意倾角地层，做三维偏移较容易实现，对网格要求较灵活。缺点：费时；难以处理速度的横向变化；偏移噪声大，“划弧”现象严重；确定偏移参数困难。 -区别：A克希霍夫积分偏移考虑了波的振幅值随传播距离和方向不同的影响，保持了波的

现代地震勘探技术作业

中国地质大学（北京）地震属性综述报告名称: 地震属性综述学生姓名：王丹学号：2010120052 所在院（系）：地球物理与信息技术学院

地震属性分类及其地质意义地震勘探是在地表激发人工震源，由震源所引起的震动以地震波的形式向地下传播，并在一定的条件下向上反射传回地表，然后由地表的仪器(检波器)记录反射回来的地震波，从而得到地震记录(也叫地震资料)；之后对地震资料进行相关的处理与解释便可以间接地反映和得到地下相关信息。由于地下介质是地震波传播的载体，所以地下介质的物理性质，如岩性、孔隙度、密度以及流体性质等都会对传播中的地震波产生影响，如地震波的能量、波形、振幅、频率、相位等将在传播过程中发生变化。而这种影响和变化又将在地震记录中保留相应的信息。所以，通过对地震记录(地震资料)的“深加工”或者特殊处理，将会从地震资料中获取更多的有用信息以为地质服务。在早期的油气资源勘探中，地震勘探的目标主要是寻找地下有利的大尺度的构造圈闭，所以只需利用有限的地震资料信息便可达到目的。但是，随着油气勘探与开发难度的加大，人们迫切地需要更多地了解地下地层的岩性、流体性质等信息。这就促使人们运用新的技术和思想去从地震资料中发掘出更多的有用信息。从而，也就推动了地震属性技术的出现与发展。地震属性技术延伸了人类的视觉,从而有助于人们发现更多的隐藏于地震资料中的信息，也有助于人们从多角度去获取和分析地下地质信息，从而实现对地下地质的充分与准确认识。 1地震属性的发展与分类随着油气勘探、开发工作的深入，也为了充分、有效地利用获取不易的地震资料，现今的地震解释人员需要从地震数据中提取越来越多的信息，然后利用这些信息综合解释地下构造、地层和岩性特征以及流体性质，最终定义精确的油藏模型，用于钻井决策、估计地质储量和可采储量。由于生成地震属性是获取所需信急的一条重要捷径，因此，长期以来地震属性技术一直是地震特殊处理和解释的主要研究内容。地震属性是叠前或者叠后地震数据，经数学变换而导出的有关地震波的几何形态、运动学特征、动力学特征和统计学特征。长期以来以来地震数据的使用仅仅局限于对地震波同相轴的拾取，以实现面对油气储集体的几何形态、构造特征的描述。但是地震数据中隐藏着更加丰富的有关岩性、物性及流体成分等相关信

地震数据处理

地震数据整体流程不同软件的地震数据处理方式不同，但是所有软件的处理流程基本是固定不变的，最多也是在处理过程中处理顺序的不同。整体流程如下： 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 静校正（野外静校正）静校正为利用测得的表层参数或利用地震数据计算静校正量，对地震道进行时间校正，以消除地形、风化层等表层因素变化时对地震波旅行时的影响。静校正是实现共中心点叠加的一项最主要的基础工作。直接影响叠加效果，决定叠加剖面的信噪比和垂向分辨率，同时影响叠加速度分析的质量。静校正方法： 1）高程静校正 2）微测井静校正-利用微测井得到的表层厚度、速度信息，计算静校正量 3）初至折射波法 4）微测井（模型法）低频+初至折射波法高频 4 叠前噪音压制干扰波严重影响叠加剖面效果。在叠前对各种干扰进行去除，为后续资料处理打好基础。常见干扰有：面波、折射波、直达波、多次波、50Hz工业电干扰及高能随机干扰等多种情况。不同干扰波有不同特点和产生原因，根据干扰波和一次反射波性质（如频率、相位、视速度等）上的不同，把干扰和有效波分离，从而达到干扰波的去除，提高地震资料叠加效

地震数据处理解释技术发展研究

地震数据处理解释技术发展研究地震数据处理解释是地震勘探的主要组成部分，是石油天然气勘探开发产业链中对油田勘探开发效益影响最大、技术含量最高的一环。…… 一、地震数据处理解释是地震勘探的主要组成部分地震勘探就是通过人工地震反射波“给地球做CT”，让油气勘探者能够“看见”地层的地质构造和油藏情况，为石油公司“找油”做出含油气评价、提出钻井位置、模拟油藏未来的生产动态以便为后续油气藏开采和开发提供技术资料。地震勘探包括地震采集、处理和解释三大部分：地震采集是利用野外地震采集系统获取地震数据处理所需的反射波数据；地震数据处理的目的是对地震采集数据做各种处理提高反射波数据的信噪比、分辨率和保真度以便于解释；地震解释分为构造解释、地层解释，岩性和烃类检测解释及综合解释，目的是利用地震反射波的地质特征和意义确定井位寻找石油。地震数据处理依赖于地震采集数据的质量，处理结果直接影响解释的正确性和精确度和找油的成功率。图1 地震勘探产业链构成地震数据处理解释是地震勘探的主要组成部分，是石油天然气勘探开发产业链中对油田勘探开发效益影响最大、技术含量最高的一环。其原因有四：1、石油勘探地震数据处理解释与井位部署成功率、油田发现、油田采收率、油田增储上产等经济效益直接相关，是寻找油气资源的关键技术； 2、石油勘探技术发展的基础主要体现在地震数据处理环节中地震成像技术的发展；3、地震数据处理解释下游钻井业务等油气开采技术均十分成熟；4、上游地震数据采集依赖于先进的仪器设备，理论简单。综合而言，地震数据处理的质量和地震成像的准确度与清晰度直接决定油气资源的发现的成败和勘探成功率，是影响后期油田生产建设最重要的环节。 BP公司北海油田日产量与地震数据处理解释新技术的关系表明，新技术尤其是地震成像技术的发展和应用对于油田产量的增加影响极大。图2 石油勘探地震数据处理解释技术对北海油田的产量的影响由此可见，地震数据处理解释是地震勘探的主要组成部分，其发展和技术进步对于解决人类能源供应问题具有十分重要的意义。二、地震数据处理解释技术发展历程地震数据处理解释技术中最核心的就是地震成像技术，因此地震数据处理解释技术的发展历程主要依据地震成像技术的发展水平进行划分。地震数据处理解释最早出现于20世纪20年代初期。随后的40年间由于是对光点记录（1920—1950）和模拟记录（1950—1965）进行处理，在这一阶段地震处理解释技术发展缓慢，也没有可实用的地震成像技术出现。

地震数据处理课程设计(报告)

《地震资料数据处理》课程设计总结报告专业班级：姓名：学号：设计时间：指导老师： 2011年5月30日

目录一、设计内容……………………………………………………………… （1）褶积滤波……………………………………………… （2）快变滤波……………………………………………… （3）褶积滤波与快变滤波的比较………………………… （4）设计高通滤波因子…………………………………… （5）频谱分析……………………………………………… （6）分析补零对振幅谱的影响…………………………… （7）线性褶积与循环褶积………………………………… （8）最小平方反滤波……………………………………… （9）零相位转换…………………………………………… （10）最小相位转换………………………………………… （11）静校正………………………………………………… 二、附录………………………………………………………………………… （1）附录1：相关程序…………………………………… （2）附录2：相关图件……………………………………

【附录1：有关程序】 1.褶积滤波 CCCCCCCCCCCCCCCCC 褶积滤波CCCCCCCCCCCCCCCCC PROGRAM MAIN DIMENSION X(100),H1(-50:50),H2(-50:50),Y_LOW(200),Y_BAND(200) PARAMETER (PI=3.141592654) CCCCCCCC H1是低通滤波因子，H2为带通滤波因子CCCCCC REAL X,H1,H2,Y_LOW,Y_BAND REAL dt,F,F1,F2 INTEGER I dt=0.002 F=70.0 F1=10.0 F2=80.0 OPEN(1,FILE='INPUT1.DA T',FORM='FORMATTED',STATUS='UNKNOWN') READ(1,*)(X(I),I=1,100) CCCCCCCCCCCCCCCCCC低通滤波器CCCCCCCCCCCCCCCCC DO 10 I=-50,50 IF (I.EQ.0)THEN H1(I)=2*F*PI/PI ELSE H1(I)=SIN(2*PI*F*I*dt)/(PI*I*dt) END IF 10 CONTINUE CCCCCCCCCCCCCCCC输出低通滤波因子CCCCCCCCCCCCCCCC OPEN(2,FILE='H1_LOW.DAT',FORM='FORMATTED',STATUS='UNKNOWN') WRITE(2,*)(H1(I),I=-50,50) CLOSE(2) CALL CON(X,H1,Y_LOW,100,101,200) CCCCCCCCCCCCCCCC输出滤波后的数据CCCCCCCCCCCCCCCC OPEN(3,FILE='Y_LOW.DA T',FORM='FORMATTED',STATUS='UNKNOWN') WRITE(3,*)(Y_LOW(I),I=51,150) CLOSE(3) CCCCCCCCCCCCCCCCCC带通滤波器CCCCCCCCCCCCCCCCCCCC DO 20 I=-50,50 IF(I.EQ.0)THEN H2(I)=140 ELSE H2(I)=SIN(2*PI*F2*I*dt)/(PI*I*dt)-SIN(2*PI*F1*I*dt)/(PI*I*dt) END IF 20 CONTINUE CCCCCCCCCCCCCCC输出带通滤波因子CCCCCCCCCCCCCCCCC OPEN(4,FILE='H2_BAND.DAT',FORM='FORMA TTED',STATUS='UNKNOWN')