纵横波微测井资料处理方法及装置

技术领域

本发明涉及石油地震勘探近地表地层参数调查技术领域，尤其涉及一种纵横波微测井资料处理方法及装置。

背景技术

岩性油气藏、复杂裂隙裂缝型油气藏是目前我国陆上油气勘探主要的对象。相对于以往的构造油气藏，这些复杂油气藏难以用常规纵波(P波)勘探方法解决。理论上，多分量地震勘探技术在解决复杂结构与各向异性上有独到的优势，因此得到越来越多的研究和应用。在三分量地震资料采集工作中常用微测井方法采集表层的纵、横波资料，为消除或减弱近地表对地震资料品质产生的不利影响而提供基础数据。其目的是进一步提高表层结构调查精度，更有力地指导地震资料采集施工以及为多分量地震资料处理提供更可靠的静校正量数据，最终达到提升地震资料质量的目的。

但现有的对微测井纵、横波资料进行处理和解释等环节存在一些问题：

以往都是纵、横波信息分开进行的，包括资料加载、班报定义、分层解释等，过程比较繁琐，解释结果常存在相互矛盾的地方，很难进行融合统一，从而导致纵横波微测井资料的处理结果不够准确，进而影响多分量地震勘探的效果。

发明内容

本发明实施例提供一种纵横波微测井资料处理方法，用以提高纵横波微测井资料处理结果的准确性，该方法包括：

对目标区域的横波微测井资料进行带通滤波后，对所述纵横波微测井资料进行初至拾取，得到纵波初至时间数据和横波初至时间数据；

利用线性拟合方法，对纵波初至时间数据进行分层处理，得到纵波速度分层；利用非线性拟合方法，对横波初至时间数据进行分层处理，得到横波速度分层；

根据纵波速度分层和横波速度分层，进行纵横波联合解释，得到纵波表层结构参数和横波表层结构参数；

根据纵波表层结构参数、横波表层结构参数、目标区域的SPS成果数据，建立目标区域的纵横波表层结构模型。

具体实施时，利用如下公式，根据横波初至时间数据，对横波时深关系进行分段非线性拟合，确定横波分层速度：

Vs＝a

其中，Vs代表横波分层速度；a

具体实施例中，根据纵波表层结构参数、横波表层结构参数、目标区域的SPS成果数据，建立目标区域的纵横波表层结构模型，包括：

在建立横波表层结构时，在纵向上，层间速度采用非线性插值函数形式的变速填充，在横向上，大层结构以纵波表层结构分层为基准，层间横波速度采用非线性插值。

在本发明一具体实施例中，为了提高目标区域的纵横波微测井资料的精度，对目标区域的横波微测井资料进行带通滤波前还包括：

将目标区域内的三分量检波器的接收道，按照接收道的道类型进行归类；

将同一类的接收道所接收的纵波和横波分类显示，在同一分量窗口下对所述纵波和横波进行干扰数据过滤，确定目标区域的纵横波微测井资料。

在本发明另一具体实施例中，还提供一种纵横波微测井资料处理方法，在建立目标区域的纵横波表层结构模型后，还包括：

根据目标区域的纵横波表层结构模型，确定纵横波表层属性信息；

根据目标区域的纵横波表层结构模型、纵横波表层属性信息和目标区域的表层岩性录井资料，得到目标区域的近地表结构模型，确定目标区域的用于地震勘探的激发井深。

在本发明又一具体实施例中，还提供一种纵横波微测井资料处理方法，在建立目标区域的纵横波表层结构模型后，还包括：

根据目标区域的纵横波表层结构模型，选择纵波和横波的基准面参数、纵波替换速度和横波替换速度；

利用纵波和横波的基准面参数、纵波替换速度和横波替换速度，进行纵波静校正量和横波静校正量的提取计算，确定目标区域的纵波静校正量和横波静校正量。

本发明实施例还提供一种纵横波微测井资料处理装置，用以提高纵横波微测井资料处理结果的准确性，该装置包括：

滤波及初至拾取模块，用于对目标区域的横波微测井资料进行带通滤波后，对所述纵横波微测井资料进行初至拾取，得到纵波初至时间数据和横波初至时间数据；

