地震资料数字处理原理及实践

地震资料数字处理原理及实践(PROMAX版)

前言

本书是基于作者多年的数据采集和处理经验写成，书中没有超出前人的任何高深理论，讲述的都是资料处理的常识性知识，在书中作者把理论和实际处理方案以及软件有机的结合在了一起，目的是希望在最短的时间内培养出一个合格的处理员。学习本书虽然不一定使读者立刻成为一个数据处理大师，但从这里起步，无疑会提高大师的成功率。

作为一本从采集到处理细节都兼顾的通用地震资料处理指导书，它既可以作为数据处理的入门书籍，又可以作为一本处理实践教材。多年来,大师们出品的地震资料处理典籍都是着重阐述物探理论,还没有一本做到知行合一,本书则力图做到物探基础理论和处理实践相结合,为读者完整介绍常规处理过程。书中的每段文字、每张图片都是作者亲手打造，前辈的经典理论不会改变，改变的是读者如何解读并把它们付诸实践的过程。

本书采取地震勘探原理和实际地震资料处理相结合的方法，详细讲解了地震资料处理环节中遇到的各种基本物探原理，提供了基础应用实例。采取理论+实例+软件实现一体化的学习方法，避免了以往勘探原理中那种只片面强调理论而不注重实践的情况，使读者在使用Landmark专业的处理软件Promax过程中，就能够掌握基础物探理论，做到学即可用，学则有成。

再吹一下牛皮，本书虽然只是讲解地震资料处理的基础理论，但绝不是以往理论书籍的纯理论讲解，本人从事过多年的一线生产实践，看了太多的理论与实际生产脱节的名著，因此决定写一本物探人员都看得懂、做得来资料处理书，让有志从事地震资料处理的技术人员看过此书后，不需要其它培训资料就可以从事此行业了。

书中的许多观点和方法都是大多数处理人员实际工作中经常面对，而又没有专业资料详细讲解的，因此本人力求以内业人员的角度解决这些问题，提供通俗易懂的解决方案。

以上的内容均与本书实体内容无关，以下正文部分如果与某些前辈的书籍雷同，就当作者是原创，同时笔者也万分高兴，终于写出和大师同样的文字了。

2013年4月2日

第一节坐标系统………………………………………………………………………………1-7 第二节野外地震数据采集……………………………………………………………………8-16 第三节野外观测系统………………………………………………………………………17-22 第四节常用地震波分析…………………………………………………………………… 23-28 第五节常规二维处理………………………………………………………………………29-146 第六节常规三维处理……………………………………………………………………147-169 第七节处理概念讨论……………………………………………………………………170-181 附录地震数据的收集、提交成果……………………………………………………… 182-183

第一节坐标系统

我们生活在地球上，据可靠情报说地球是一个两极稍扁，赤道稍鼓的椭球体，平均半径在6371km。地球表面有山有水有河流，疙疙瘩瘩高低起伏，这些起伏相对于地球的半径来说，真是小菜一碟，如果把地球捏成一个半径1m的球，8800多米高的珠穆朗玛峰的高度也就1.4毫米，所以高低对地球直径来说，不是问题，问题是我们生活在地球上。

那么我们如何知道自己处在地球上的那个位置呢？那就要研究科学术语坐标了。坐标有很多种，传说有直角坐标、极坐标和球面坐标系等，虽然品种较多，但最终目的还是为了确定某一个点的位置，回答“我在那”的问题，至于回答“我是谁”的问题，那是哲学家的事。

1.经纬度坐标

现在我们研究一下球面坐标系，在一个球面上怎样定义我家老屋的位置呢?这很简单，使用经纬度啊。在地球表面，平行于赤道的一系列圆构成纬线，从赤道开始到南极或北极在球面上正好是90°，所以南半球就叫南纬，北半球呢？自己猜猜看吧。连接地球两极和赤道上某个点可以有唯一的一个平面，过这三个点做一个圆，这是半径最大的圆了，这些大圆中有一个正好通过格林威治天文台，有人把那个半个圆弧就叫做0°经线。0°经线以西叫西经，以东叫东经，球体只有360度，所以西经和东经最大也就是180°。这样地球上每个点的经纬度再加上地球半径就对应了唯一的一个物理位置了,这也太复杂了吧。

现在，当你告诉客人你家的屋顶坐标是北纬29°58'44.79"、东经31°8'4.91"，来访者就会想了，如果我从大英帝国的格林威治天文台开始出发，我要沿着赤道向东大约走1/12个圆周，然后再沿着经线也走1/12个圆周，差不多就找到你了。如果你的客人对地球不了解，或者缺少空间想象能力，他肯定会迷路的。

