专利名称：一种地球物理勘探地震资料处理方法
技术领域：
本发明涉及石油地球物理勘探地震资料处理技术，是一种二维地震资料连片处理中的速度文件拼接方法，可有效缩短地震资料处理周期，提高地震资料处理效率。
背景技术：
为整体认识盆地或宏观构造的构造、沉积特征、勘探潜力，进而优选勘探部署方向，需要对以往采集的二维地震资料进行拼接处理，进而进行构造解释、层位描述等研究工作。连片拼接处理可节约巨大的耗资，避免重新进行地震资料的采集。二维地震资料连片处理中，地震资料的拼接处理研究较多，在分块地震资料的预处理后，如静校正、噪音压制、振幅补偿、反褶积、速度分析等，再进行子波整形和面元均化处理，即可实现叠加地震资料的连片拼接，在进行偏移时，速度文件是必须的，传统方法是采用在叠前将地震资料进行拼接后重新抽道集，进而形成整个拼接线的速度谱文件，在重新提取速度谱后，可将拾取的速度文件作为偏移的速度输入，进而实现连片资料的叠后偏移成像处理。该类方法需要重新提取速度谱，增加了很大工作量的同时，使资料处理周期加长，影响勘探部署和勘探节奏。地震资料处理中的成像精度和处理效率是地球物理工作追求的两大目标，提高地震成像精度，是当前地球物理研究者的重要研究目标，在保证成像精度的情况下，提高地震资料的处理效率，同样具有重要意义。本发明可对二维地震资料中的各区域速度直接进行高效率拼接处理，无需重新拾取速度谱，可有效缩短地震资料处理周期，本发明给出的速度拼接方法是当前二维地震资料连片处理中的强有力工具之一。
发明内容
针对传统二维地震资料连片处理中的叠后偏移成像处理需要对连片后的地震数据重新拾取速度的不足之处，本发明提供了一种二维地震资料连片处理中的速度文件拼接方法。本发明通过以下技术方案实现一种地球物理勘探地震资料处理方法，具体为二维地震资料连片处理中的速度文件拼接方法，其采用以下步骤I)地震资料获取，应用检波器记录人工震源激发的反射地震波，并记录到磁带上；同时记录检波器和震源的空间大地坐标、地表高程、低速带和降速带厚度及速度参数，并记录到文本文件上；2)对每块地震数据，进行叠前预处理；对上面得到的检波器和震源的空间大地坐标、地表高程、低速带和降速带厚度及速度参数的反射地震数据进行静校正、噪音压制、振幅补偿、反褶积处理；在静校正处理中应用炮集折射波初至拾取方法得到的炮点和检波点静校正量，实现静校正处理；根据研究区地震资料的特点确定具体噪音压制方法，振幅补偿采取地表一致性振幅补偿方法，反褶积处理采用地表一致性反褶积方法；3)分块地震资料的速度分析
对上述叠前预处理后的地震数据，按CMP号自小到大进行数据重排，进而按CMP间隔(10、20或40个)形成相干速度谱道集，按不同时刻能量最大值为依据，对每一个速度谱道集按时间方向自上而下拾取时间速度对，对一条二维测线按CMP号自小而大拾取完成后，按文本格式将时间-速度对的值输出，形成地震速度文件，完成二维地震资料的速度分析。4)速度文件拼接处理对不同区块的地震速度文件，按连片处理中地震资料大地坐标自小而大规则进行速度文件中CMP号的修改，而同一空间位置的时间-速度对在修改前后的值保持不变；如不同区块间存在CMP间距不同的情况，需要进行地震资料的面元均化处理，将参加拼接处理的所有地震区块的资料的面元统一(即CMP距相同)；然后按照不同区块的拼接顺序，以第一区块的文本格式的速度文件为标准，CMP号修改后的其它时间-速度对文本文件依次复
制到第一区块的文件中，形成一个新的文本文件，该文件就是拼接后的速度文件。5)拼接后的速度文件能够用于不同工区二维连片地震资料的动校正处理、叠加处理和叠后偏移处理，得到连片的拼接处理最终成果，为宏观区域构造解释和地质认识提供基础资料。其中，对于噪音类型为面波的地震资料，所采取的相应的压制方法为自适应面波压制方法；对于噪音类型为随机噪音的地震资料，所采取的相应的压制方法为中值滤波法压制随机噪音方法；对于噪音类型为高能噪音的地震资料，所采取的相应的压制方法为区域异常滤波法压制高能噪音方法。优选地，步骤2)中噪音采用自适应面波衰减方法压制低频面波，压制的结果sjt)可表不为sr (t) = s ⑴· f (t)其中s(t)为含面波的地震记录，f(t)为滤波函数；/(O = 1( 7 > O)，
F( η)为一单调递增函数，n = s (tj) -1A (tj) I, s (tj)为tj时刻的振幅值，A (tj)为动校正后的道集上后一半偏移距的振幅的平均值。当n <0时候，面波压制的滤波函数f(t)= LO0步骤2)中采用中值滤波法压制随机干扰，中值滤波法的基本原理是在地震数据体中滤波后的一个点可用其邻域的点来表示，设原来的函数用f(x，y)表示，滤波后的函数g(x，y)可表述为g(u) = ^^(/G，>0)，w为滤波计算的区域。也就是说，滤波后的点是邻域里的中值。优选地，振幅补偿采取地表一致性振幅补偿方法，反褶积处理采用地表一致性反褶积方法；包括地表一致性振幅补偿方法和地表一致性反褶积方法的一致性处理是地震资料处理的基础工作，设在i点激发、在j点接收的地震记录为Xij (t)，该记录可分解为和炮点、检波点、中心点和偏移距点的褶积；
权利要求
1.一种地球物理勘探地震资料处理方法，具体为二维地震资料连片处理中的速度文件拼接方法，其采用以下步骤1)地震资料获取，应用检波器记录人工震源激发的反射地震波，并记录到磁带上；同时记录检波器和震源的空间大地坐标、地表高程、低速带和降速带厚度及速度参数，并记录到文本文件上；2)对每块地震数据，进行叠前预处理；对上面得到的检波器和震源的空间大地坐标、 地表高程、低速带和降速带厚度及速度参数的反射地震数据进行静校正、噪音压制、振幅补偿、反褶积处理；在静校正处理中应用炮集折射波初至拾取方法得到的炮点和检波点静校正量，实现静校正处理；根据研究区地震资料的特点确定具体噪音压制方法，振幅补偿采取地表一致性振幅补偿方法，反褶积处理采用地表一致性反褶积方法；3)分块地震资料的速度分析对上述叠前预处理后的地震数据，按CMP号自小到大进行数据重排，进而按CMP间隔 (10,20或40个)形成相干速度谱道集，按不同时刻能量最大值为依据，对每一个速度谱道集按时间方向自上而下拾取时间速度对，对一条二维测线按CMP号自小而大拾取完成后， 按文本格式将时间-速度对的值输出，形成地震速度文件，完成二维地震资料的速度分析。4)速度文件拼接处理对不同区块的地震速度文件，按连片处理中地震资料大地坐标自小而大规则进行速度文件中CMP号的修改，而同一空间位置的时间-速度对在修改前后的值保持不变；如不同区块间存在CMP间距不同的情况，需要进行地震资料的面元均化处理，将参加拼接处理的所有地震区块的资料的面元统一(即CMP距相同)；然后按照不同区块的拼接顺序，以第一区块的文本格式的速度文件为标准，CMP号修改后的其它时间-速度对文本文件依次复制到第一区块的文件中，形成一个新的文本文件，该文件就是拼接后的速度文件。5)拼接后的速度文件能够用于不同工区二维连片地震资料的动校正处理、叠加处理和叠后偏移处理，得到连片的拼接处理最终成果，为宏观区域构造解释和地质认识提供基础资料。
2.根据权利要求1所述的地球物理勘探地震资料处理方法，其特征在于，对于噪音类型为面波的地震资料，所采取的相应的压制方法为自适应面波压制方法；对于噪音类型为随机噪音的地震资料，所采取的相应的压制方法为中值滤波法压制随机噪音方法；对于噪音类型为高能噪音的地震资料，所采取的相应的压制方法为区域异常滤波法压制高能噪音方法。
3.根据权利要求2所述的地球物理勘探地震资料处理方法，其特征在于，所述步骤2) 中噪音采用自适应面波衰减方法压制低频面波，压制的结果sjt)可表示为Sr (t) = s(t) · f (t)其中S(t)为含面波的地震记录，f(t)为滤波函数;/(0=.、(7>0)，F(n)为1.0 + F(77)一单调递增函数，η = S (tj)1-1A (tj) I, S (tj)为tj时刻的振幅值,A(tj)为动校正后的道集上后一半偏移距的振幅的平均值。当η < O时候，面波压制的滤波函数f(t) = 1.0。
4.根据权利要求2所述的地球物理勘探地震资料处理方法，其特征在于，所述步骤2) 中采用中值滤波法压制随机干扰，中值滤波法的基本原理是在地震数据体中滤波后的一个点可用其邻域的点来表示，设原来的函数用f(x，y)表示，滤波后的函数g(x，y)可表述为
5.根据权利要求1所述的地球物理勘探地震资料处理方法，其特征在于，所述振幅补偿采取地表一致性振幅补偿方法，反褶积处理采用地表一致性反褶积方法；包括地表一致性振幅补偿方法和地表一致性反褶积方法的一致性处理是地震资料处理的基础工作，设在i点激发、在j点接收的地震记录为Xij(t)，该记录可分解为和炮点、检波点、中心点和偏移距点的褶积；
6.根据权利要求1所述的地球物理勘探地震资料处理方法，其特征在于所述步骤4)的面元均化处理在连片的二维地震资料存在面元不均时应用，经过面元均化处理后，可得到相同CMP距的地震资料。
7.根据权利要求1所述的地球物理勘探地震资料处理方法，其特征在于所述步骤4)的速度拼接的具体公式可表示为
8.根据权利要求1所述的地球物理勘探地震资料处理方法，其特征在于，所述步骤4) 的面元均化处理前后的参考位置的确定方法如下面元均化处理前参考CMP位置(⑶Prel)和面元均化处理后参考CMP位置(⑶Pm2)的确定方法如下对拼接的任意两个速度文件，两者在边界处出现重叠现象，按照空间坐标的按自小到大的方向进行拼接处理，这时依据地震资料拼接时面元均化处理前后该空间位置处所对应的CMP号分别记为⑶Pm1和⑶Pm2。
9.根据权利要求8所述的地球物理勘探地震资料处理方法，其特征在于，对于任意地震资料，如果其面元大小相等，则CDPm1等于CDPm2。
全文摘要
本发明公开一种地球物理勘探地震资料处理方法，具体为二维地震资料连片处理中的速度拼接方法。该二维地震资料连片处理中的速度文件拼接方法包括地震资料获取、对每块地震数据进行叠前预处理、分块地震资料的速度分析及速度文件拼接处理、拼接后的速度文件能够用于不同工区二维连片地震资料的动校正处理、叠加处理和叠后偏移处理。使用该方法依据不同地震资料的拼接空间位置、不同面元均化前后的CMP距、拼接点所对应的CMP号等参数，可实现二维速度资料的连片拼接，拼接后的速度文件和前述的地震资料作为输入，可实现拼接地震资料的叠后偏移。该技术可大大缩短二维地震资料拼接处理的周期，可为地下宏观构造的整体认识和勘探部署提供技术支撑。
文档编号G01V1/28GK102998704SQ201210536609
公开日2013年3月27日 申请日期2012年12月13日 优先权日2012年12月13日
发明者石颖, 王维红, 孙丽艳, 刘淑芬, 孙永河, 高宁 申请人:石颖