绪论 1、地球物理勘探的概念 （1）简称“物探”，是通过观察存在地球及其周围的地球物理场的特征和岩石的各种物理特性来研究地质规律和勘查各种矿产的各种方法的总称。（2）是以物理学原理为基础，利用电子学、计算机的数字处理、信息论等科学技术中的新技术所建立起来的一整套勘探地下矿产的方法。（3）是借助于各种物探仪器在地面观测地下岩石的各种物理参数，从而解释和推断地下岩石的构造特点、岩石性质等，从而到达勘查地下矿产（金属非金属矿产、煤、油气等等)的目的。 2、地球物理勘探的分类，不同勘探方法的优缺点。 重力勘探：利用岩石的密度差异 磁法勘探：利用岩石的磁性差异 电法勘探：利用岩石的电性差异 地震勘探：利用岩石的弹性差异 放射性勘探：利用岩石的放射性差异 地震勘探的优点：精度高，分辨率高，穿透深度大，能较详细地了解由浅至深一整套地层的地质规律。缺点：成本高 3、地震勘探的概念、分类，目前地震勘探以何种方法为主。 概念：利用岩石的弹性差异来进行矿产勘察。是通过人工激发地震波，研究地震波在弹性不同的地下地层中传播的规律，以查明地下的地质构造，达到油气或其他勘探目的的一种物探方法。 分类：地质法（优：在找油初期，可以起到一个指向作用，避免了盲目性，成本低。缺：野外地质方法很难准确了解地下地质情况！）；钻探法（优点：精度最高，缺点：一孔之见，而采用大量的钻井，不仅成本高，而且效率低）；物探方法（优点：精度高于地质法，成本低于钻探法；不足：精度低于钻探法，成本高于地质法）。 应用最多的方法：物探方法 4、地震勘探的三个阶段 地震资料野外采集、地震资料室内处理、地震资料解释。 第一章 各种介质的概念 重点：①物体是否为弹性、塑性介质与受力大小、时间及温度有关。②均匀介质与各向同性介质的关系。 （1）理想弹性介质：当介质受外力后立即发生形变，而外力消失后能立即完全恢复为原状的介质； （2）粘弹性介质：当外力消失后不是立即恢复原状，而是过一段时间后才恢复原状的介质称为粘弹性介质。 （3）塑性介质：当外力消失后不能完全恢复原状,保留了一部分形变的介质称为塑性介质。（4）各向同性介质：凡介质的弹性性质与空间方向无关的介质称为各向同性介质 （5）各向异性介质：凡介质的弹性性质与空间方向有关的介质称为各向同性介质 （6）均匀介质：弹性性质（波速）不随空间坐标的变化而变化，是常数。 （7）非均匀介质：弹性性质（波速）随着空间坐标的变化而变化，不是定值。 （8）层状介质：如果非均匀介质的物理性质呈层状分布，则称这种介质为层状介质。层状介质中各层的弹性系数是不变的。层状介质模型已经成为地震勘探中常用的物理 模型。 （9）连续介质：层状介质的层数无限增加，每层的厚度无限减小时，层状介质就可以视
安徽理工大学 一、名词解释（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反滤波
说明 一、起草单位与起草人 本合同由中国石油化工集团公司法律事务部组织起草，由胜利石油管理局负责编写，起草人：加庆段清海，联系：0，传真：0，电邮：duanqhslof.。 二、注意事项 1、本合同适用基地场景,发包方、承包方均为系统且在同一基地的单位所发生的地震资料处理业务。 2、本合同的修改。修改本合同不影响甲方实质性权利义务的，应经甲方兼职合同管理员审查同意。修改本合同影响甲方实质性权利义务的，应经甲方专职合同管理员审查同意。 3、具体条款使用说明。 （1）地震资料处理合同示文本作为一个整体，其部的各条款容之间是具有关联性的，在实际应用过程中如对个别条款做出变动，那么其相对应的条款也要做出相应的调整。如：要调整双方权利义务的条款容，在与之相对应的违约责任条款中也要改动相应的容。 （2）文本中质量标准和技术要求条款的规定，应结合实际针对不同井的具体情况，选择、引用明确的标准，并把该质量标准详细列明作为本合同的组成部分。 （3）文本中HSE条款对甲方、乙方在安全、环保、健康方面做出了原则性的要求和规定，在实际操作中可以引用HSE方面的法规或相关规定执行，或双方另行签订HSE责任书将容细化，并作为合同的附件双方共同遵守。 （4）文本中的价款支付方式和费用的调整，可根据具体项目的不同和本单位的习惯性做法，在与乙方协商一致后做出调整。 （5）违约金的约定在文本中都是以“空格”的形式列出的，在实践过程中应根据具体情况协商做出约定。 （6）违约责任条款中关于赔偿限额的规定，参考国同行业在此问题上的惯例，制定出一个客观的、合理的赔偿额度。 （7）文本中有关“时间”、“期限”的要求，在实际填写中应结合生产实际，按照地震资料处理服务的工序、要求制定出合理的时间和期限。 合同编号：
1地震时间剖面的显示方式为波形加变面积图。 2在地震资料预处理中，数据解编是把按时分道排列的野外磁带数据重新排列转换成按道分时排列的形式，方便对数据作实质性处理。 3数字滤波处理时地震资料处理中提高地震记录信噪比的重要方法。 4沿测线闭合追踪同一反射层位时，To应该闭合，当闭合差超过0.5△T时，就认为没能闭合。 5将时间剖面转换成深度剖面，一般采用的是平均速度，也可用层速度逐层求层厚度的方法。 6要想保留有效波，滤去干扰波，需要设计理想滤波器的类型为带通滤波。 7水平叠加处理所需要的道集排列形式是共中心点排列。 8不整合面是地壳升降运动引起的，与油气聚集有密切关系。 9在均匀介质中，水平叠加时间剖面上的一个脉冲，对应地质空间的反射界面是，而地质空间的一个点，对应水平叠加时间剖面的。 10地震子波为零相位子波。 11根据作图等值线的性质不同，地震构造图可分为等深度构造图和等To构造图两大类。 12波动方程偏移是以波动理论为基础的偏移处理方法，主要由波场延拓和成像两部分组成。 13背斜构造断开后，下降盘等值线的范围比同深度上升盘的小。 14地震层序是沉积层序在地震剖面上的反应。 1什么是偏移现象？ 当反射界面水平时，反射波同相轴与地下界面形态一致，当反射界面倾斜时，反射波到同相轴与地下界面形态不一致，若将反射时间作时深转换，所得视界面与地下真实反射界面比较，不论是界面长度、界面位置及界面倾角两者均不一致，视界面相对于界面，想界面下倾方向偏移，而且倾角变小。 2简述断层在时间剖面上的识别标志。 ①标准层反射同相轴发生错断，是断层在地震剖面上表现的基本形式。由于断层规模不 同，可表现为波组或波系的错断，这是中小型断层的反应。 ②标准层反射波同向轴数目突然增减或消失，波组间隔发生突变，断层下降盘地层加厚， 上升盘地层变薄。 ③反射同相轴形状和产状发生突变，这往往是断层作用所致。因断层的屏蔽作用，造成 下盘反射同相轴凌乱，甚至出现空白反射带。 ④标准层反射波同相轴发生分叉、合并、扭曲及强相位转换。 ⑤断面波、绕射波等异常的出现，是判别断层的主要标准。 3简述一维合成地震记录的制作过程。 4反滤波效果的好坏主要受哪些因素的影响？ ①各种反滤波方法都必须有若干假设条件。 ②反射地震记录褶积模型问题 ③噪声干扰的影响 ④原始地震资料的质量问题 5反射波同向轴的主要识别标志有哪些？ ①相位相同②能量增强③波形相似④连续性 6爆炸反射界面形成的原理是什么？ 叠加剖面相当于自激自收剖面，，若将剖面中时间除以2，或传播速度减半，就可将自激自
地震资料的处理方法和结果在很大程度上受野外采集参数的影响。 地震剖面的“三高”：高信噪比、高分辨率和高保真度。 地震资料处理主要有三个阶段；每一个阶段都是为了提高地震分辨率，即分离出两个无论在空间上还是时间上都非常相近的同相轴的能力。 ●（a）反褶积是通过压缩基本地震子波成为尖脉冲并压制交混 回响，沿着时间方向提高时间分辨率； ●（b）叠加是沿着偏移距方向压缩，把地震资料的数据量压缩 成零偏移距剖面，以提高信噪比； ●（c）偏移是一个使绕射收敛，并将叠加剖面上的倾斜同相轴 归到它们地下的真实位置上，通常在叠加剖面（接近于零偏移 距剖面）上做偏移，来提高横向分辨率。 ●几何扩散校正：通过给数据加一增益恢复函数以校正波前（球面）扩散对振幅的影响。 ●建立野外观测系统：把所有道的炮点和接收点位置坐标等测量信息都储存于道头中以保证各道的正确叠加。 ●野外静校正：对陆上资料，把所有炮点和接收点位置均校正到一个公共基准面上以消除高程、低降速带和井深对旅行时的影响。
关于分辩率的讨论: 有一种普遍的误解，认为要增加时间分辨率只需要高频，这是不 真实的。 只有低频或只有高频不能改善时间分辨率。要增加时间分辨率低 频和高频两者都需要。 时间分辨率取决于有效信号的频带宽度. 最小平方法---根据误差的平方和最小来设计滤波器； 最小相位信号是具有对相同振幅谱的物理可实现信号中相位最 小的信号，或者说能量延迟最小的信号。 最小相位滤波器是具有同样振幅响应的一切可能的滤波器簇中 能量延迟最小的滤波器，也称最小延迟滤波器。 若最小相位滤波器的输入是最小相位，则其输出也是最小相位, 对于地震子波，除了零相位子波外，最小相位子波的分辨率最高。 下面的四个子波中哪一个是最小相位的： 子波A ：（4，0，-1） 子波B ：（2，3，-2） 子波C ：（-2，3，2） 子波D ：（-1，0，4） 频率、视波数和视速度的关系为: **=k f V
地震资料格式说明
————————————————————————————————作者：————————————————————————————————日期:
§3 资料处理流程说明: 资料处理的基本流程如下图所示: 解编 预处理 (建立工区，切除,振幅处理等) 一次静校正 一、二维数字滤波 抽道集 高精度速度分析 剩余静校正高精度动校正 水平迭加 滤波、反滤波 (倾斜相干加强) 迭后偏移 一维数字滤波 振幅均衡、输出 在资料的处理过程中,应根据资料的信噪比和分辨率情况选择模块，组合流程，以达到事半功倍的效果。在处理过程中,应首先根据野外电子观测班报和测量电子班报建立工区基本参数文件（由建立工区模块完成)，若无测量结果，可根据模块提示完成建立工区基本参数文件的工作。本系统适合于有或无测量资料的情况;同时也适合于变观资料处理。文件格式参见相关模块说明。 §４处理资料文件格式说明: ４.1 SEＧ-Y 记录格式(标准） (1）卷头: ３６00字节
(a)AＳCII 区域： 3200字节(40条记录 x 8０字节/每条 记录）。 (b)二进制数区域: 400字节(32０1~３６00)。 3213~３214 字节—每个记录的数据道数(每炮道数或总道数)。 3２17~3218 字节—采样间隔(μs)。 3２２1～3222 字节—样点数／每道(道长)。 ３225~３２２６字节—数据样值格式码1－浮点; 3255~32５6 字节—计量系统：1－米, 2-英尺。 ３261~3262*字节—文件中的道数(总道数)。 3２6９~327０*字节—数据域(性质）:０-时域，1－振幅,2-相位谱 “ * “ 号字为非标准定义。 (２)道记录块： (a)道头字区: 含: 60个字/4字节整或12０个字／2字节整, 共24０个字节，按二进制格式存放。 ·SEG—Ｙ格式道头说明: 字号(4字节) 字号(2字节) 字节号内容说明 1 1— 2 1—４一条测线中的道顺序号，如果一条测线有 若干卷磁带，顺序号连续递增。 2 3—4 5—8 在本卷磁带中的道顺序号。每卷磁带的道 顺序号从l开始。 3 5—６ 9—12 * 原始的野外记录号（炮号）。 4 7—８ 13—1６在原始野外记录中的道号。 ５９—１0 17—20 测线内炮点桩号(在同一个地面点有多于 一个记录时使用)。 ６ 1１—12 ２1—24 CMP号（或CDＰ号)。(弯线=共反射面元号) 7 １３—１4 25—28 在CMＰ道集中的道号（在每个CＭP道集 中道号从1开始）。 8—1 1５ 29—３０* 道识别码： l=地震数据； 4=爆炸信号; 7＝计时信号; 2=死道； 5=井口道；８=水断信号; ３＝无效道(空道)；6＝扫描道；9…N ＝选择使用