这种定义坐标的方法真是太复杂了，而大多数人一生都生活在一个固定的区域，都认为自己的那片土地是平的。例如你问一个老大爷，去旮旯胡同怎么走，他会告诉你向东走到第二个街口左转再走100米，这时他老人家使用的肯定不是经纬度坐标。我认为他使用的是直角坐标系，也就是欧几里得坐标系统，在平面上以相对于原点的距离XY来定位的。地球很大，地球表面一定范围内的区域都可以看成一个平面的。即使是山区，也是在一个平面上隆起的高地。实际生活中，多数人也是这样做的，都是依赖某个点为起点算距离的，考虑高程的不太多，除非跳楼的和跳水的。

事实上，地球确实是一个不规则球体。宇航员们已经把拍到的照片发到网上了，我还是有点不相信自己生活在一个球上，这个球一转动我岂不是大头朝下了？然后他们又把陆地分为高原、山地、丘陵、平原、盆地等；海底又可以分为大陆架、大陆坡、深海平原、海底山脉等；这么复杂的地球如何以数字形式描述其上某一点的高程呢？

2.椭球体

一提到复杂难解的问题，数学家们就会眼睛发亮了，他们总会把简单的问题搞得更为复

杂，而地球坐标这个问题，那更要命了。首先要把地球模型简化，地球表面十分接近于一个规则的旋转椭球体，即一个椭圆绕着它的短轴（南北两极）旋转而成的旋转椭球体，我们称之为地球椭球，实际上就是一个大肚子的扁扁球。人为构造理论的椭球体就叫做参考椭球体，它忽略了地球表面的高低不平状态，这为描述地球表面的形态提供了一个有效的方法。大地测量学中, 参考椭球体是一个数学上定义的地球表面，它近似于大地水准面。大地水准面吗？大地水准面是一个假想的由地球自由静止的海水平面扩展延伸而形成的闭合曲面。通常是被认为是地球真实轮廓。其形状和大小非常接近自然地球的形状和大小，并且位置比较稳定，因此在大范围的区域内，一般选取大地水准面作为野外测量成果的共同基准面。但不要把大地水准面想象成一个光滑的球面，实际上大地水准面不是一个规则的曲面。因为重力线方向并非恒定的指向地心，导致处处与重力线方向正交的大地水准面也不是一个规则的曲面。大地水准面实际上是一个起伏不平的重力等位面，如果大地水准面上面没有摩擦力，你在大地水准面上滑动就不需要克服重力做功。

由于参考椭球体相对简单，大地控制网的计算和显示点坐标（如纬度，经度和海拔）首选就是这个几何模型。通常所说地球的形状和大小，实际上就是以参考椭球体的长半轴、短半轴和扁率来表示的。目前最常用的参考椭球，是美国国防部制图局(DMA)在1984年构建的WGS-84。椭球体作为一个描述地表的数学模型，历史上曾经使用过多种参考模型，它们的长轴和短轴稍有差异。

下表列出了一些最常见的参考椭球：

有了一个椭球体，如何定义地表某个点的坐标位

置呢？在球面上首先考虑经纬度坐标，也就是地理坐

标系统。这好像又回到前面研究的经纬度坐标上了，

对了，球体上的坐标都和经纬度坐标有着直系的亲缘

关系，不搞清经纬度坐标，其他的坐标怎么来的很难

说清。对以参考椭球面为基准面的坐标，地面点P 的

位置用大地经度L 、大地纬度B 和该点沿着法线方向到

椭球面的距离(大地高)H 表示。当点在参考椭球面上

时，仅用大地经度和大地纬度表示。啥叫大地经度呢？

图1-2-1 参考椭球体

通过地球表面上某个点与两极构成的平面叫做子午面，P点大地的经度就是通过该点的子午面与本初子午面(通过英国格林威治天文台的子午面)之间的夹角；P点的纬度就是通过该点大地椭球体的法线与赤道面的夹角；大地高是地面点沿法线到参考椭球面的距离。所以地表上任意一个位置都可以使用(纬度B，经度L，大地高H)表示，这就和立体几何中(x,y,z)表示空间任意点的坐标类似了。也就是笛卡尔直角坐标系统。

那我们现在是否有办法使用笛卡尔直角坐标系统表示地球上任意点的坐标呢？答案是肯定可以的，但道路依然曲折。

