一种震源横波勘探的近地表静校正处理方法
【专利摘要】本发明提供了一种震源横波勘探的近地表静校正处理方法，属于地球物理勘探地震资料处理领域。本方法包括：(1)数据输入：输入加载了空间属性的地震勘探原始数据；(2)利用横波层析静校正法获得横波静校正量；利用比例系数法获得横波静校正量；(3)判断优化获得最优横波静校正量。本发明解决了震源横波的近地表静校正处理问题，将其结果应用于实际地震资料处理，获得了较高信噪比的横波剖面。
【专利说明】一种震源横波勘探的近地表静校正处理方法
【技术领域】
[0001]本发明属于地球物理勘探地震资料处理领域，具体涉及一种震源横波勘探的近地表静校正处理方法，可应用于横波勘探资料处理以及特殊条件下炸药震源激发三分量接收的震源横波地震资料的处理。
【背景技术】
[0002]由于横波勘探的成本高，生产应用不普遍，多数是利用转换波进行多波勘探。在某些特定的地震地质条件下，在近地表各向异性介质中炸药震源激发，多分量接收的多波勘探中，转换波不如震源横波强，即需要进行震源横波的处理。由于横波速度较纵波速度低，纵横向变化大，震源横波的近地表静校正量往往变化剧烈，且横波的初至时间难以准确拾取，可见解决震源横波静校正的难度更大，较纵波静校正问题更为突出。
[0003]当前多波勘探的资料处理中，近地表静校正处理技术还不是很成熟，多波资料的近地表静校正处理问题仍然是有待研究解决。
【发明内容】
[0004]本发明的目的在于解决震源横波的近地表静校正量难题，提供一种震源横波勘探的近地表静校正处理方法，改善震源横波剖面品质，进而提高处理成果解释的可靠性。
[0005]本发明是通过以下技术方案实现的:
[0006]一种震源横波勘探的近地表静校正处理方法，包括:
[0007](I)数据输入:输入加载了空间属性的地震勘探原始数据；
[0008](2)利用横波层析静校正法获得横波静校正量；利用比例系数法获得横波静校正量;
[0009](3)判断优化获得最优横波静校正量。
[0010]在步骤(3)之后进一步包括:
[0011](4)静校正应用:将最优横波静校正量应用于地震数据上；
[0012](5)数据输出:输出应用了最优横波静校正量后的地震数据。
[0013]所述横波层析静校正法包括:
[0014](Al)对原始地震数据进行预处理，得到改善的初至横波；所述预处理包括道编辑、极性校正、优势频带滤波和偏振能量归位；
[0015](A2)采用自动初至拾取与人工交互修正相结合的办法，进行横波的初至拾取；
[0016](A3)利用拾取的横波初至，进行横波层析静校正计算，获得横波低速带的等效速度及横波静校正量。
[0017]所述比例系数法包括:
[0018](BI)首先完成纵波分量的初至拾取、层析静校正计算，得到纵波低速带的等效速度及纵波静校正量；
[0019](B2)利用纵、横波微测井资料，根据处理基准面和替换速度，分别求取浅层的纵波静校正量和横波静校正量，然后将纵波静校正量除以横波静校正量的值作为比例系数；
[0020]所述纵、横波微测井资料包括层速度、厚度和地表高程信息；
[0021]所述处理基准面和替换速度是由用户给定的；
[0022](B3)利用所述比例系数求取横波近地表静校正量:将比例系数乘上纵波静校正量获得横波的层析静校正量。
[0023]所述步骤(B2)中，如果纵波层析静校正和横波层析静校正结果均可靠，则将步骤(BI)获得的纵波低速带的等效速度除以步骤(A3)获得的横波低速带的等效速度的值作为比例系数。
[0024]所述步骤(3)具体如下:
[0025]将利用比例系数法和利用横波层析静校正法获得的横波静校正量分别应用于地震数据上，然后对两个应用效果进行比较，如果利用比例系数法得到的横波静校正量的应用效果好于利用横波层析静校正法得到的横波静校正量的应用效果，则将利用比例系数法得到的横波静校正量作为最优横波静校正量，反之则将利用横波层析静校正法得到的横波静校正量作为最优横波静校正量；
[0026]如果两种方法得到的横波静校正量的应用效果都不好，则利用择优综合利用法获得最优横波静校正量；
[0027]所述应用效果是指构造形态和信噪比。
[0028]所述择优综合利用法包括:
[0029]步骤①:炮点采用比例系数法得到炮点静校正量，检波点分别采用比例系数法和横波层析静校正法，得到2套初叠加剖面；
[0030]步骤②:比较步骤①得到的2套初叠加剖面的优劣，找出构造形态合理且信噪比高的地段所对应的检波点号的范围；
[0031]步骤⑧:输出步骤②找出的那部分检波点的静校正量；
[0032]步骤④:检波点采用步骤⑧的输出静校正量，炮点分别采用比例系数法和横波层析静校正法，得到2套初叠加剖面；
[0033]步骤⑤:比较步骤④得到的2套剖面优劣，找出构造形态合理且信噪比高对应的炮点号范围；
[0034]步骤⑥:输出步骤⑤找出的那部分炮点的静校正量；
[0035]步骤⑦:同时调用步骤③和步骤⑥的输出，两者构成最优横波静校正量。
[0036]与现有技术相比，本发明的有益效果是:该发明在一定程度上(取决于资料条件)解决了震源横波的近地表静校正处理问题，应用于实际地震资料处理，获得了较高信噪比的横波剖面，且横波剖面品质远比该区的转换波剖面好，横波剖面与纵波剖面面貌整体一致性较好，反射波组特征明显，可用于地质解释。
【专利附图】
【附图说明】
[0037]图1为原始X分量横波单炮记录。
[0038]图2为原始X分量横波单炮初至拾取记录。
[0039]图3-1为采用“比例系数法”得到的炮点静校正量。
[0040]图3-2为采用“纵波层析静校正法”得到的炮点静校正量。
[0041]图3-3为采用“横波层析静校正法”得到的炮点静校正量。
[0042]图4-1为“比例系数法”的检波点静校正量。
[0043]图4-2为“纵波层析静校正法”的检波点静校正量。
[0044]图4-3为“横波层析静校正法”的检波点静校正量。
[0045]图5为“横波层析静校正法”的应用效果，与图6对比可以看出:方框所示部分是“横波层析静校正法”的效果好。
[0046]图6为“比例系数法”的应用效果，与图5对比可以看出:方框所示部分是“比例系数法”的效果好。
[0047]对比图5与图6可知:两种方法各有利弊，应综合利用。
[0048]图7a为对震源横波采用“比例系数法”得到的单炮记录。
[0049]图7b为对震源横波采用“横波层析静校正法”得到的单炮记录。
[0050]图7c为对震源横波做常量时移静校正的单炮记录。
[0051]图8为本发明方法的步骤框图。
[0052]图9a为震源横波未做静校正处理的单炮记录。
[0053]图9b为震源横波采用“比例系数法”静校正后的单炮记录。
【具体实施方式】
[0054]下面结合附图对本发明作进一步详细描述:
[0055]如图8所示,本发明横波近地表静校正处理方法包括:
[0056]根据处理资料是否具备良好品质的横波初至、纵横波近地表结构调查的具体情况，决定采用下列方法之一:⑴横波层析静校正方法；⑵比例系数法；⑶择优综合利用法。三种静校正方法分别阐述如下:
[0057](I)横波层析静校正法
[0058]偏振能量归位及波场净化横波初至拾取，横波层析静校正计算。利用震源横波初至，求取横波层析静校正量。在原始横波初至的品质较高(主要是指原始横波初至的信噪比高、能量强，易于软件自动拾取时间。能否获得准确的横波初至时间就是衡量标准)的情况下，可在震源横波记录上精确拾取横波初至，采用层析静校正技术建立表层横波低速带模型，计算震源横波静校正量。具体步骤如下:
[0059](I)采用商业软件(如可采用两方奇科公司的OMEGA处理系统，比如TRACE_SELECT、TOOLKIT、TV_FILTER等)完成预处理(包括:道编辑、极性校正、优势频带滤波、偏振能量归位)，得到改善的初至横波。
[0060](2)采用商业软件中的“自动初至拾取”与“人工交互修正”(是软件中的两个功能。)相结合的办法，进行横波的初至拾取。
[0061](3)与纵波层析静校正类同，利用拾取的横波初至，采用商业软件进行横波层析静校正计算，获得横波低速带的速度信息及横波静校正量。
[0062](2)比例系数法
[0063]纵波层析静校正和微测井表层结构调查，求取横波近地表静校正量。根据浅层纵、横波静的速度比例关系，求取横波近地表静校正量。具体步骤如下:
[0064](I)首先要认真细致地解决好纵波静校正问题，完成纵波分量的初至拾取、层析静校正计算，得到纵波低速带的速度信息及纵波静校正量(利用商业处理系统实现)。
[0065](2)利用纵、横波微测井资料(包括:层速度、厚度、地表高程信息)，根据处理基准面和替换速度(一般由用户给定)，分别求取浅层的纵、横波静校正量，进而求出比例系数；另外，在具有良好的纵、横波层析静校正结果的条件情况下，也可利用纵、横波层析静校正输出的纵、横波低速带的等效速度信息，求出比例系数。
[0066](3)利用纵、横波静校正量比例系数，求取横波近地表静校正量。即将比例系数乘上纵波的层析静校正量，即可获得横波的层析静校正量。
[0067](3)择优综合利用法
[0068]当静校正方法(I)和(2)的效果都不够理想时，且各有利弊时，可考虑采用该静校正方法。即在对比前面两种静校正效果的基础上，择优、综合利用横波静校正量，下面是一种简单的实施方式:
[0069]步骤①:炮点采用比例系数法的炮点静校正量，检波点另分别采用比例系数法、横波层析静校正法，得到2套初叠加剖面(通过采用处理系统里的常规的CDP道集分选、速度分析和叠加处理即可)；
[0070]步骤②:比较上述2套剖面优劣，找出构造形态合理(通过与纵波剖面中的构造形态对比可知)且信噪比(两者比较波组连续的信噪比高)高的地段所对应的检波点号范围；
[0071]步骤③:从2套检波点静校正量中择优输出一套构造形态合理且信噪比高对应的检波点静校正量；
[0072]步骤④:检波点采用⑧的输出静校正量，炮点另分别采用比例系数法、横波层析静校正法，得到2套初叠加剖面；
[0073]步骤⑤:比较2套剖面优劣，找出构造形态合理且信噪比高对应的炮点号范围；
[0074]步骤⑥:从2套炮点静校正量中择优输出一套构造形态合理且信噪比高对应的炮点静校正量；
[0075]步骤⑦:同时调用③和⑥的输出(即横波静校正量，它由炮点静校正量和检波点静校正量两部分组成)进行后续的处理。
[0076]横波近地表静校正是决定资料处理剖面品质的基础，只有在横波近地表静校正和地表一致性剩余静校正问题解决后，去噪问题才能得到较好的解决，进而才能提高速度精度和剖面品质。
[0077]本发明可较好地克服横波近地表静校正技术难题，为充分发挥震源横波在特殊地区地震勘探中的作用奠定了基础，可改善震源横波剖面品质，进而提高处理成果解释的可靠性。另外，该发明也可间接的应用到转换波的近地表静校正处理中，对检波点进行横波静校正。因此，该发明在多波勘探的处理中具有很奸的应用前景。
[0078]本发明的一个实施例如下:
[0079]根据工区内已知的纵、横波的测井资料，求得了比例系数为2，故图3-1和图4-1比例系数法的横波静校正量是通过纵波层析静校正量乘以2获得。图3-1与图3-3为不同近地表静校正方法的炮点静校正量对比，图4-1与图4-3为不同近地表静校正方法的检波点静校正量对比，可见它们之间的相似性较好，说明这两种方法有一定的可靠性。图7a与图7b对比，可见两种静校正方法都取得了明显效果，较没有进行横波静校正(图7c)的波组连续。
[0080]图5为“横波层析静校正法”的应用效果，与图6对比可以看出:方框所示部分是“横波层析静校正法”的效果好(主要看反射波组的形态是否合理和波组是否连续)。图6为采用“比例系数法”的应用效果，与图5对比可以看出:方框所示部分是“比例系数法”的效果好。对比图5与图6可知:两种方法各有利弊，应采用“择优综合利用法”。
[0081]图9a和图9b为一炮近地表静校正前后单炮记录的对比，其中，图9a为震源横波未做静校正处理的单炮记录，图％为震源横波采用“比例系数法”静校正后的单炮记录，比较两图可见静校正后的波组连续性改善明显。
[0082]本发明通过对多种近地表静校正处理方法的研究，设计了震源横波资料的近地表静校正方法，具体有:(I)横波层析静校正法；(2)比例系数法；(3)择优综合利用法。通过该发明应用，可以提高横波资料的品质和可解释性，充分发挥横波在油气勘探中的作用。
[0083]本发明在具备良好品质的横波初至或横波近地表结构调查资料条件的情况下，均可较好的解决横波的近地表静校正处理问题，但受原始记录中横波初至的能量和品质的限制，在难以获得理想的横波初至情况下，横波层析静校正方法的应用将受到限制，横波层析静校正的作用不能得到充分发挥，时好时坏，将来通过技术进步，在原始横波初至拾取精度能得到改进的情况下，利用大炮横波初至进行层析静校正，可获得更好的横波低速带模型和近地表静校正量；另受纵横波的微测井等资料的限制，在没有横波近地表结构调查资料时，就不能利用纵横波静校正量的比例关系，求取横波近地表静校正量；为此，在条件许可(即:既有好品质的横波初至，义有横波近地表结构调查资料)情况下，可通过综合利用两种静校正方法(横波层析静校正法和比例系数法)，获得更为可靠的结果，以适应复杂的静校正问题。
[0084]上述技术方案只是本发明的一种实施方式，对于本领域内的技术人员而言，在本发明公开了应用方法和原理的基础上，很容易做出各种类型的改进或变形，而不仅限于本发明上述【具体实施方式】所描述的方法，因此前面描述的方式只是优选的，而并不具有限制性的意义。
【权利要求】
1.一种震源横波勘探的近地表静校正处理方法，其特征在于:所述方法包括以下步骤: (1)数据输入:输入加载了空间属性的地震勘探原始数据； (2)利用横波层析静校正法获得横波静校正量；利用比例系数法获得横波静校正量； (3)判断优化获得最优横波静校正量。
2.根据权利要求1所述的震源横波勘探的近地表静校正处理方法，其特征在于:在步骤(3)之后进一步包括: (4)静校正应用:将最优横波静校正量应用于地震数据上； (5)数据输出:输出应用了最优横波静校正量后的地震数据。
3.根据权利要求1所述的震源横波勘探的近地表静校正处理方法，其特征在于:所述横波层析静校正法包括: (Al)对原始地震数据进行预处理，得到改善的初至横波；所述预处理包括道编辑、极性校正、优势频带滤波和偏振能量归位； (A2)采用自动初至拾取与人工交互修正相结合的办法，进行横波的初至拾取； (A3)利用拾取的横波初至，进行横波层析静校正计算，获得横波低速带的等效速度及横波静校正量。
4.根据权利要求3所述的震源横波勘探的近地表静校正处理方法，其特征在于:所述比例系数法包括: (BI)首先完成纵波分量的初至拾取、层析静校正计算，得到纵波低速带的等效速度及纵波静校正量； (B2)利用纵、横波微测井资料，根据处理基准面和替换速度，分别求取浅层的纵波静校正量和横波静校正量，然后将纵波静校正量除以横波静校正量的值作为比例系数； 所述纵、横波微测井资料包括层速度、厚度和地表高程信息； 所述处理基准画和替换速度是由用户给定的； (B3)利用所述比例系数求取横波近地表静校正量:将比例系数乘上纵波静校正量获得横波的层析静校正量。
5.根据权利要求4所述的震源横波勘探的近地表静校正处理方法，其特征在于:所述步骤(B2)中，如果纵波层析静校正和横波层析静校正结果均可靠，则将步骤(BI)获得的纵波低速带的等效速度除以步骤(A3)获得的横波低速带的等效速度的值作为比例系数。
6.根据权利要求1所述的震源横波勘探的近地表静校正处理方法，其特征在于:所述步骤(3)具体如下: 将利用比例系数法和利用横波层析静校正法获得的横波静校正量分别应用于地震数据上，然后对两个应用效果进行比较，如果利用比例系数法得到的横波静校正量的应用效果好于利用横波层析静校正法得到的横波静校正量的应用效果，则将利用比例系数法得到的横波静校正量作为最优横波静校正量，反之则将利用横波层析静校正法得到的横波静校正量作为最优横波静校正量； 如果两种方法得到的横波静校正量的应用效果都不好，则利用择优综合利用法获得最优横波静校正量； 所述应用效果是指构造形态和信噪比。
7.根据权利要求6所述的震源横波勘探的近地表静校正处理方法，其特征在于:所述择优综合利用法包括: 步骤①:炮点采用比例系数法得到炮点静校正量，检波点分别采用比例系数法和横波层析静校正法，得到2套初叠加剖面； 步骤②:比较步骤①得到的2套初叠加剖面的优劣，找出构造形态合理且信噪比高的地段所对应的检波点号的范围； 步骤⑧:输出步骤②找出的那部分检波点的静校正量； 步骤④:检波点采用步骤⑧的输出静校正量，炮点分别采用比例系数法和横波层析静校正法，得到2套初叠加剖画； 步骤⑤:比较步骤④得到的2套剖面优劣，找出构造形态合理且信噪比高对应的炮点号范围； 步骤⑥:输出步骤⑤找出的那部分炮点的静校正量； 步骤⑦:同时调用步骤③和步骤⑥的输出，两者构成最优横波静校正量。
【文档编号】G01V1/36GK104133244SQ201310157790
【公开日】2014年11月5日 申请日期:2013年5月2日 优先权日:2013年5月2日
【发明者】邬达理, 谌艳春, 余波, 朱蓓蓓 申请人:中国石油化工股份有限公司, 中国石油化工股份有限公司石油物探技术研究院