速度分层模块，用于利用线性拟合方法，对纵波初至时间数据进行处理，得到纵波速度分层；利用非线性拟合方法，对横波初至时间数据进行处理，得到横波速度分层；

联合解释模块，用于根据纵波速度分层和横波速度分层，进行纵横波联合解释，得到纵波表层结构参数和横波表层结构参数；

纵横波表层结构模型建立模块，用于根据纵波表层结构参数、横波表层结构参数、目标区域的SPS成果数据，建立目标区域的纵横波表层结构模型。

具体实施例中，速度分层模块，具体用于：

利用如下公式，根据横波初至时间数据，对横波时深关系进行分段非线性拟合，确定横波分层速度：

Vs＝a

其中，Vs代表横波分层速度；a

具体实施过程中，纵横波表层结构模型建立模块，具体用于：

在建立横波表层结构模型时，在纵向上，层间速度采用非线性插值函数形式的变速填充，在横向上，大层结构以纵波表层结构分层为基准，层间横波速度采用非线性插值。

在本发明一具体实施例中，为了提高目标区域的纵横波微测井资料的精度，还提供一种纵横波微测井资料处理装置，还包括：

资料优选模块，用于：

对目标区域的横波微测井资料进行带通滤波前，将目标区域内的三分量检波器的接收道，按照接收道的道类型进行归类；

将同一类的接收道所接收的纵波和横波分类显示，在同一分量窗口下对所述纵波和横波进行干扰数据过滤，确定目标区域的纵横波微测井资料。

在本发明另一具体实施例中，还提供一种纵横波微测井资料处理装置，还包括：

激发井深确定模块，用于：

根据目标区域的纵横波表层结构模型，确定纵横波表层属性信息；

根据目标区域的纵横波表层结构模型、纵横波表层属性信息和目标区域的表层岩性录井资料，得到目标区域的近地表结构模型，确定目标区域的用于地震勘探的激发井深。

在本发明又一具体实施例中，还提供一种纵横波微测井资料处理装置，还包括：

静校正量计算模块，用于：

根据目标区域的纵横波表层结构模型，选择纵波和横波的基准面参数、纵波替换速度和横波替换速度；

利用纵波和横波的基准面参数、纵波替换速度和横波替换速度，进行纵波静校正量和横波静校正量的提取计算，确定目标区域的纵波静校正量和横波静校正量。

本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述纵横波微测井资料处理方法。

本发明实施例也提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述纵横波微测井资料处理方法的计算机程序。

本发明实施例中，通过对目标区域的横波微测井资料进行带通滤波后，对纵横波微测井资料进行初至拾取，得到纵波初至时间数据和横波初至时间数据；利用线性拟合方法，对纵波初至时间数据进行处理，得到纵波速度分层；利用非线性拟合方法，对横波初至时间数据进行处理，得到横波速度分层；根据纵波速度分层和横波速度分层，进行纵横波联合解释，得到纵波表层结构参数和横波表层结构参数；根据纵波表层结构参数、横波表层结构参数、目标区域的SPS成果数据，建立目标区域的纵横波表层结构模型。通过对目标区域的横波微测井资料进行带通滤波，能够去除横波中的干扰波信息，提高横波初至拾取时的精度和效率；通过将纵波和横波微测井资料同步进行处理，得到纵波速度分层和横波分层，并根据纵波速度分层和横波速度分层，进行纵横波联合解释，相较于现有技术中纵波和横波各自解释的方法，能够消除纵波和横波解释结果之间相互矛盾的情况，使得解释结果更加准确；由于横波初至拾取时的精度和纵横波资料解释结果准确度的提高，能够建立更加准确地纵横波表层结构模型，从而提高纵横波微测井资料处理结果的准确性。

附图说明

为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例中纵横波微测井资料处理方法的示意图。

图2为本发明一具体实施例中纵横波微测井资料处理方法的示意图。

图3为本发明另一具体实施例中纵横波微测井资料处理方法的示意图。

图4为本发明又一具体实施例中纵横波微测井资料处理方法的示意图。

图5为本发明具体实例中所提供的纵、横波微测井联合调查技术流程示意图。