知道了椭球体是什么东西，还要知道另一个概念：基准面。基准面是利用特定椭球体对特定地区地球表面的逼近。上面的定义中出现了两个特定，第一个是指存在很多参考椭球体，每个大国有能力都要搞一个自己的椭球，这样大概能使敌人摸不清自己的藏身之处吧；第二个特定是指在给定的椭球体下产生的自己关心的地球表面区域的近似地表数学模型。

椭球体与基准面之间的关系是一对多的关系，也就是基准面是在椭球体基础上建立的，同样的椭球体能定义不同的基准面，如前苏联的Pulkovo 1942、非洲索马里的Afgooye基准面都采用了克拉索夫斯基椭球体，但它们的基准面是不同的。我们通常称谓的北京54坐标系和西安80坐标系实际上指的是我国的两个大地基准面。我国参照前苏联从1953年起采用克拉索夫斯基(Krassovsky)椭球体建立了我国的北京54坐标系，1978年采用国际大地测量协会推荐的1975地球椭球体建立了我国新的大地坐标系，西安80坐标系。目前最长使用的就是WGS84坐标系统，WGS1984基准面采用WGS84椭球体，它是一地心坐标系，即以地心作为椭球体中心，目前GPS测量数据多以WGS1984为基准。

3.地图投影

上述分析费了很大的气力，空间点的坐标还是在球面上，现在就要把球面上的数据转换到直角坐标系上来了。地图投影就是干这个事的，其科学的定义是：运用一定的数学法则，将地球椭球面的经纬线网相应地投影到平面上的方法。即将椭球面上各点的地球坐标变换为平面相应点的直角坐标的方法。定义表明了坐标投影实质上是将地图从球面转换到平面的数学变换，数学变化当然就有很多种，其中最常用的有墨卡托投影、高斯克里格投影和兰伯特等角投影等。

3.1墨卡托(Mercator)投影

墨卡托(Mercator)投影，是一种"等角正切圆柱投影”，荷兰地图学家墨卡托在1569年拟定, 假设地球被围在一中空的圆柱里，其标准纬线（赤道）与圆柱相切接触，然后再假想地球中心有一盏灯，把球面上的图形投影到圆柱体上，再把圆柱体展开，这就是一幅选定标准纬线上的“墨卡托投影”绘制出的地图。

墨卡托投影特点有两个特点：

A、是没有角度变形，由每一点向各方向的长度比相等，它的经纬线都是平行直线，且相交成直角，经线间隔相等，纬线间隔从标准纬线向两极逐渐增大。

B、墨卡托投影的地图上长度和面积变形明显，但标准纬线无变形，从标准纬线向两极变形逐渐增大，但因为它具有各个方向均等扩大的特性，保持了方向和相互位置关系的正确。

在地图上保持方向和角度的正确是墨卡托投影的优点，因此墨卡托投影地图常用作航海图和航空图，如果循着墨卡托投影图上两点间的直线航行，方向不变可以一直到达目的地，因此它对船舰在航行中定位、确定航向都具有有利条件，给航海者带来很大方便。

3.2高斯-克吕格投影和UTM投影

3.2.1高斯-克吕格投影

此投影是一种“等角横切圆柱投影”。德国数学家、物理学家、天文学家高斯（Ｇauss，1777～1855）于十九世纪二十年代拟定，后经德国大地测量学家克吕格（Johannes Kruger，1857～1928）于1912年对投影公式加以补充，故得名。

高斯是一个数学天才，他被认为是最重要的数学家之一，有数学王子的美誉，被誉为历史上最伟大的数学家之一，和阿基米德、牛顿、欧拉同享盛名。18岁的高斯就发现了质数分布定理和最小二乘法，他还独立发现了二项式定理的一般形式，发现概率统计中的标准正态分布（或高斯分布）函数，严格证明了每一个N阶的代数方程必有N个复数解。

高斯是现代大地测量学的奠基人，除了当数学家，1818年至1826年间，高斯找到了一份大地测量工作，白天他到野外测量数据，夜晚计算。五六年间，经他亲自计算过的大地测量数据，超过100万次。作为测量组长，当野外观测走上正轨后，高斯就在屋里处理观测成果，并写出了近20篇对现代大地测量学具有重大意义的论文。在这些论文中，推导了由椭圆面向圆球面投影时的公式，并做出来详细证明。当然，这些论文也为他赚了不少稿费。后来，德国大地测量学家克里格于 1912年对高斯投影公式加以补充，就成为现在的高斯-克里格投影系统。

其原理如图（图1-3-1）设想用一个圆柱横切于球面上投影带的中央经线，按照投影带中央经线投影为直线且长度不变和赤道投影为直线的条件，将中央经线两侧一定经差范围内的球面正形投影于圆柱面。然后将圆柱面沿过南北极的母线剪开拉平，即获高斯一克吕格投影平面。

