66第5章地震波处理2010

2020/1/291地震数据处理过程一般分三个阶段:预处理、参数提取和分析、资料处理。处理的最终结果是得到供解释用的水平叠加时间剖面或叠加偏移时间剖面。第一节预处理一、数据解编野外磁带记录数据是按时序排列的，即依次记下每道的第一个采样值，2020/1/292各道记完后，再依次记下各道的第二个采样值，依此类推。•在数据处理中，将按时序排列的形式转换为按道序排列(即第一道的所有数据都排在第二道之前，使同一道数据都排放在一起)这种预处理称为数据解编或重排。2020/1/293•二、编辑•在地震数据采集中，由于施工现场复杂，外界干扰大，难免出现一些不正常道和共炮点记录，这些记录信噪比低，如果参与叠加处理会严重影响处理效果。2020/1/294在正式处理之前，需要对这些不正常的记录进行编辑处理，例如对信噪比很低的不正常道进行充零处理，发现极性反转的工作道对它们进行改正等。另外，还要显示有代表性的记录并观察初至同相轴，以便进行初至切除。2020/1/295切除是为了消除包括噪声的记录开始部分所存在的高振幅，这样做对避免以后处理时出现的叠加噪声有好处。切除的方法就是用零乘需要切除的记录段。•2020/1/296抽道集也叫共深度点选排，是把具有相同炮检距点的记录道排成一组，按共深度点号次序排在一起。抽道集处理后，磁带上记录的次序是以共深度点号（CDP)为次序的记录，以后所有的处理都将方便地以共深度点格式进行。•三,抽道集2020/1/297•在野外数据采集过程中，为了使来自不同深度信号的能量能够以一定的水平记录在磁带上，数字地震仪采用了增益控制，对浅层信号放大倍数低，深层信号放大倍数高。对经过增益控制的地震记录恢复到地面检波器接收到的振幅值的处理称为增益恢复。2020/1/298数字仪对信号进行增益控制时的增益指数己记录在记录格式的阶码上，因此增益恢复的公式为A=A0/2n其中A0为记录到的采样值，A为地面检波器接收到的增益控制前的振幅值，n为阶码(即增益指数)。2020/1/299球面扩散是当波离开震源时由于波前扩散造成的振幅衰减，能量发生扩散，波的强度减小，而波场的总能量不变。如果介质是各向同性的，则能量衰减与传播距离的平方成反比。通常速度都是随深度的增加而增加.非弹性衰减是弹性能量由于摩擦而耗散为热的吸收的结果，波动能量消失。•2020/1/2910参数提取与分析的目的是为寻找在常规处理或其他处理中常用的最佳处理参数，以及有用的地震信息，如频谱分析、速度分析、相关分析等。这类数字处理还可为校正与偏移及各种滤波等处理提供速度和频率信息，并可以自成系统处理出相应的成果图件，如频谱、速度谱，通过相关分析进行相关滤波等。2020/1/2911•地震勘探所得到的记录中包含有效波和干扰波，这些波之间在频谱特征上存在很大差别。为了解有效波和干扰波的频谱分布范围，需要对随时间变化的地震记录讯号进行傅里叶变换，得到随频率而变化的振幅和相位的函数，（地震记录的频谱—振幅谱和相位谱）。对地震波形函数进行傅里叶变换求取频谱的过程叫频谱分析。2020/1/2912•1．面波频谱的峰值低于有效波，声波频谱峰值偏高，与有效波的频谱范围有较宽的重叠;•2．微震干扰波的频带较宽;•3．有些规则干扰波与有效波频谱差异不大，如浅层记录中的外界相干干扰波和多次波•4.横波与纵波相比频谱峰值低，频带窄;•5高速薄层反射波频谱相对厚层要偏高;2020/1/29136．大炸药量激发比小炸药量激发频谱要偏低，小炸药量激发比锤击频谱要宽;7.反射波的频率随着低降速带厚度的增加而降低。当低速带较薄或表层速度较高时，获得的反射波频率较高。2020/1/2914•速度参数在浅层地震资料的数据处理和解释中是非常重要的参数，例如校正、叠加和偏移都需要知道速度。•2020/1/2915图4.6-36可控震源野外记录的相关处理2020/1/2916另外速度参数可提供关于构造和岩、土性质有价值的信息，例如构造探查要了解地下反射界面的分布，实质上是波阻抗参数的地下分布，岩性探查要得到地下岩性的分布，更与各种岩性参数(例如速度、吸收系数、泊松比等)的提取有关。2020/1/2917•由于地下介质的复杂性，故参数提取是一个十分复杂而艰巨的任务，只能用一些简化的方法和近似的假设条件来求取。•所用方法不同，可得到不同定义的速度参数。地震速度分析中普遍采用速度谱分析和速度扫描技术，得到平均速度、均方根速度、层速度等速度参数。•2020/1/2918设共反射点道集内有N个记录道，其炮检距分别为x1，x2，…xn，各个记录道对应的正常时差分别为t1，t2，…tn，炮检距为xi的第i个记录道的正常时差为•2020/1/2919)16.2.4(02220tvxttsii可见对于给定的炮检距，正常时差ti是垂直反射时间t0和叠加速度vs的函数，因而从反射波正常时差t的分析中可以提供均方根速度的信息，这就是速度分析的基础。02202tVxttt2020/1/2920在速度分析中，常用的方法是速度谱和速度扫描。速度谱的概念是从频谱的概念借用而来的，频谱表示波的能量相对频率的变化规律，人们就将地震波的能量相对速度的变化称为速度谱。根据动校正原理，选取一系列试验速度分别代入(4·2·16)求取正常时差ti，2020/1/2921••波时距曲线进行动校正，看其校正以后双曲线形状的同相轴是否变成t=t0的水平同相轴，如果变成水平同相轴，则所采用的速度就是最佳叠加速度.•衡量同相轴是否被拉平，可以选择不同的判别准则，如果以共反射点波组叠加波形的能量来判别，则当选用速度合适时，同相轴被拉平。2020/1/2922形成同相位叠加，叠加波形的能量为最大。按这一原理设计的速度谱称为叠加速度谱。为求取地震波的叠加速度，可以根据地区的地质情况和有关资料，大致估计在某一t0时的速度变化范围为v1,v2,…vm,将每一速度代入(4.2.16)式计算每道的正常时差ti,做动校正，然后把动校正后的各道记录振幅按下式叠加，得相应的平均振幅•2020/1/2923图4.2-4用多次覆盖资料计算速度谱原理图2020/1/2924图4.2-5计算叠加速度谱的网络2020/1/2925图4.2-6三维显示形式的速度谱2020/1/2926•速度扫描是用一组试验速度分别对单张(CDP)道集记录或单次覆盖共炮点记录作速度扫描动校正，即一次用一个试验速度对整张记录上的所有波组进行动校正(恒速动校正)，得到一张校正后的记录。•当所用的某一试验速度正好与某一t0时间所对应的真实速度一致时，此t0时刻的同相轴2020/1/2927•会变得平直，其他同相轴或者上弯(速度过高，校正不足)，或者下弯(速度过低，校正过量)。寻找各试验速度校正记录上的平直同相轴，可以得到不同t0时间处反射波的速度。由速度扫描获得的速度是叠加效果较好的速度.•它适用于地质条件复杂得不到好速度谱的地区，但处理费时长，成本高。2020/1/2928•由于通过速度扫描或速度谱求出的速度反映了叠加效果的好坏，一般称之为叠加速度vs实质上它表示用双曲线拟合有效波时距曲线时，拟合效果最好的速度，故也称之为双曲线拟合速度。•2020/1/2929(三)各种速度之间的关系在水平层状介质中，波沿某一条射线传播时，它传播的总路径与总时间之比就是射线速度vr,)18.2.4()1()1(1212212122niiiiniiirvpvhvphv2020/1/2930•式中hi为第i层厚度，p为射线参数。这是沿一条射线取平均算出的速度。射线不同，vr也不同。由于射线速度不仅考虑了波在界面的“偏折效应”，同时也考虑了横向不均匀的影响，•1.射线速度比其他速度更精确；•2.当射线参数p为零(或炮检为零)时的射线速度即平均速度，2020/1/2931•因此射线速度大于等于平均速度;•3.炮检距为无穷大时的射线速度等于水平层状介质中最高速度层的速度。•4.由于均方根速度是等效均匀层的最佳射线速度，它也考虑了射线在界面的“偏折效应”，它总是大于平均速度。2020/1/29325.水平层状介质情况下，炮检距不太大时的叠加速度就是均方根速度;6.对于单层均匀介质，叠加速度就是介质的真速度；在倾斜界面情况下，叠加速度是等效速度（均方根速度除以界面倾角之余炫）2020/1/2933••1．迪克斯（Dix）公式•水平层状介质情况下，叠加速度就是均方根速度，因此迪克斯（Dix）公式就是由均方根速度求层速度的公式。••2020/1/29342/11,0,01,21,0,2,0nnnnnnnttvtvtv对于倾斜地层，叠加速度就是等效速度，Dix公式写为2/11,0,021,1,02,,0)cos()cos(nnnennennttvtvtv2020/1/2935NoImage第五节校正和叠加处理一、野外静校正包括井深校正、地形校正、低速带校正。见图2020/1/2936NoImage静校正示意图2020/1/2937进行野外静校正处理时，首先对共深度点道进行选排，找出每一道的炮点和检波点的位置，求出相应炮点和检波点的t静校正值，可以从井口记录道上直接读出，称为井口值。因波从O向下传播少用了时间T，校正时要把此值加到波的旅行时间中.2020/1/29381．井深校正井深校正是将激发源0的位置由井底校正到地面0j(见图)，其方法有二:2020/1/2939NoImage1)在井口埋置一井□检波器，记录直达波由0传至地面0j的时间j，即井深校正值，又称为井口时间。2)用已知的表层参数及井深数据，按下式计算井深校正量)1)(100hvhhvjj2020/1/2940NoImage式中V0是低速带波速，V为基岩波速，h0+hj为炮井中低速带厚度，h是基岩中炸药埋置深度。因为井深校正总是向时间增大的方向校正，故此式前面取负号。2020/1/2941NoImage2．地形校正地形校正是将测线上位于不同地形处的炮点和检波点校正到基准面上。如图6.4-5所示，炮点地形校正量为0001hv2020/1/2942NoImage而检波点地形校正量是sshv01此道(第j炮第l道)总的地形校正量为)(100sijhhv2020/1/2943NoImage地形校正量有正有负，通过h0、hs的正负体现出来。通常规定当测点高于基准面时为正，低于基准面时为负2020/1/2944NoImage3．低速带校正低速带校正校正是将基准面下的低速层速度用基岩速度代替。求取低速带校正量的公式在炮点处为:2020/1/2945NoImage)11(0'vvhjj在检波点处为)11(0'vvhll故此道总的低速带校正量为))(11(0'ljjihhvv2020/1/2946NoImage因为基岩速度大于低速带速度，低速带校正量总是正的,第j炮第l道的总的校正量为)(1)(10ljlshhhvhhvt静2020/1/2947NoImage(3)静校正的实现由于静校正值有正有负，校正时则使记录道样值可能向前(向小时间方向)或向后(向大时间方向)在计算机内存中搬家。2020/1/2948NoImage如果记录校正值为+4个采样间隔，则该道记录全部样值要向前般动四个单元，搬动时要从小时间的样值开始并依次搬，记录道最前面的四个样值校正后被冲掉，2020/1/2949地震道的记录样式2020/1/2950NoImage结果第一个样值就是原记录道上的第五个样值，然后把尾部的四个单元冲零。当静校正值是-4个采样间隔时，则该记录道的全部样值要依次向后面搬动四个单元，2020/1/2951NoImage开始搬时将倒数第五个样值搬到最后，并倒序依次向后搬，结果原记录上的最后四个样值冲掉，把最前面的四个单元充零。对于不满一个采样间隔的校正量则由插值滤波实现校正。由于对于同一道