图6为本发明具体实例中设计的1701测线的微测井点位完成图。

图7为本发明具体实例中激发点与接收点垂向关系图。

图8为本发明具体实例中采集到的纵、横波微测井资料示例图。

图9为本发明具体实例中接收道类型定义窗口示意图。

图10为本发明具体实例中道类型剔除分选定义窗口示意图。

图11为本发明具体实例中横波工作道增强识别后的示意图。

图12为本发明具体实例中初至拾取及初至时间显示情况。

图13为本发明具体实例中纵、横波的联合解释处理成果图。

图14为本发明具体实例中采用非线性拟合速度分层方法解释的横波微测井成果的时距图。

图15为本发明具体实例中建立的纵横波表层结构模型示意图。

图16为本发明具体实例中1701测线中的某一段测线采用线性分层与非线性分层方法解释的横波微测井资料计算的静校正量成果对比示意图。

图17为本发明具体实例中应用计算得到的静校正量单炮效果对比图。

图18为本发明具体实例中2221点应用老方法存档解释成果图。

图19为本发明具体实例中2457点应用老方法存档解释成果图。

图20为本发明实施例中纵横波微测井资料处理装置的示意图。

图21为本发明一具体实施例中纵横波微测井资料处理装置的示意图。

图22为本发明另一具体实施例中纵横波微测井资料处理装置的示意图。

图23为本发明又一具体实施例中纵横波微测井资料处理装置的示意图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明实施例提供了一种纵横波微测井资料处理方法，用以提高纵横波微测井资料处理结果的准确性，如图1所示，该方法包括：

步骤101：对目标区域的横波微测井资料进行带通滤波后，对纵横波微测井资料进行初至拾取，得到纵波初至时间数据和横波初至时间数据；

步骤102：利用线性拟合方法，对纵波初至时间数据进行处理，得到纵波速度分层；利用非线性拟合方法，对横波初至时间数据进行处理，得到横波速度分层；

步骤103：根据纵波速度分层和横波速度分层，进行纵横波联合解释，得到纵波表层结构参数和横波表层结构参数；

步骤104：根据纵波表层结构参数、横波表层结构参数、目标区域的SPS成果数据，建立目标区域的纵横波表层结构模型。

由图1所示流程可以得知，本发明实施例中，通过对目标区域的横波微测井资料进行带通滤波后，对纵横波微测井资料进行初至拾取，得到纵波初至时间数据和横波初至时间数据；利用线性拟合方法，对纵波初至时间数据进行处理，得到纵波速度分层；利用非线性拟合方法，对横波初至时间数据进行处理，得到横波速度分层；根据纵波速度分层和横波速度分层，进行纵横波联合解释，得到纵波表层结构参数和横波表层结构参数；根据纵波表层结构参数、横波表层结构参数、目标区域的SPS成果数据，建立目标区域的纵横波表层结构模型。通过对目标区域的横波微测井资料进行带通滤波，能够去除横波中的干扰波信息，提高横波初至拾取时的精度和效率；通过将纵波和横波微测井资料同步进行处理，得到纵波速度分层和横波分层，并根据纵波速度分层和横波速度分层，进行纵横波联合解释，相较于现有技术中纵波和横波各自解释的方法，能够消除纵波和横波解释结果之间相互矛盾的情况，使得解释结果更加准确；由于横波初至拾取时的精度和纵横波资料解释结果准确度的提高，能够建立更加准确地纵横波表层结构模型，从而提高纵横波微测井资料处理结果的准确性。

具体实施时，首先获取目标区域的横波微测井资料，对目标区域的横波微测井资料进行带通滤波，对目标区域的横波微测井资料进行带通滤波后，对目标区域的纵横波微测井资料进行初至拾取，得到纵波初至时间数据和横波初至时间数据。

由于野外采集环境比较复杂，而且在三分量检波器对微测井资料的采集过程中，对三分量检波器的耦合要求比较高，在同等激发条件下，由于环境噪音、地表因素、施工因素的影响对采集资料造成的干扰，压制和干扰了横波。通过采用带通滤波的方法，能够去除掉干扰波信息，提高横波初至拾取的精度并提高工作效率。