通俗的讲就是把地球视为球体,假想一个平面卷成一个横圆柱面并把它套在球体外面,使横轴圆柱的轴心通过球的中心,球面上一根子午线与横轴圆柱面相切。这样,该子午线在圆柱面上的投影为一直线,赤道面与圆柱面的交线是一条与该子午线投影垂直的直线。赤道面在圆柱面上的投影也是一条直线，将横圆柱面展开成平面,由这两条正交直线就构成高斯-克吕格平面直角坐标系。展开的投影带中赤道投影是Y轴，连接南北两极的经线构成X轴。

高斯-克吕格投影的特点是除中央经线和赤道为直线外，其他经线均为对称于中央经线的曲线。投影没有角度变形，在长度和面积上变形也很小，中央经线无变形，自中央经线向投影带边缘，变形逐渐增加，变形最大处在投影带内赤道的两端。由于其投影精度高，变形小，而且计算简便（各投影带坐标一致，只要算出一个带的数据，其他各带都能应用），因

此在大比例尺地形图中广泛应用，可以满足军事上各种需要，并能在图上进行精确的量测计算。

按一定经差将地球椭球面划分成若干投影带,这是高斯投影中限制长度变形的最有效方法。分带时既要控制长度变形使其不大于测图误差，又要使带数不致过多以减少换带计算工作，据此原则将地球椭球面沿子午线划分成经差相等的瓜瓣形地带,以便分带投影。通常按经差6度或3度分为六度带或三度带。

6°带分法是从格林威治零度经线起，每6°分为一个投影带,全球共分为60个投影带。东半球从东经0°-6°为第一带，中央经线为3°，依此类推,投影带号为1-30。其投影代号n 和中央经线经度L0的计算公式为：L0=(6n-3)°。

西半球投影带从180°回算到0°，编号为31-60,投影代号n 和中央经线经度L0的计算公式为：L0=360-(6n-3)°。

三度带是在六度带的基础上分成的，它的中央子午线与六度带的中央子午线和分带子午线重合，即自 1.5度子午线起每隔经差3度自西向东分带，带号依次编为三度带第 1、2…120带。

图1-3-1 高斯-

克吕格投影示意

图1-3-2 高斯-克吕格投影带分带示意图

地震勘探复习资料

绪论 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、剩余时差：当采用一次波的正常时差公式进行动校正之后，除了一次反射波之外，其他类型的波仍存在一定量的时差，我们将这种进过动校正后残留的时差叫做剩余时差。

地震等自然灾害应急预案及处理流程

地震应急预案及处理流程为加强我院安全生产工作，做好安全生产和灾害事故应急处理工作，保护人民的生命和财产安全，根据《中华人民共和国安全生产法》和《灾害事故医疗救援工作管理办法》、参照《全国救灾防病预案》、《国家突发公共事件医疗卫生救援应急预案》和《医疗卫生机构灾害事故防范和应急处置指导意见》有关规定，结合我院实际，制定本预案：一、指导思想根据有关规定和我院安全工作的总体部署，切实做好地震等灾害事故各项准备工作，当破坏性地震发生后迅速启动本预案，统一部署，紧急处置，迅速全面地做好各项抗震救灾准备，高效、有序地开展应急自救工作，以最快速度恢复医疗工作正常开展，最大限度减轻地震灾害，减少人员伤亡和经济损失。二、组织机构 1、指挥部总指挥：院长（党支部书记）副总指挥：业务副院长成员：保卫科、后勤科、医务科、护理部、各临床科室主任职责：

(1)统一领导，健全组织，强化工作职责，加强对破坏性地震及防震减灾工作的研究，完善各项应急预案的制定和各项措施的落实。 (2)充分利用各种渠道进行地震灾害知识的宣传教育，组织、指导医院防震抗震知识的普及教育，广泛开展地震灾害中的自救和互救训练，不断提高广大医务人员防震抗震的意识和基本技能。 (3)认真做好各项物资保障，严格按预案要求积极筹储、落实食品饮用水、防冻防雨、医疗器械、抢险设备等物资，强化管理，使之始终保持良好状态。 (4)破坏性地震发生后，采取一切必要手段，组织各方面力量全面进行抗震减灾工作，把地震灾害造成的损失降到最低点。 (5)调动一切积极因素，迅速恢复正常医疗秩序，全面保证和促进社会安全稳定。指挥部设在院办，电话： 2、疏散组：组长：保卫科科长组员：各临床、医技科室主任、护士长职责：平时负责全院地震等自然灾害培训演练的具体工作，保持疏散通道畅通。 (1)现场指挥，迅速组织医务人员指导患者及其家属按照

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

