声波测井校正方法、装置、设备、存储介质及程序产品与流程

1.本发明涉及石油勘探技术领域，尤其涉及一种声波测井校正方法、装置、设备、存储介质及程序产品。


背景技术：

2.目前，在油气勘探领域，受仪器本身与地下复杂地质条件的共同影响，在裸眼井的测井过程中，声波时差的测量值与真实值往往存在一定的误差，从而导致基于声波测井的人工合成记录与其后续的井震标定工作受不同程度的影响，可能得到一些错误的结论。
3.现有的校正方法有井眼补偿法，这种校正方法主要校正由井眼塌垮导致的测量误差。这种方法具有一定的局限性，声波测井校正的效率和准确率还有待于提高。


技术实现要素：

4.本发明提供一种声波测井校正方法、装置、设备、存储介质及程序产品，以解决现有校正技术的一些局限性以及提高声波测井校正的准确率。
5.第一方面，本发明提供一种声波测井校正方法，所述方法包括：
6.获取目标区域的多个被测井的垂直地震剖面测井速度曲线以及声波速度曲线；
7.对于每一被测井，根据所述被测井的垂直地震剖面测井速度曲线以及声波速度曲线进行拟合操作，得到所述被测井对应的拟合参数；
8.根据多个被测井对应的拟合参数，得到所述目标区域的垂直地震剖面测井速度曲线与声波速度曲线的拟合参数；
9.根据所述目标区域对应的拟合参数，对所述目标区域除所述多个被测井以外的其它井的声波速度曲线进行校正。
10.可选的，根据所述被测井的垂直地震剖面测井速度曲线以及声波速度曲线进行拟合操作，得到所述被测井对应的拟合参数，包括：
11.建立垂直地震剖面测井速度曲线与声波速度曲线的回归函数，所述回归函数包括拟合参数；
12.在所述被测井的垂直地震剖面测井速度曲线与声波速度曲线上分别选取多个点，根据所选取的点对应的速度值进行拟合操作，为所述被测井对应的拟合参数赋值。
13.可选的，回归函数为
[0014][0015]
拟合操作的收敛条件为取最小值；
[0016]
其中，表示被测井对应的回归函数，a0，a1，
…
，an表示被测井的拟合参数，xi表示声波速度曲线上的第i个点对应的速度值，yi表示实测垂直地震剖面测井速度曲线上的第i个点对应的速度值，表示拟合公式计算值，ξ表示单井误差值，n表示为大于等于2
的正整数，p表示为选取的点数。
[0017]
可选的，根据多个被测井对应的拟合参数，得到所述目标区域的垂直地震剖面测井速度曲线与声波速度曲线的拟合参数，包括：
[0018]
将所述多个被测井对应的拟合参数代入以下公式：
[0019][0020]
其中，δ
l
表示第l组拟合参数下m个被测井的整体误差，表示第l组拟合参数对第j个被测井进行计算得到的单井误差，ξj表示第j个被测井取得最佳拟合参数的最小单井误差，m表示为被测井的个数；
[0021]
根据收敛条件，从所述多个被测井对应的拟合参数中选择所述目标区域对应的拟合参数；
[0022]
其中，所述收敛条件为δ
l
取最小值。
[0023]
可选的，对于每一被测井，根据所述被测井的垂直地震剖面测井速度曲线以及声波速度曲线进行拟合操作，得到所述被测井对应的拟合参数，包括：
[0024]
建立所述目标区域的垂直地震剖面测井速度曲线与声波速度曲线的交互图，所述交互图的横轴为声波测井速度，纵轴为垂直地震剖面测井速度，所述交互图中的交互点的位置由多个被测井的垂直地震剖面测井速度曲线以及声波速度曲线确定；
[0025]
根据所述交互图中的各个交互点的位置，确定所述各个交互点的集中程度是否满足预设条件；
[0026]
若满足，则对于每一被测井，根据所述被测井的垂直地震剖面测井速度曲线以及声波速度曲线进行拟合操作，得到所述被测井对应的拟合参数。
[0027]
可选的，对所述目标区域除所述多个被测井以外的其它井的声波速度曲线进行校正，包括：
[0028]
获取所述目标区域中待分析的井的声波速度曲线；
[0029]
根据所述待分析的井的声波速度曲线以及所述目标区域对应的拟合参数，确定所述待分析的井对应的拟垂直地震剖面测井速度曲线；
[0030]
根据所确定的拟垂直地震剖面测井速度曲线，对所述待分析的井进行井震标定。
[0031]
第二方面，本发明提供一种声波测井校正装置，包括：
[0032]
获取模块，用于获取目标区域的多个被测井的垂直地震剖面测井速度曲线以及声波速度曲线；
[0033]
第一得到模块，用于对于每一被测井，根据所述被测井的垂直地震剖面测井速度曲线以及声波速度曲线进行拟合操作，得到所述被测井对应的拟合参数；
[0034]
第二得到模块，用于根据多个被测井对应的拟合参数，得到所述目标区域的垂直地震剖面测井速度曲线与声波速度曲线的拟合参数；
[0035]
校正模块，用于根据所述目标区域对应的拟合参数，对所述目标区域除所述多个被测井以外的其它井的声波速度曲线进行校正。
[0036]
第三方面，本发明提供一种声波测井校正设备，包括：
[0037]
至少一个处理器和存储器；
[0038]
所述存储器存储计算机执行指令；
[0039]
所述至少一个处理器执行所述存储器存储的计算机执行指令，使得所述至少一个处理器执行如第一方面任一项所述的声波测井校正方法。
[0040]
第四方面，本发明提供一种计算机可读存储介质，所述计算机可读存储介质中存储有计算机执行指令，所述计算机执行指令被处理器执行时用于实现如第一方面任一项所述的声波测井校正方法。
[0041]
第五方面，本发明提供一种计算机程序产品，包括计算机程序，所述计算机程序被处理器执行时实现如第一方面任一项所述的声波测井校正方法。
[0042]
本发明提供了一种声波测井校正方法、装置、设备、存储介质和程序产品，所述方法包括：通过获取目标区域的多个被测井的垂直地震剖面测井速度曲线以及声波速度曲线，对于每一被测井，根据所述被测井的垂直地震剖面测井速度曲线以及声波速度曲线进行拟合操作，得到所述被测井对应的拟合参数，根据多个被测井对应的拟合参数，得到所述目标区域的垂直地震剖面测井速度曲线与声波速度曲线的拟合参数，根据所述目标区域对应的拟合参数，对所述目标区域除所述多个被测井以外的其它井的声波速度曲线进行校正，从而可以通过统计学手段分析隐含在目标区块的内在规律，根据多个被测井建立目标区域的校正方案，实现声波测井校正的定量描述，在一定程度上弥补了现有校正技术的局限性，大幅提升声波测井的标定精度与可靠性，提高校正效率，降低勘探成本。
附图说明
[0043]
为了更清楚地说明本发明或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
[0044]
图1为本发明实施例提供的一种应用场景示意图；
[0045]
图2为本发明实施例提供的一种声波测井校正方法的流程示意图；
[0046]
图3为本发明实施例提供的一种校正前的多井声波测井速度与垂直地震剖面测井速度交互图；
[0047]
图4为本发明实施例提供的一种校正后多井声波测井速度与垂直地震剖面测井速度交互图；
[0048]
图5a为本发明实施例提供的校正前多井声波测井速度与垂直地震剖面测井速度概率分布图；
[0049]
图5b为本发明实施例提供的校正后多井声波测井速度与垂直地震剖面测井速度概率分布图；
[0050]
图6a为本发明实施例提供的校正前声波速度曲线与垂直地震剖面速度曲线叠合图；
[0051]
图6b为本发明实施例提供的校正后声波速度曲线与垂直地震剖面速度曲线叠合图；
[0052]
图7a为本发明实施例提供的校正前声波测井与垂直地震剖面测井时深关系图；
[0053]
图7b为本发明实施例提供的校正后声波测井与垂直地震剖面测井时深关系图；
[0054]
图8为本发明实施例提供的校正前后某地区2井综合标定剖面的对比图；
[0055]
图9为本发明实施例提供的校正前后某地区1井综合标定剖面的对比图；
[0056]
图10为本发明实施例提供的声波速度校正图版；
[0057]
图11为本发明实施例提供的声波测井校正装置的结构示意图；
[0058]
图12为本发明实施例提供的声波测井校正设备的硬件结构示意图。
具体实施方式
[0059]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0060]
本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”、“第三”“第四”等(如果存在)是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本发明的实施例例如能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
[0061]
图1为本发明实施例提供的一种应用场景示意图，在对油田的某些复杂区块进行油气勘探时，当地震勘探确定了区域的地质构造后，需要对地层中被测井可能的含油情况进行一定的了解。通常是采用声波测井技术进行测井作业，直接测井得到的数据具有一定的误差，因此需要对测井结果进行校正。
[0062]
其中，在一些技术中，可以采用井眼补偿对声波测速进行校正，这种校正方法通常是校正由井眼垮塌导致的测量误差。具体地，先求取垮塌指数来描述井眼的垮塌程度，再根据垮塌指数的大小对声波测井进行一定程度的补偿。但是当垮塌严重时，这种方法不适用。
[0063]
在另一些技术中，采用泥浆侵入对声波时差进行校正，校正受钻进液侵入比较严重的地层。具体地，可利用电阻率或者其他测井曲线求取一个蚀变指数来描述钻井液侵入原状地层深度及影响程度，再根据蚀变指数对声波时差进行校正。但是这种方法比较局限，只考虑了泥浆侵入这一种影响因素。
[0064]
在一些综合技术中，可以将两种校正方法结合到一起，同时利用垮塌指数和蚀变指数对声波时差进行校正。但是，这些方法只是考虑了井眼垮塌和泥浆侵入对声波测井带来的影响，并将其定量化。因此，这些方法并不适应于一些复合盆地。
[0065]
在某些复杂区块，由于声波测井的误差较大，常规校正方法又无法满足需求，往往借助可靠性更高的垂直地震剖面测井来完成井震标定工作，这样一来，很大程度上增加了勘探的成本，为解决异常复杂地质条件下声波时差的校正问题、尽可能消除多种影响因素的对声波测井的综合影响，本发明实施例提供一种声波测井校正方法，将某一个已测得几个井垂直地震剖面测井数据的区域进行统计分析，以速度为切入点，通过将垂直地震剖面测井速度曲线与声波测井速度曲线二者之间进行拟合，拟合得到声波速度的校正公式，从而实现对该区域其他井的声波测井的校正，在一定程度上弥补了现有校正技术的局限性和
不足。
[0066]
下面以具体地实施例对本发明的技术方案进行详细说明。下面这几个具体的实施例可以相互结合，对于相同或相似的概念或过程可能在某些实施例不再赘述。
[0067]
图2为本发明实施例提供的一种声波测井校正方法的流程示意图，如图2所示，本实施例提供的方法，可以包括：
[0068]
步骤201、获取目标区域的多个被测井的垂直地震剖面测井速度曲线以及声波速度曲线。
[0069]
可选的，目标区域可以为待处理的多个被测井所在的区域，所述多个被测井所处的位置可以具有相同或相似的地质条件，横向上没有太大的跨度。具体地，目标区域可以包含处在同一个二级构造单元的区域。其中，二级构造单元为盆地中由一系列相似的单一构造所组成的构造带。
[0070]
在本实施例中，选取目标区域中多个井作为被测的井，根据垂直地震剖面(vetical seismic profile，简称vsp)测井和声波测井技术，得到多个井对应的两条速度曲线，一条是基于垂直地震剖面测井的速度曲线，一条是基于声波测井计算得出的声波速度曲线。
[0071]
其中，基于垂直地震剖面测井的速度曲线图与声波速度曲线图横轴均为井深，纵轴为速度。
[0072]
可选的，声波速度曲线是基于声波时差曲线计算得到的。声波时差速度与声波测井速度互为倒数。
[0073]
其中，声波时差曲线，又可称为声波曲线和声波测井曲线，是声波测井时由测得的声波时差速度得到的曲线。
[0074]
步骤202、对于每一被测井，根据所述被测井的垂直地震剖面测井速度曲线以及声波速度曲线进行拟合操作，得到所述被测井对应的拟合参数。
[0075]
可选的，可对每一个井基于最小二乘法对垂直地震剖面测井速度曲线以及声波速度曲线进行拟合，得到每一个井的拟合参数，从而得到每一个井的回归函数。
[0076]
步骤203、根据多个被测井对应的拟合参数，得到所述目标区域的垂直地震剖面测井速度曲线与声波速度曲线的拟合参数。
[0077]
可选的，根据得到的多个被测井的拟合参数，可以从中选取一个合适的拟合参数为该目标区域的拟合参数。
[0078]
步骤204、根据所述目标区域对应的拟合参数，对所述目标区域除所述多个被测井以外的其它井的声波速度曲线进行校正。
[0079]
可选的，得到目标区域对应的拟合参数后，可得到该目标区域对应的校正公式。根据得到的校正公式，可对该目标区域其他井进行声波测井校正，利用声波测井得到其他井的声波速度曲线，将声波速度带入校正公式中，可实现对其它井的声波速度曲线进行校正。
[0080]
具体地，当计算出某一目标区域的拟合参数后，可得到该区域的校正函数。选取该目标区域多个被测井以外的其它井，对其它井进行声波测井，得到其它井的声波速度曲线，用校正函数对该声波速度曲线进行校正，可得到较正后的声波速度曲线。
[0081]
在本发明实施例中，通过获取目标区域的多个被测井的垂直地震剖面测井速度曲线以及声波速度曲线，对于每一被测井，根据所述被测井的垂直地震剖面测井速度曲线以
及声波速度曲线进行拟合操作，得到所述被测井对应的拟合参数，根据多个被测井对应的拟合参数，得到所述目标区域的垂直地震剖面测井速度曲线与声波速度曲线的拟合参数，根据所述目标区域对应的拟合参数，对所述目标区域除所述多个被测井以外的其它井的声波速度曲线进行校正，从而可以通过统计学手段分析隐含在目标区块的内在规律，根据多个被测井建立目标区域的校正方案，实现声波测井校正的定量描述，在一定程度上弥补了现有校正技术的局限性，大幅提升声波测井的标定精度与可靠性，提高校正效率，降低勘探成本。
[0082]
在上述实施例的基础上，可选的，根据所述被测井的垂直地震剖面测井速度曲线以及声波速度曲线进行拟合操作，得到所述被测井对应的拟合参数，包括：
[0083]
建立垂直地震剖面测井速度曲线与声波速度曲线的回归函数，所述回归函数包括拟合参数；在所述被测井的垂直地震剖面测井速度曲线与声波速度曲线上分别选取多个点，根据所选取的点对应的速度值进行拟合操作，为所述被测井对应的拟合参数赋值。
[0084]
针对被测井对应的垂直地震剖面测井速度曲线和声波速度曲线，建立被测井对应的两条曲线的回归函数，回归函数中包含有拟合参数。
[0085]
具体地，对某一被测井，可根据深度不同选取被测井两条曲线对应的速度值，选取不同深度的速度值，对速度值进行拟合，得到被测井的拟合参数。
[0086]
可选的，可遍历每一被测井两条速度曲线上的速度值进行拟合，得到每一个被测井对应的拟合参数，从而得到多个拟合参数。
[0087]
在本实施例中，通过对每一个被测井两条曲线上的速度进行拟合，得到每一个被测井的拟合参数，量化目标区域声波测井的误差，为声波测井的校正方法提供了一种新的思路。
[0088]
可选的，所述回归函数为
[0089][0090]
拟合操作的收敛条件为取最小值；其中，表示被测井对应的回归函数，a0，a1，
…
，an表示被测井的拟合参数，xi为声波速度曲线上的第i个点所对应的速度值，yi表示垂直地震剖面测井速度曲线上的第i个点所对应的速度值，为拟合公式计算值，ξ表示单井误差值，n表示为大于等于2的正整数，p表示为选取的点数。
[0091]
其中，两条曲线上相同序号的点对应的深度相同。举例来说，在选取某一被测井的两条速度曲线上的速度值时，可在两条速度曲线上选取i＝3时对应深度的速度值，例如，当i＝3时选取的深度为200m，两条速度曲线对应的同一深度的速度值分别选取为3400m/s和4500m/s。当选取50个深度对应的速度值时，可对该被测井的两条速度曲线进行拟合。
[0092]
根据建立的每一个井对应的两条曲线的回归函数，求得每一个井对应的拟合参数。其中，可以根据如下表达式求得：
[0093][0094]
其中，m表示为被测井的个数，n表示为大于等于2的正整数。
[0095]
可选的，当n＝2时，可以得到回归函数为大多数井的声波速度曲线和垂直地震剖面测井速度曲线关系由二次回归函数就可以概述清楚。
[0096]
当收敛条件取得最小值时，可得到每一个被测井能够使自身拟合误差最小的拟合参数j组。其中，被测井的个数与拟合参数的个数可以是相同的。
[0097]
在本实施例中，有效的将地层流体性质、井筒温压条件及孔隙裂缝等这些对声波时差测量有明显影响却难以量化的地质因素综合考虑，得到每一口井的拟合参数，从而保留每一口井的单井特性。
[0098]
可选的，根据多个被测井对应的拟合参数，得到所述目标区域的垂直地震剖面测井速度曲线与声波速度曲线的拟合参数，包括：
[0099]
将所述多个被测井对应的拟合参数代入以下公式：
[0100][0101]
其中，δ
l
表示第l组拟合参数下m个被测井的整体误差，表示第l组拟合参数对第j个被测井进行计算得到的单井误差，ξj表示第j个被测井取得最佳拟合参数的最小单井误差，m表示为被测井的个数；根据收敛条件，从所述多个被测井对应的拟合参数中选择所述目标区域对应的拟合参数；其中，所述收敛条件为δ
l
取最小值。
[0102]
可将每个井对应的拟合参数依次带入公式(3)中，得到每个井在该公式下的误差δ
l
，通过求取最小的误差δ
l
，得到该目标区域对应的拟合参数。
[0103]
在本实施例中，通过代入法找到目标区域的拟合参数，使得找到的拟合参数更加符合该目标区域的地质特征，拟合出声波测井的校正，得到于垂直地震剖面测井速度相一致的声波测井曲线，大幅提升声波测井标定精度和可靠性。在此基础上，可建立多目标区域的校正公式，最大程度减少垂直地震剖面测井次数，降低勘探成本。
[0104]
可选的，对于每一被测井，根据所述被测井的垂直地震剖面测井速度曲线以及声波速度曲线进行拟合操作，得到所述被测井对应的拟合参数，包括：
[0105]
建立所述目标区域的垂直地震剖面测井速度曲线与声波速度曲线的交互图，所述交互图的横轴为声波测井速度，纵轴为垂直地震剖面测井速度，所述交互图中的交互点的位置由多个被测井的垂直地震剖面测井速度曲线以及声波速度曲线确定；根据所述交互图中的各个交互点的位置，确定所述各个交互点的集中程度是否满足预设条件；若满足，则对于每一被测井，根据所述被测井的垂直地震剖面测井速度曲线以及声波速度曲线进行拟合操作，得到所述被测井对应的拟合参数。
[0106]
根据目标区域所测井对应的垂直地震剖面测井速度曲线和声波速度曲线，利用jason软件做两个曲线的交互图，判断交互图上的交互点的位置是否满足预设条件，其中预设条件可以根据实际情况进行设定。具体地，可在横轴为声波测井速度，纵轴为垂直地震剖
面测井速度的交互图上，划定一条y＝x的直线，将平面划定为上下两部分，设定上面部分为第一区域，下面部分为第二区域，判定当该目标区域在交互图上第一区域的交互点占总体的百分比超过预设阈值时，该地域的声波速度比垂直地震剖面测井速度大的值是具有某种关联的，得到该目标区域的被测井是满足条件的，可以构建该目标区域的校正函数，对两条曲线进行拟合，得到对应的拟合参数。
[0107]
图3为本发明实施例提供的一种校正前的多井声波测井速度与垂直地震剖面测井速度交互图。如图3所示，该交互图的横轴表示为声波速度，单位为m/s，纵轴表示为垂直地震剖面测井速度，单位为m/s。有一条y＝x的直线将交互图划分为两部分，分别表示为第一区域和第二区域，其中，第一区域表示为vsp速度》声波速度，意为该区域的交互点的垂直地震剖面测井速度大于声波测井速度。第二区域表示为vsp速度《声波速度，意为该象限处的交互点垂直地震剖面测井速度小于声波测井速度。当交互图中第一区域的交互点占总体交互点的一定百分比时，相对集中于y＝x曲线的上侧，可得到该目标区域的被测井的交互点中垂直地震剖面测井速度大于声波测井速度的交互点占有较大的比例，两者之间具有关联关系，可对该目标区域的每一个井进行曲线拟合，得到拟合参数。
[0108]
具体地，可以设定百分比为百分之八十，当第一区域的交互点占总体交互点的百分之八十时可进行进一步的曲线拟合。
[0109]
由于在声波测井时，受到的影响因素比较多，通常得到的声波速度相较于受影响因素较少的地震剖面测井速度偏小。因此，当声波测井速度相比较于地震剖面测井速度偏小的交互点占一定百分比时，可知两者的速度曲线满足一定的关系，可进行进一步的校正。
[0110]
在本实施例中，对被测井的声波测井速度与垂直地震剖面测井速度进行交互分析，判断交互图上交互点的位置关系，得到该目标区域的被测井满足统计学基础，得到后续校正工作具有可行性。
[0111]
可选的，对所述目标区域除所述多个被测井以外的其它井的声波速度曲线进行校正，包括：
[0112]
获取所述目标区域中待分析的井的声波速度曲线；根据所述待分析的井的声波速度曲线以及所述目标区域对应的拟合参数，确定所述待分析的井对应的垂直地震剖面测井速度曲线；根据所确定的垂直地震剖面测井速度曲线，对所述待分析的井进行井震标定。
[0113]
根据得到的该目标区域对应的拟合参数，得到该目标区域的校正函数。可通过校正函数对该目标区域其他需要分析的井进行校正。先利用声波测井获取待分析井的声波时差速度，再根据声波时差速度得到声波速度，将声波速度带入校正函数中，得到校正后的声波速度曲线，也就是作为拟合后的该待分析井的拟垂直地震剖面测井速度曲线。可将校正后的拟垂直地震剖面测井速度曲线应用于井震标定工作。
[0114]
在本实施例中，通过计算该目标区域的拟合参数，对声波测井数据进行校正补偿，在节省复杂区块垂直地震剖面测井速度曲线测井成本的同时，在一定程度上弥补了现有校正手段的局限性和不足。
[0115]
图4为本发明实施例提供的一种校正后多井声波测井速度与垂直地震剖面测井速度交互图。如图4所示，根据图3交互图分析，该目标区域交互图上的交互点满足预设的条件，可对该目标区域的被测井进行拟合，得到该目标区域的拟合参数，从而得到校正公式为f(x)＝6.16
×
10-5
x2+0.38x+1548。通过对图3中的声波速度进行校正后，可得到交互点是均
匀分布在y＝x曲线两侧的。当交互点均匀分布在y＝x曲线两侧时，表明校正后的声波速度同垂直地震剖面测井速度相差较小，校正后的声波速度可等同于垂直地震剖面测井速度，可知该校正公式具有一定的可行性。
[0116]
图5a为本发明实施例提供的校正前多井声波测井速度与垂直地震剖面测井速度概率分布图。如图5a所示，用一种立体方式显示声波速度和地震剖面测井速度的关系，在校正前，图中的数据点更集中于图中的右半部份，在图中右侧第二区域可以得到vsp速度》声波速度，左侧第一区域为vsp速度《声波速度。显然，该目标区域的被测井的数据点在立体图上是靠右半部分相对集中，不分散的。因此，满足预设条件，可进行曲线拟合，得到该目标区域对应的拟合参数。
[0117]
图5b为本发明实施例提供的校正后多井声波测井速度与垂直地震剖面测井速度概率分布图。如图5b所示，通过对该目标区域进行校正，得到校正后图中的数据点分布到中间部分，相当于两个速度基本持平。相比较于图5a可知，在进行校正后的概率分布图中样点数具有明显的变化，并且位置比较集中于第一区域与第二区域之间。
[0118]
在一种可实现的实施例中，对西部某地层进行标定统层时发现利用声波测井制作的合成记录与地震数据不匹配，误差较大，无法满足标定工作。在用声波测井数据对比垂直地震剖面测井数据时，发现两者速度基线有明显的分离现象，从时深关系上对比发现，垂直地震剖面时深关系斜率更大，二者在底部出现明显分离，标定结果相差较大，因此单纯利用声波测井数据进行标定会出现较大的问题。
[0119]
在对标定误差进行综合分析时，发现复杂的地质条件与地层是导致声波测量产生误差的重要因素，这样的影响因素是综合的，常见的校正方法通常是针对单一影响因素进行校正，因此在校正原理上无法满足要求，需要构建新的校正方法。在本发明实施例的情况下，先对构造单元a的所有井的垂直地震剖面与声波测井进行了统计分析，在大数据、多样本的统计下，发现所有井的声波测井速度均小于垂直地震剖面测井速度，且样点分布集中，满足本发明的基础条件，通过可行性验证。并对构造单元a区域内的所有井进行拟合求取，得到该区域的校正公式。再对该区域的所有井进行声波测井校正。验证得到校正后的测井速度曲线与垂直地震剖面速度曲线基线高度重合，校正后的样本点在分布上明显转移到第一区域与第二区域之间，分布密度相对均匀。
[0120]
图6a为本发明实施例提供的校正前声波速度曲线与垂直地震剖面速度曲线叠合图。如图6a所示，在校正前，从时深曲线上来看，三个被测井的地震剖面速度相比较于声波速度偏右，说明地震剖面速度偏大，声波速度偏小。
[0121]
图6b为本发明实施例提供的校正后声波速度曲线与垂直地震剖面速度曲线叠合图。如图6b所示，在校正后，从时深曲线上来看，三个被测井的地震剖面速度与声波速度二者基本重合。说明根据校正公式得到的声波速度曲线是可以很好的拟合垂直地震剖面曲线。
[0122]
图7a为本发明实施例提供的校正前声波测井与垂直地震剖面测井时深曲线关系图。如图7a所示，某地区2井，某地区101井和某地区3井三个被测井的时深曲线从顶部到底部的差距是明显变大的，由于穿过相同深度的岩层时，垂直地震剖面的测井速度相较于声波速度较快，时间较短，随着深度的累加，底部的差距越来越大。其中，时深关系曲线表示将速度曲线转换成时间和深度的关系。
[0123]
图7b为本发明实施例提供的校正后声波测井与垂直地震剖面测井时深曲线关系图。如图7b所示，某地区2井，某地区101井和某地区3井在校正后，声波速度和垂直地震剖面测井速度近似相等，则两者时间和深度上的对比也相同，因此，声波时深曲线与垂直地震剖面测井时深曲线具有高度的重合。
[0124]
图8为本发明实施例提供的校正前后某地区2井综合标定剖面的对比图。如图8所示，校正后的合成记录比校正前的合成记录较短，校正后制作的合成记录其波组特征、全井段的波组关系与地震剖面匹配良好，校正之后的速度曲线用来作标定的时候，标定效果更好，标定准确性与可信度得到大幅提升。
[0125]
图9为本发明实施例提供的校正前后某地区1井综合标定剖面的对比图。如图9所示，校正后的合成记录比校正前的合成记录较短，校正后制作的合成记录其波组特征、全井段的波组关系与地震剖面匹配良好，校正之后的速度曲线用来做标定的时候，标定效果更好，标定准确性与可信度得到大幅提升。
[0126]
图10为本发明实施例提供的声波速度校正图版。如图10所示，根据本发明提供的校正方法，可对构造单元a、构造单元b和构造单元c建立校正图版，得到每一地区的校正函数。该校正函数可以用于对相关地区的声波速度曲线进行校正。
[0127]
图11为本发明实施例提供的声波测井校正装置的结构示意图。如图11所示，本实施例的图11为本发明实施例提供的声波测井校正装置的结构示意图装置可以包括：获取模块1101、第一得到模块1102、第二得到模块1103和校正模块1104。
[0128]
获取模块1101，用于获取目标区域的多个被测井的垂直地震剖面测井速度曲线以及声波速度曲线；
[0129]
第一得到模块1102，用于对于每一被测井，根据所述被测井的垂直地震剖面测井速度曲线以及声波速度曲线进行拟合操作，得到所述被测井对应的拟合参数；
[0130]
第二得到模块1103，用于根据多个被测井对应的拟合参数，得到所述目标区域的垂直地震剖面测井速度曲线与声波速度曲线的拟合参数；
[0131]
校正模块1104，用于根据所述目标区域对应的拟合参数，对所述目标区域除所述多个被测井以外的其它井的声波速度曲线进行校正。
[0132]
可选的，所述第一得到模块1102在根据所述被测井的垂直地震剖面测井速度曲线以及声波速度曲线进行拟合操作，得到所述被测井对应的拟合参数时，具体用于：
[0133]
建立垂直地震剖面测井速度曲线与声波速度曲线的回归函数，所述回归函数包括拟合参数；
[0134]
在所述被测井的垂直地震剖面测井速度曲线与声波速度曲线上分别选取多个点，根据所选取的点对应的速度值进行拟合操作，为所述被测井对应的拟合参数赋值。
[0135]
可选的，回归函数为
[0136][0137]
拟合操作的收敛条件为取最小值；
[0138]
其中，表示被测井对应的回归函数，a0，a1，
…
，an表示被测井的拟合参数，xi表示声波速度曲线上的第i个点所对应的速度值，yi表示垂直地震剖面测井速度曲线上的第i
个点对应的速度值，表示拟合公式计算值，ξ表示单井误差值，n表示为大于等于2的正整数，p表示为选取的点数。
[0139]
可选的，所述第二得到模块1103具体用于：
[0140]
将所述多个被测井对应的拟合参数代入以下公式：
[0141][0142]
其中，δ
l
表示第l组拟合参数下m个被测井的整体误差，表示第l组拟合参数对第j个被测井进行计算得到的单井误差，ξj表示第j个被测井取得最佳拟合参数的最小单井误差，m表示为被测井的个数；
[0143]
根据收敛条件，从所述多个被测井对应的拟合参数中选择所述目标区域对应的拟合参数；
[0144]
其中，所述收敛条件为δ
l
取最小值。
[0145]
可选的，所述第一得到模块1102具体用于：
[0146]
建立所述目标区域的垂直地震剖面测井速度曲线与声波速度曲线的交互图，所述交互图的横轴为声波测井速度，纵轴为垂直地震剖面测井速度，所述交互图中的交互点的位置由多个被测井的垂直地震剖面测井速度曲线以及声波速度曲线确定；
[0147]
根据所述交互图中的各个交互点的位置，确定所述各个交互点的集中程度是否满足预设条件；
[0148]
若满足，则对于每一被测井，根据所述被测井的垂直地震剖面测井速度曲线以及声波速度曲线进行拟合操作，得到所述被测井对应的拟合参数。
[0149]
可选的，所述校正模块1104在对所述目标区域除所述多个被测井以外的其它井的声波速度曲线进行校正时，具体用于：
[0150]
获取所述目标区域中待分析的井的声波速度曲线；
[0151]
根据所述待分析的井的声波速度曲线以及所述目标区域对应的拟合参数，确定所述待分析的井对应的垂直地震剖面测井速度曲线；
[0152]
根据所确定的垂直地震剖面测井速度曲线，对所述待分析的井进行井震标定。
[0153]
本发明实施例提供的声波测井校正装置，可以实现上述如图1-图10所示的实施例的声波测井校正方法，其实现原理和技术效果类似，此处不再赘述。
[0154]
图12为本发明实施例提供的一种声波测井校正设备的结构示意图。如图12所示，本实施例提供的设备，可以包括：至少一个处理器121和存储器122；
[0155]
所述存储器122存储计算机执行指令；
[0156]
所述至少一个处理器121执行所述存储器122存储的计算机执行指令，使得所述至少一个处理器121执行上述任一实施例所述的方法。
[0157]
其中，存储器122和处理器121可以通过总线123连接。
[0158]
本实施例提供的设备的具体实现原理和效果可以参见图1-图10所示实施例对应的相关描述和效果，此处不做过多赘述。
[0159]
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质中存储有计算机执行指令，所述计算机执行指令被处理器执行时用于实现本发明任一实施例所述
的声波测井校正方法。
[0160]
本发明实施例还提供一种计算机程序产品，包括计算机程序，所述计算机程序被处理器执行时实现本发明任一实施例所述的声波测井校正方法。
[0161]
在本发明所提供的几个实施例中，应该理解到，所揭露的装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，模块的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个模块或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或模块的间接耦合或通信连接，可以是电性，机械或其它的形式。
[0162]
所述作为分离部件说明的模块可以是或者也可以不是物理上分开的，作为模块显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案。
[0163]
另外，在本发明各个实施例中的各功能模块可以集成在一个处理单元中，也可以是各个模块单独物理存在，也可以两个或两个以上模块集成在一个单元中。上述模块成的单元既可以采用硬件的形式实现，也可以采用硬件加软件功能单元的形式实现。
[0164]
上述以软件功能模块的形式实现的集成的模块，可以存储在一个计算机可读取存储介质中。上述软件功能模块存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)或处理器执行本发明各个实施例所述方法的部分步骤。
[0165]
应理解，上述处理器可以是中央处理单元(central processing unit，简称cpu)，还可以是其它通用处理器、数字信号处理器(digital signal processor，简称dsp)、专用集成电路(application specific integrated circuit，简称asic)等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合发明所公开的方法的步骤可以直接体现为硬件处理器执行完成，或者用处理器中的硬件及软件模块组合执行完成。
[0166]
存储器可能包含高速ram存储器，也可能还包括非易失性存储nvm，例如至少一个磁盘存储器，还可以为u盘、移动硬盘、只读存储器、磁盘或光盘等。
[0167]
总线可以是工业标准体系结构(industry standard architecture，简称isa)总线、外部设备互连(peripheral component interconnect，简称pci)总线或扩展工业标准体系结构(extended industry standard architecture，简称eisa)总线等。总线可以分为地址总线、数据总线、控制总线等。为便于表示，本发明附图中的总线并不限定仅有一根总线或一种类型的总线。
[0168]
上述存储介质可以是由任何类型的易失性或非易失性存储设备或者它们的组合实现，如静态随机存取存储器(sram)，电可擦除可编程只读存储器(eeprom)，可擦除可编程只读存储器(eprom)，可编程只读存储器(prom)，只读存储器(rom)，磁存储器，快闪存储器，磁盘或光盘。存储介质可以是通用或专用计算机能够存取的任何可用介质。
[0169]
一种示例性的存储介质耦合至处理器，从而使处理器能够从该存储介质读取信息，且可向该存储介质写入信息。当然，存储介质也可以是处理器的组成部分。处理器和存储介质可以位于专用集成电路(application specific integrated circuits，简称asic)中。当然，处理器和存储介质也可以作为分立组件存在于电子设备或主控设备中。
[0170]
本领域普通技术人员可以理解：实现上述各方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成。前述的程序可以存储于一计算机可读取存储介质中。该程序在执行时，执行包括上述各方法实施例的步骤；而前述的存储介质包括：rom、ram、磁碟或者光盘等各种可以存储程序代码的介质。
[0171]
本领域技术人员在考虑说明书及实践这里公开的发明后，将容易想到本发明的其它实施方案。本发明旨在涵盖本发明的任何变型、用途或者适应性变化，这些变型、用途或者适应性变化遵循本发明的一般性原理并包括本发明未公开的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的，本发明的真正范围和精神由下面的权利要求书指出。
[0172]
应当理解的是，本发明并不局限于上面已经描述并在附图中示出的精确结构，并且可以在不脱离其范围进行各种修改和改变。本发明的范围仅由所附的权利要求书来限制。