具体实施例中，从目标区域的三分量检波器中直接接收的纵横波微测井资料，可能存在不正常数据或者错误数据等干扰数据，为了在纵横波微测井资料处理前，对三分量检波器中直接接收的纵横波微测井资料进行优选，如图2所示，本发明具体实施例还提供一种纵横波微测井资料处理方法，在图1的基础上，还包括：

步骤201：将目标区域内的三分量检波器的接收道，按照接收道的道类型进行归类；

步骤202：将同一类的接收道所接收的纵波和横波分类显示，在同一分量窗口下对纵波和横波进行干扰数据过滤，确定目标区域的纵横波微测井资料。

三分量检波器的接收道的道类型分为X分量、Y分量和Z分量，将同一类型的接收道接收到的纵波和横波分类显示，在同一分量窗口下显示，并进行优选，一方面便于数据优选，能够提高目标区域的横波微测井资料的质量；另一方面，同窗口显示还能够提高数据过滤的速度，提高工作效率。

对目标区域的纵横波微测井资料进行初至拾取，得到纵波初至时间数据和横波初至时间数据后，利用线性拟合方法，对纵波初至时间数据进行处理，得到纵波速度分层；利用非线性拟合方法，对横波初至时间数据进行处理，得到横波速度分层。

由于现有技术中，横波表层结构速度分层解释方法采用的是线性分层法，进行横波速度分层具有以下缺点：分层节点不清晰，分层随意性很大，具多解性，并且精度较低，不利于横向标定层位追踪，体现不出近地表沉积的物性变化和结构变化特征。

具体实施时，鉴于横波速度连续性特征，其拟合分层的方法直接关系到横波速度求取的合理性和精度，为此采用非线性拟合方法，非线性拟合方法是通过拟合时深关系来控制分层精度，利用如下公式，根据横波初至时间数据，对横波时深关系进行分段非线性拟合，确定横波分层速度：

Vs＝a

其中，Vs代表横波分层速度；a

得到纵波速度分层和横波速度分层后，根据纵波速度分层和横波速度分层，进行纵横波联合解释，得到纵波表层结构参数和横波表层结构参数。具体实施例中，通过在纵、横波微测井资料处理窗口下，分左(纵波初至)右(横波初至)窗口同时加载纵横波初至数据，改变了以往只能单独将纵波和横波分别加载及处理的模式，极大地提高了资料处理的精度和效率。

得到纵波表层结构参数和横波表层结构参数后，根据纵波表层结构参数、横波表层结构参数、目标区域的SPS成果数据，建立目标区域的纵横波表层结构模型。其中，SPS(Shell Processing Support)成果数据是地震队采集的地震勘探辅助数据，按接收点文件、炮点文件、关系文件和注释文件组织并存储的一种记录格式。具体实施时，在建立横波表层结构时，在纵向上层间速度采用非线性插值函数形式的变速填充，在横向上，大层结构以纵波表层结构分层为基准，层间横波速度采用非线性插值。而纵波表层结构与现有模型构建方式一致，但横波的层位划分要参照纵波结构的主要分层界面，即在横波结构的层位确定上要建立与纵波结构的分层之间的关联性，最终实现纵波结构和横波结构模型的联合建立。

建立目标区域的纵横波表层结构模型后，本具体实施例还提供一种纵横波微测井资料处理方法，如图3所示，在图1的基础上，还包括：

步骤301：根据目标区域的纵横波表层结构模型，确定纵横波表层属性信息；

步骤302：根据目标区域的纵横波表层结构模型、纵横波表层属性信息和目标区域的表层岩性录井资料，得到目标区域的近地表结构模型，确定目标区域的用于地震勘探的激发井深。

建立目标区域的纵横波表层结构模型后，另一具体实施例中还提供一种纵横波微测井资料处理方法，如图4所示，在图1的基础上，还包括：

步骤401：根据目标区域的纵横波表层结构模型，选择纵波和横波的基准面参数、纵波替换速度和横波替换速度；

步骤402：利用纵波和横波的基准面参数、纵波替换速度和横波替换速度，进行纵波静校正量和横波静校正量的提取计算，确定目标区域的纵波静校正量和横波静校正量。

步骤402具体实施时，最终的横波静校正量的计算公式为：