地震资料处理/解释大作业（处理部分）专业：勘查技术与工程班级：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.道均衡：是指在不同或同一地震记录道建立振幅平衡。 2.数字信号：相对于模拟信号，记录瞬间信息的离散的信号。模拟信号：随时间连续变化的信号. 有效信号：能为我们所利用的信号就叫有效信号。 3.最小相位：能量集中在序列前部。 4.反射波：在波速突变的分界面上，波的传播方向要发生改变，入射波的一部分被反射，形成反射波。折射波：滑行波在传播过程中也会反过来影响第一种介质，并在第一种介质中激发新的波。这种由滑行波引起的波，叫折射波。 5.共深度点：CDP。地下界面水平时，在共中心点下方的点，界面倾斜时无共深度点。 6.解编：地震数据是按各道同一时刻的样点值成列排放的，解编就是将数据重排成行。 12. 最大相位：能量集中在序列后部。 16.地震波：地震波是在岩石中传播的弹性波。多次波：在地下经过多次反射接收到的波叫多次波。 17. 切除：地震信号经动校正后被拉伸畸变，目前处理动校正拉伸畸变的方法是切除，即把拉伸严重部分的记录全部充零。 18. 混合相位：能量集中在序列中部。自相关：一个时间信号与自身的互相关。互相关：一个时间信号与另一个时间信号的相关。 21.环境噪音：交流电、人、风吹草动等环境因素所引起的对地震波有干扰的信号。随机噪音：交流电、人、风吹草动等随机因素所引起的对地震波有干扰的信号。 22.反射系数：反射振幅与入射振幅的比值。 28.模拟记录：把地面振动情况，以模拟的方式录制在磁带上。二、简答题 1、地震资料数字处理主要流程？地震资料的现场处理主要包括哪些内容? 地震勘探资料数据处理中的预处理主要包括哪些内容? 简述地震资料数据中有哪些目标处理方法? 地震资料数字处理如何分类? 地震资料数字处理质量控制有哪些? 地震资料数字处理主要流程：输入→定义观测系统→数据预处理（废炮道、预滤波、反褶积）→野外静校正→速度分析→动校正→剩余静校正→叠加→偏移→显示。地震资料的现场处理主要有：预处理、登录道头、道编辑、切除初至、抽道集、增益恢复、设计野外观测系统、实行野外静校正、还可以进行频谱分析、速度分析、水平叠加等(2分)。地震勘探资料数据处理中的预处理主要包括登录道头、废炮道编辑、切除初至、抽道集(4分)、增益恢复、预滤波、反褶积等. 地震资料数据中目标处理方法有高分辨率地震资料处理、三维地震资料处理、叠前深度偏移处理、井孔地震资料处理(4分)、多波多分量地震资料处理、时间推移地震资料处理等地震资料数字处理分类有数据预处理、数据校正、叠加和偏移归位、振幅处理、滤波、分析、正反演、复地震道技术等。（3分）地震资料数字处理质量控制包括野外原始资料检查与验收、处理流程及主要参数确定、

石油地震勘探资料处理

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

地震数据处理重点

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.野外观测系统定义：观测系统就是以野外文件号和

《地震资料数字处理》复习

《地震资料数字处理》复习地震资料数字处理围绕以下三方面工作： 1、提高信噪比； 2、提高分辨率； 3、提高保真度。一、提高信噪比的处理 1、原理利用噪声和信号在时间、空间、频率和其他变换域中的分布差异，设计滤波因子，将噪声进行压制。 2、处理顺序提高信噪比包含消除噪声和增强信号两部分内容。消除噪声一般在叠前的各种道集上进行，主要针对规则干扰如多次波和面波等，增强信号一般在叠后剖面上进行，主要针对随机噪声。 3、随机噪声是指没有固定的频率、时间、方向的振幅扰动和震动，其成因大致是来自环境因素、次生因素和仪器因素，其中次生干扰的强度与激发能量有关。随机噪声在记录上表现为杂乱无章的波形或脉冲，在频率上分布宽而不定，在空间上没有确定的视速度。随机噪声的随机性与道间距有关，如果道间距减小到一定程度，许多随机噪声表现出道间的相干性，当道距大于随机噪声的相干半径才表现出随机性。 4、一维滤波器（伪门、Gibbs现象）频率滤波器是根据信号和噪声在频率分布上的差异而设计时域或频域一维滤波算子。它压制通放带以外的频率成分，保留通放带以内的频率成分。 Gibbs现象是由于频率域的不连续或截断误差引起的，通放带和压制带之间设置过渡带可克服此现象，设计滤波器就是控制过度带的形状和宽度。 5、二维滤波器二维滤波是根据有效信号和相干噪声在视速度分布上的差异，来压制噪声或增强信号。通常用来压制低视速度相干噪声，在f-k平面上占据低频高波数区域。二维滤波比较容易产生蚯蚓化现象，而且混波相现象明显，在空间采样条件不满足或陡倾角的情况下受到空间假频的影响，一般常用于压制一些规则干扰，如面波和多次波等。 6、频率-波数域二维滤波实现步骤：（1）把时间和空间窗口里的数据变换到f-k域；（2）在f-k域，通过外科切除，按径向扇形划分压制区C(乘振幅置零)、过渡区S（乘振幅置0至1变化）、通放区P (乘振幅置1) ；（3）从f-k域反变换到t-x域。 8、数字滤波有两个特殊性质：（1）数字滤波由于时域离散化会带来伪门现象，