其中，Δt代表横波在某一个深度段的传播时间；

a代表截距；b代表横波速度变化梯度；

H代表深度；

Δt

i代表横波表层模型高速顶以上的层数；

h

vs代表横波替换速度，单位m/s；

T代表横波检波点静校正量，单位m/s。

具体实施例中，为了保证对目标区域进行合理细致的勘察，需要获取到的纵横波微测井资料能够全面且精细化，因而也就需要三分量检波器位置布置合理。因此，在纵横波微测井资料处理前，优选地需要对目标区域的表层结构进行初步调查并设计施工，具体包括地表资料采集和施工设计，根据采集的目标区域的地表资料，调整布置三分量检波器的点位密度。

同时，在三分量检波器工作时，为了保证资料采集的质量，具体实施例中，严格监控井中三分量检波器进行微测井纵、横波资料采集的实施情况，并对资料采集过程中的质量问题要及时处理，其中，三分量检波器工作采用的是同步激发地面三分量检波器同步接收的方式。

下面给出一具体实例说明本发明实施例如何进行纵横波微测井资料处理。本例应用于海拉尔盆地贝尔凹陷贝39井工区南侧的1701测线。该测线位于该工区的南侧，横向上地形起伏较大，表层结构比较复杂，该线共设计采集了18口三分量微测井作为实例分析，选取的桩号范围为2047-2668一段。

如图5所示，为本具体实例中所提供的纵、横波微测井联合调查技术流程图。主要包括如下步骤：

首先，对纵、横波微测井调查点位设计及调整：依据地形变化布置有控制点，且点位间距不能相差太大，还要考虑便于施工作业等原则，进行纵、横波微测井调查点位设计，对个别点位进行必要的调整，最终设计的1701测线的微测井点位完成图如图6所示。

在设计好的调查点位上进行纵、横波微测井资料采集，如图7所示为本具体实例中激发点与接收点垂向关系图，本具体实例将井筒深度设定为30m～80m(施工中可根据实际情况进行微测井的井深调整)，微测井激发电缆上的激发点密度由浅到深为0.5m～2m；采用单侧8只三分量检波器，排列方式为每1只检波器一组，呈线形排列，每只三分量检波器分别按着X、Z、Y的顺序进行连接，井检距分别为2m、4m、6m、8m、10m、12m、14m，16m，且X分量方向指向井口。还可以根据需要改变检波器接收图形，如扇形和L形，如图8所示，为采集到的纵、横波微测井资料示例图。

为了更好地选择纵、横波工作道，把每个接收道都按道类型(X、Y、Z分量)分别进行定义，这样可以对各分量的道进行归类，如图9所示，定义类型分为三类，即Z分量纵波接收道0、X分量横波接收道2、需要剔除的工作道1，把所有的纵、横波分类显示，在同一个分量窗口下进行优选，通过此设置可以极为方便地分选各分量。如图10所示，是道分选显示窗口，经过纵、横波微测井记录每一道的分量类型定义后，可以在显示波形资料的窗口选择需要显示的各分量道的记录，对纵、横波的集中优选提供了很大的便利，可以保证应用最优的工作道进行资料解释，以保证资料质量。

具体地，对观测系统定义：接收道数24道，200200200200200200200200，首先是定义接收道数，之后是定义偏移距即每个接收道与井口的平面距离。对道类型定义：210210210210210210210210，数字“0”代表纵波分量，数字“1”代表横波“Y”分量，数字“2”代表横波“X”分量。

选出最优的工作道后，对横波微测井资料进行带通滤波，以消除环境噪音、地表因素、施工因素的影响对采集资料造成的干扰。如图11所示，为本具体实例中横波工作道增强识别后的示意图，图的上半部分是未进行分频滤波的波列图，图的下半部分进行了分频滤波的波列图，对横波信息的识别以及初至拾取效果显著。

对横波微测井资料进行带通滤波后，进行纵、横波初至时间拾取：首先估算出一个完整横波波形的频率，保留估算值一定范围的有效频率，给定通放频率范围：0，10，50，60。

由于井中激发地面接收所得到的横波不是初至波，而是续至波，波形起跳的位置不好分辨，所以把波峰或者波谷作为拾取初至的标准位置。纵波初至拾取是把纵波起跳初至位置作为初至拾取的标准位置，同时可结合人工干预或编辑，完成初至拾取，例如，图12所示为初至拾取及初至时间显示情况。

具体实施时，完成初至拾取后，在同尺度显示窗口下，首先进行纵、横波预处理速度分层，纵波还是应用线性拟合方法，横波速度采用非线性函数拟合分层方法，然后通过各自的影响因素来完成纵、横波的联合解释，处理成果如图13所示，图中左侧纵波初至显示窗口，右侧为横波初至显示窗口，中间显示的是岩性柱。

鉴于横波速度连续性特征，其拟合分层的方法直接关系到横波速度求取的合理性和精度，为此采用非线性拟合方法，非线性拟合方法是通过拟合时深关系来控制分层精度，如图14所示，为采用非线性拟合速度分层方法解释的横波微测井成果的时距图，从图中可以清晰地看出分辨岩性介质变化规律。

进行纵横波联合解释，得到如图13所示的纵波时距图和横波时距图后，根据纵波表层结构参数、横波表层结构参数、目标区域的SPS成果数据，并加载表层岩性资料，建立纵横波表层结构模型。具体过程，包括：在建立横波表层结构时，在纵向上，层间速度采用非线性插值函数形式的变速填充，在横向上，大层结构以纵波表层结构分层为基准，层间横波速度采用非线性插值。而纵波表层结构与现有模型构建方式一致，但横波的层位划分要参照纵波结构的主要分层界面，即在横波结构的层位确定上要建立与纵波结构的分层之间的关联性，最终实现纵波结构和横波结构模型的联合建立。最终建立的纵横波表层结构模型如图15所示，图15中的上半部分显示的是所建立的纵波表层结构模型，图15中的下半部分显示的是所建立的横波表层结构模型。

建立纵横波表层结构模型后，选择基准面参数以及各自的替换速度来完成各自的静校正量的提取计算。

本具体实例中，这一条测线的纵、横波微测井点位设计、资料采集、资料加载、横波信息的增强识别、初至拾取、纵横波速度拟合分层、纵横波结构模型构建及静校正量提取等各环节的实际应用目的性明确、各环节衔接顺畅、反馈及时，提高了整个流程的效率、并明显提高了最终资料处理结果的精度。

为了验证本发明所提供的纵横波微测井资料处理方法所能达到的效果，本发明具体实施例中还在大庆探区利用上述具体实例中所提供的近地表结构微测井纵、横波联合同步调查技术，完成了三块三分量地震勘探三维采集工程项目，合计处理三分量微测井1237口，纵、横波静校正量各自提交物理点数是648063个，取得了非常理想的应用效果，具体包括：

(1)资料解释效率和精度的提高：通过采用横波速度拟合分层方法和新的纵、横波处理方法，工作效率提高4倍以上，具体数据详见表1，资料解释精度平均提高16.9％，具体数据详见表2。

表1老方法解释与新方法时效分析表

表2老方法解释与新方法解释精度分析表

(2)模型构建和静校正量计算的效率和精度的提高：在提高解释结果的精度和效率的基础上，又运用了联合构建纵、横波表层结构模型的模式，使得表层模型的精度提高至少16.9％以上，静校正量的精度提高也至少在16.9％以上，如图16和图17所示，图16为上述具体实例中1701测线中的某一段测线采用线性分层与非线性分层方法解释的横波微测井资料计算的静校正量成果对比，图的最上部分显示的是地表高程线，图的中间部分是两种分层方法计算的静校正量曲线，图的最下部分显示的是相应的误差分布，由于在计算横波静校正量时，都是将解释成果中最深的拐点作为高速顶界面，箭头所指的位置是差值最大的区域，该区域内的那口微测井的两种解释方法确定的高速顶界面时间相差近40ms，而本发明提供的新方法解释成果计算的静校正量曲线与地形变化相似，符合检波点静校正量的变化规律，由此看见本申请提供的新方法是比较可信的。图17中的a部分是某区块横波原始单炮，b部分是采用老处理方法计算的静校正量单炮，c部分是采用本发明提供的新处理方法计算的静校正量单炮，从图上不难看出采用新方法成果计算的静校正量单炮的同相轴连续性较好，应用效果提高较为明显。