地震勘探资料处理

本科生实验报告实验课程基于 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 初至切除流程模块

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

《地震资料数据处理》课程设计总结报告专业班级：姓名：学号：设计时间：指导老师： 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')

地震勘探资料数字处理

中国地质大学（北京）课程名称：应用地震学教师：段云卿第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 =??= ?

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

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克希霍夫积分偏移考虑了波的振幅值随传播距离和方向不同的影响，保持了波的

(完整版)地震资料数字处理复习题答案

地震资料数字处理复习题一、名词解释（20分） 1、速度谱把地震波的能量相对于波速的变化关系的曲线称为速度谱。在地震勘探中，速度谱通常指多次覆盖技术中的叠加速度谱。 2、反滤波又称反褶积，是指为提高纵向分辨率，去掉大地滤波器的作用，把延续几十至100ms 的地震子波b（t）压缩成原来的震源脉冲形式，地震记录变成反映反射系数序列的窄脉冲组合。 3、地震资料数字处理就是利用数字计算机对野外地震勘探所获得的原始资料进行加工、改进，以期得到高质量的、可靠的地震信息，为下一步资料解释提供可靠的依据和有关的地质信息。 4、数字滤波数字滤波就是指用数学运算的方式用数字电子计算机来实现滤波。对离散化后的信号进行滤波，输入、输出都是离散数据。 5、水平叠加将不同接收点受到的来自地下同一反射点的不同激发点的信号，经动校正叠加起来。 6、叠加速度在一般情况下，都可将共中心点反射波时距曲线看作双曲线，用一个同样的式子来表示：t2=t02+x2/Vα2，其中，Vα就是叠加速度。 7、静校正把由于激发和接收时地表条件变化所引起的时差找出来，再对其进行校正，使畸变了的时距曲线恢复成双曲线，以便能够正确地解释地下的构造情况，这个过程叫做静校正。 8、动校正消除由于接受点偏离炮点所引起的时差的过程，又叫正常时差校正。 9、假频一个连续信号用过大的采样得到的离散序列实际上含有连续信号中高频成分的贡献。这些高频成分折叠到离散时间序列中较低的频率。这种现象是由连续信号采样不足引起的，称作假频。 10、亮点技术所谓“亮点”狭义地说是指地震反射剖面上由于地下油气藏存在所引起的地震反射波振幅相对增强的“点”。利用地震反射波的振幅异常，同时也利用反射波的极性反转、水平反射的出现、速度的降低及吸收系数的增大等一系列亮点标识综合指示地下油、气藏的存在，进而直接寻找油、气藏的技术。 11、相关定量地表示两个函数之间相似程度的一种数学方法。 12、自相关表示波形本身在不同相对时移值时的相关程度。（一个时间信号与自身的互相关） 13、环境噪音由自然条件或环境（如风吹草动、工业交流电的干扰等）造成的对地震波有效信号的干扰。 14、有效信号野外地震工作想要得到的含有地下地质信息的地震信号。 15、振幅振动物体离开平衡位置的最大距离，在数值上等于最大位移的大小。 16、共中心点在不同激发点、不同接收点的记录中具有公共炮检的中点。（野外采用多次覆盖工作方法时，如界面水平，则每次都能接受来自界面上同一点的反射。该点在地面上的投影称为共中心点。） 17、共深度点地下界面水平时，在共中心点下方的点，界面倾斜时无共深度点。 18、反褶积同2反滤波

地震资料解释大作业报告

新疆准噶尔盆地西北缘夏子街地区地震资料解释报告中国石油大学（北京）地球科学学院肖鸿宇