从图16和图17不难看出，前期成果的有效性，相应的减少了结构模型构建的对比计算次数，效率上虽有提高但不明显。

(3)通过采用微测井纵、横波联合调查设计、采集、处理一体化技术流程，对资料的精度及有效性提高显著：

表3纵、横波联合调查处理技术应用效果综合分析表

通过两个工区的应用分析，具体测试结果统计数据见表3。近地表结构微测井纵、横波联合调查一体化技术流程已经成功解决了施工设计、资料采集、横波速度拟合分层方法问题和微测井纵、横波联合解释的效率和精度问题，通过微测井纵、横波联合处理和近地表微测井纵、横波结构模型联合解释向上级流程形成有效的反馈，更好的控制了单点资料质量和静校正量计算精度，在近地表结构微测井纵、横波调查各环节衔接、反馈顺畅，在整体环节上提高的效率比较显著。

(4)资料效果对比

根据纵、横波联合观测表层调查的应用技术流程，从计算精度和工作效率两方面阐述微测井纵、横波联合同步调查技术的应用效果，选取的实例资料是海拉尔盆地贝尔凹陷贝39井工区南侧的1701测线，该线横向上地形起伏较大，表层结构比较复杂，该线共设计采集了18口三分量微测井作为实例分析，测线位于该工区的南侧，横向上地形起伏较大，表层结构比较复杂，选取1701测线的桩号2047-2668一段，用新老方法分别计算静校正量，从图16上的两个静校正量曲线对比可以看出，在2147-2287区间，静校正量差异较大，从2301到最后，其变化形态与之间的差异较小，吻合度较好。为了分析产生差异的原因，在两个不同区段分别选取代表点位2221和2457。

下面对2221与2457两个点的新老处理方法的时距图和处理结果进行对比分析：

对于2221点，图18是2221点应用老方法存档解释成果图，资料分层非常简单，对横波速度变化规律描述不清，精度较低。交互解释资料高速顶位置较存档资料的高速顶位置低了近16米、40ms，差异非常明显。因此，从图18和表4对比分析可以看出，采用老处理方法的横波分层为两层，采用新处理方法的横波分层为4层，在计算横波静校正量时，都是将解释成果中最深的拐点作为高速顶界面，两种解释方法确定的高速顶界面时间相差近40ms，而新处理方法的成果计算的静校正量曲线与地形变化相似，符合检波点静校正量的变化规律，是比较可信的。

表4该物理点采用新处理技术前后数据对比

对于2457点，图19和表5给出了新老方法解释的时距图。该点新处理方法的横波分层为3层，而图19显示的老方法处理的横波分层为2层，分层非常简单，本申请提供的新方法较老方法的分层比较详细，处理精度高。但是最下面的高速层拐点时间仅相差不到2ms，而且两种方法确定的高速顶界面相差也小(新23.1m-老23.98m＝-0.88m)，所以两种方法计算的静校正量曲线相差不大，均与地形曲线变化相似，可信度均较高。

表5该物理点采用新处理技术前后数据对比

三分量表层结构调查的最终目的是构建精细且准确的横波表层结构模型。构建精细准确的横波表层结构模型的条件是：既要在纵向上能体现横波的结构分层及速度变化规律来，又要在横向上体现出速度变化规律和追踪连层，以保证纵、横波模型的层位的合理性、规律性及可追踪性，由此才能计算出精准的静校正量，而本发明所建立的横波表层结构模型正是基于以上思路所建立的，所以能够构建精细且准确的横波表层结构模型，进而计算出精准的静校正量。

上述具体应用的实施仅为举例，其余实施方式不再一一赘述。

基于同一发明构思，本发明实施例还提供一种纵横波微测井资料处理装置，由于纵横波微测井资料处理装置所解决问题的原理与纵横波微测井资料处理方法相似，因此纵横波微测井资料处理装置的实施可以参见纵横波微测井资料处理方法的实施，重复之处不再赘述，具体结构如图20所示：

滤波及初至拾取模块2001，用于对目标区域的横波微测井资料进行带通滤波后，对纵横波微测井资料进行初至拾取，得到纵波初至时间数据和横波初至时间数据；

速度分层模块2002，用于利用线性拟合方法，对纵波初至时间数据进行处理，得到纵波速度分层；利用非线性拟合方法，对横波初至时间数据进行处理，得到横波速度分层；

联合解释模块2003，用于根据纵波速度分层和横波速度分层，进行纵横波联合解释，得到纵波表层结构参数和横波表层结构参数；

纵横波表层结构模型建立模块2004，用于根据纵波表层结构参数、横波表层结构参数、目标区域的SPS成果数据，建立目标区域的纵横波表层结构模型。

具体实施例中，速度分层模块2002，具体用于：

利用如下公式，根据横波初至时间数据，对横波时深关系进行分段非线性拟合，确定横波分层速度：

Vs＝a

其中，Vs代表横波分层速度；a

具体实施时，纵横波表层结构模型建立模块2004，具体用于：

在建立横波表层结构时，在纵向上，层间速度采用非线性插值函数形式的变速填充，在横向上，大层结构以纵波表层结构分层为基准，层间横波速度采用非线性插值。

具体实施例中，还提供一种纵横波微测井资料处理装置，如图21所示，在图20的基础上，还包括：

资料优选模块2101，用于：

对目标区域的横波微测井资料进行带通滤波前，将目标区域内的三分量检波器的接收道，按照接收道的道类型进行归类；

将同一类的接收道所接收的纵波和横波分类显示，在同一分量窗口下对纵波和横波进行干扰数据过滤，确定目标区域的纵横波微测井资料。

具体实施时，如图22所示的纵横波微测井资料处理装置，在图20的基础上，还包括：

激发井深确定模块2201，用于：

根据目标区域的纵横波表层结构模型，确定纵横波表层属性信息；

根据目标区域的纵横波表层结构模型、纵横波表层属性信息和目标区域的表层岩性录井资料，得到目标区域的近地表结构模型，确定目标区域的用于地震勘探的激发井深。

具体实施例中，如图23所示为另一具体实施例中所提供的纵横波微测井资料处理装置，在图20的基础上，还包括：

静校正量计算模块2301，用于：

根据目标区域的纵横波表层结构模型，选择纵波和横波的基准面参数、纵波替换速度和横波替换速度；

利用纵波和横波的基准面参数、纵波替换速度和横波替换速度，进行纵波静校正量和横波静校正量的提取计算，确定目标区域的纵波静校正量和横波静校正量。

本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述纵横波微测井资料处理方法。

本发明实施例还提供一种计算机可读存储介质，存储有执行上述纵横波微测井资料处理方法的计算机程序。

综上所述，本发明实施例提供的纵横波微测井资料处理方法及装置具有如下优点：

通过对目标区域的横波微测井资料进行带通滤波后，对纵横波微测井资料进行初至拾取，得到纵波初至时间数据和横波初至时间数据；利用线性拟合方法，对纵波初至时间数据进行处理，得到纵波速度分层；利用非线性拟合方法，对横波初至时间数据进行处理，得到横波速度分层；根据纵波速度分层和横波速度分层，进行纵横波联合解释，得到纵波表层结构参数和横波表层结构参数；根据纵波表层结构参数、横波表层结构参数、目标区域的SPS成果数据，建立目标区域的纵横波表层结构模型。通过对目标区域的横波微测井资料进行带通滤波，能够去除横波中的干扰波信息，提高横波初至拾取时的精度和效率；通过将纵波和横波微测井资料同步进行处理，得到纵波速度分层和横波分层，并根据纵波速度分层和横波速度分层，进行纵横波联合解释，相较于现有技术中纵波和横波各自解释的方法，能够消除纵波和横波解释结果之间相互矛盾的情况，使得解释结果更加准确；由于横波初至拾取时的精度和纵横波资料解释结果准确度的提高，能够建立更加准确地纵横波表层结构模型，从而提高纵横波微测井资料处理结果的准确性。

本领域内的技术人员应明白，本发明的实施例可提供为方法、装置、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

本发明是参照根据本发明实施例的方法、装置和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。