一、概述 1、研究区背景信息本区位于新疆准噶尔盆地西北缘夏子街地区。区域地质研究表明本区主要发育三个大的区域不整合面：白垩系/侏罗系不整合面、侏罗系/三叠系不整合面和三叠系/二叠系不整合面。其基本特征为：①均为角度不整合面②界面之上均为厚层砂砾岩，界面之下岩性较细，甚至于为泥岩。本区砂砾岩波阻抗高，泥岩波阻抗低。本区在区域构造单元上属克——乌逆掩断裂带，发育众多断裂，主断裂与区域构造的走向基本一致，断面一般北西倾斜，它们和方向不一、倾向不同的派生断裂一起将本区切割为若干断块，形成了三叠系以背斜——断块为主的构造圈闭。 2 地震资料概况资料采用新疆准噶尔盆地西北缘夏子街地区部分三维测线，共21条。其中主测线（inline）10条，联络测线（crossline）11条，线距均为500m，比例尺为1:25000。本区剖面经过三维叠加偏移处理。采用SEG标准显示，负反射系数界面为波峰，正反射系数界面为波谷。 3 任务与要求 1)标志层和不整合面的识别要求：对所有剖面进行观察，确定本区存在几个标志层。识别出三个大的区域不整合面，确定其是波峰还是波谷。确定在三叠系内部和侏罗系内部有几个标志层。确定各标志层的地层接触关系和性质。 2)断层解释和主要界面的追踪对比要求：对21条剖面进行断层解释，对三大区域不整合面和三叠系内部、侏罗系下部标志层共5个界面进行追踪对比。要求断层和界面在所有的剖面交点上闭合

3)三叠系底界T0构造图编制和构造特征分析要求：剖面上断点投影到平面底图上，并且标出断层的性质、倾向。断层组合。T0图勾绘，等值线间距为20ms。分析主要断层和褶皱的性质及形成时期 4)地震相分析（编制侏罗系下部层序地震相平面图）要求：对侏罗系下部层序进行剖面地震相划分。编制沉积相平面图。进行的初步的沉积相解释。二、地层格架 1 标志层及其特征地震反射标志层是指波形特征突出、稳定且分布广泛的同相轴或波组。通过对地震资料的分析，确定以下4个标志层，以剖面Y 360（附图1）为例，标志层为T1、T3、T5、T6，标志层T1、T3、T5、T6均延伸较远，几乎全区都有分布，而且特征十分明显突出，横向变化稳定。 T1、T3与上覆、下伏同相轴均呈整一接触关系。T5与上覆同相轴均呈上超接触关系，于下伏同相轴呈整一接触关系，T6与上覆同相轴呈整一接触关系，与下伏同相轴呈削截接触关系。 2 不整合界面及其特征不整合面是指具有不整合接触关系的两套地层之间的接触面。以剖面Y 360（附图1）为例，不整合面T1、T4、T6将岩层分别划分为白垩系、侏罗系、三叠系和二叠系，T1、T6同时也为标志层。均为全区分布，而且不整合面延伸较远，特征十分明显。如果通过横剖面进行观察，以剖面X433（附图2）为例，则可以发现在横剖面上T1、T4延伸较远，T6则较短。T1与上覆同相轴呈上超接触关系，与下伏同相轴呈整一接触关系。T4与上覆同相轴呈整一接触关系，与下伏同相轴呈削截接触关系。T6与上覆同相轴呈整一接触关

地震资料处理合同(基地内)

说明一、起草单位与起草人本合同由中国石油化工集团公司法律事务部组织起草，由胜利石油管理局负责编写，起草人：加庆段清海，联系：0，传真：0，电邮：duanqhslof.。二、注意事项 1、本合同适用基地场景,发包方、承包方均为系统且在同一基地的单位所发生的地震资料处理业务。 2、本合同的修改。修改本合同不影响甲方实质性权利义务的，应经甲方兼职合同管理员审查同意。修改本合同影响甲方实质性权利义务的，应经甲方专职合同管理员审查同意。 3、具体条款使用说明。（1）地震资料处理合同示文本作为一个整体，其部的各条款容之间是具有关联性的，在实际应用过程中如对个别条款做出变动，那么其相对应的条款也要做出相应的调整。如：要调整双方权利义务的条款容，在与之相对应的违约责任条款中也要改动相应的容。（2）文本中质量标准和技术要求条款的规定，应结合实际针对不同井的具体情况，选择、引用明确的标准，并把该质量标准详细列明作为本合同的组成部分。（3）文本中HSE条款对甲方、乙方在安全、环保、健康方面做出了原则性的要求和规定，在实际操作中可以引用HSE方面的法规或相关规定执行，或双方另行签订HSE责任书将容细化，并作为合同的附件双方共同遵守。（4）文本中的价款支付方式和费用的调整，可根据具体项目的不同和本单位的习惯性做法，在与乙方协商一致后做出调整。（5）违约金的约定在文本中都是以“空格”的形式列出的，在实践过程中应根据具体情况协商做出约定。（6）违约责任条款中关于赔偿限额的规定，参考国同行业在此问题上的惯例，制定出一个客观的、合理的赔偿额度。（7）文本中有关“时间”、“期限”的要求，在实际填写中应结合生产实际，按照地震资料处理服务的工序、要求制定出合理的时间和期限。合同编号：

地震资料处理复习总结

数字地震记录中，每个地震到是按一个按一定时间采样间隔排列的时间序列数字滤波，每一个地震道都可以用一系列具有不同频率和不同振幅、相位的简谐曲线叠加而成。应用一维傅里叶变换可以得到每个地震道德简谐成分；应用傅里叶反变换可以将简谐成分合成为原来的地震道的时间序列函数。通常由傅里叶变换得到的频谱为一个复函数，称为复数谱。它可以写成指数形式式中复数的模，称为振幅谱; 复数的幅角，称为相位谱。离散情况下和这个差不多一维频谱的特征： 1．傅里叶变换的几个基本性质线性翻转共轭时移褶积相关（功率谱） 2．假频尼奎斯特频率二维谱分析二维波场函数X(x,t)的二维傅里叶变换° X(,)ωκ 表明了二维波场函数X(x,t)的各个频率f 一波数简谐成分的频一波谱。由°X(,)ωκ这些频率f 一波数的简谐成分叠加即可恢复原来的波场函数X(x,t)。二维傅里叶变换X(w,k)称为二维函数X(x,t)的频一波谱。其模量为函数X(x,t)的振幅谱。如果有效波和干扰波得平面简谐波成分有差异，有效波的平面简谐波成分与干扰波的平面简谐波成分不同的视速度传播，则可以用二维视速度滤波将他们分开，达到压制干扰，提高性噪比的目的。二维频谱的特征：空间假频 ~~ () ()()()()i w i w X w X w e A w e ??==)(ωA ()?ω1()()tan () i r x w w x w ?- =()A w =t f ?=21N o k o k ~ X(,)k ω

在地震勘探中，用数字仪器记录地震波时，为了保持更多的波得特征，，通常利用宽频带进行记录，因此在宽频带范围内记录了各种反射波的同时，也记录了各种干扰波。有效波和干扰波得差异表现在多个方面（频谱、传播方向、能量……）。利用频谱特征的不同来压制干扰波，以突出有效波的方法就是数字滤波。滤波器的响应特性：对滤波器能力的最普遍度量是其响应特性滤波器的频率特性：其滤波器时间函数或滤波因子的频谱称为滤波器的频率特性，滤波器的时间特性（单位脉冲响应）：在时间域的表示方法中，令一个单位脉冲通过一个滤波器，然后观侧滤波器的输出，这个滤波器输出的自然过程曲线称为滤波器的时间特性。也称滤波器的“脉冲响应” 频率响应函数应该就是时间和频率响应函数合起来应该就是就是响应特征滤波机理：输出信号的振幅谱等于输入信号的振幅谱与滤波器的振幅频率特性的乘积，输出信号的相位谱等于输入信号的相位潜与滤波器相位特性之和。（频率）时间域上就是褶积褶积滤波的物理意义：它相当于把地震信息分解为起始时间、极性、幅度各不相同的脉冲序列，令这些脉冲按时间顺序依次通过滤波器，这样在滤波器的输出端就得到对输入脉冲序列的脉冲响应，这些脉冲响应有不同的起始时间，不同的极性和不同的幅度(这个幅度是与引起它的输入脉冲幅度成正比的)，将它们叠加起来就得到滤波后的输出分。频率域滤波的步骤 (1)对已知地震记录道进行频谱分析。 (2)设计合适的滤波器 (3)进行滤波运算 (4)对输出信号谱进行傅里叶反变换褶积滤波的具体计算褶积滤波的具体计算步骤如下： (1)对地震记录进行频谱分析，确定通频带中心频率和带宽。 (2)确定滤波因子长度N 。 )()()(~ ~~w H w X w X =)(t x ∧ )(~ w H )(t h )(~ w H ) ()()(w w w H x x Φ+Φ=Φ∧)()()(~ ~w H w X w X ?=∧)(~ w X )(t x