本公开涉及但不限于地层评价领域，特别是涉及到一种地层波走时的确定方法和装置。
背景技术：
在地层评价领域，从声波测井阵列声波信号中提取地层波的走时对地层评价具有重要意义，如获取地层波走时才可以对地层波进行幅度或衰减计算，从而评价地层的岩性或岩石中孔隙流体的性质；地层波的走时还可以用于地层的径向速度剖面成像，从而评价岩石的脆性，指导储层的压裂作业。
在电缆声波测井中，地层波的走时计算较为容易，采用能量比法或利用地层波的速度计算走时。在随钻声波测井中，采用传统的电缆声波计算地层波走时的方法存在精确度不高的问题，不再适用。
技术实现要素：
以下是对本文详细描述的主题的概述。本概述并非是为了限制权利要求的保护范围。
本公开实施例提供了一种地层波走时的确定方法和装置，应用于随钻声波测井仪，显著提升了随钻声波测井中计算随钻声波地层波走时的精确度。
本公开提供一种地层波走时的确定方法，应用于随钻声波测井仪，包括，
获取随钻声波测井仪在设定的测量深度测量到的多个接收器对应的多个原始随钻声波单极波形，并根据所述多个原始随钻声波单极波形，确定阵列滤波波形；
根据所述阵列滤波波形，确定地层波的纵波慢度和横波慢度；
根据所述地层波的纵波慢度，对所述阵列滤波波形进行叠加处理后，形成实测阵列波形；
根据所述阵列滤波波形、所述纵波慢度、所述横波慢度和设定的计算参数，计算理论阵列波形；其中，所述计算参数包括：声源距离、流体声速、地层密度和井眼半径；
根据所述实测阵列波形和所述理论阵列波形，判断是否满足收敛条件，如果满足，则根据所述理论阵列波形提取地层纵波走时；如果不满足，则修改所述计算参数中的地层密度和/或井眼半径，重新计算所述理论阵列波形，并判断是否满足收敛条件。
一些示例性实施例中，所述根据所述多个原始随钻声波单极波形，确定阵列滤波波形，包括：
根据所述随钻声波测井仪的钻铤波隔声阻带，设置钻铤波带通滤波区间；
根据所述钻铤波带通滤波区间，对所述多个原始随钻声波单极波形分别进行有限冲击响应fir带通滤波；
根据滤波后的多个波形，确定所述阵列滤波波形。
一些示例性实施例中，所述根据所述阵列滤波波形，确定地层波的纵波慢度和横波慢度，包括：
根据所述阵列滤波波形，采用时间-慢度相关分析，确定所述地层波的纵波慢度和所述地层波的横波慢度。
一些示例性实施例中，所述根据所述地层波的纵波慢度，对所述阵列滤波波形进行叠加处理后，形成实测阵列波形，包括：
对所述阵列滤波波形中的每一个滤波后波形，分别在时间上根据所述地层波的纵波慢度采用向前和向后叠加处理后，形成所述实测阵列波形。
一些示例性实施例中，所述多个原始随钻声波单极波形包括n个原始随钻声波单极波形，n为大于1的整数；
所述对所述阵列滤波波形中的每一个滤波后波形，分别在时间上根据所述地层波的纵波慢度采用向前和向后叠加处理后，形成所述实测阵列波形，包括：
对第l个滤波后波形，根据所述地层波的纵波慢度采用向前和向后叠加处理后得到的波形如下：
其中，zl表示第l个接收器距离声源的距离。wn(t)表示第n个接收器的滤波后的波形，
根据n个叠加处理后得到的波形，确定所述实测阵列波形。
一些示例性实施例中，所述根据所述阵列滤波波形、所述纵波慢度、所述横波慢度和设定的计算参数，计算理论阵列波形，包括：
从所述阵列滤波波形中截取预设周期的地层纵波信号，根据所述纵波慢度、所述横波慢度、所述设定的计算参数和所截取的预设周期的地层纵波信号，采用实轴积分法计算所述理论阵列接收波形。
一些示例性实施例中，所述根据所述实测阵列波形和所述理论阵列波形，判断是否满足收敛条件，包括：
计算所述实测阵列波形和所述理论阵列波形的均方根误差；
如果所述均方根误差小于预设的误差阈值，则确定满足收敛条件；如果大于或等于预设的误差阈值，则确定不满足收敛条件。
一些示例性实施例中，所述根据所说理论阵列波形提取地层纵波走时，包括：
根据所述理论阵列波形，采用能量比法或阈值法提取所述地层纵波走时。
本公开还提供一种地层波走时的确定方法，应用于随钻声波测井仪，包括，
在待测量井的深度范围内，选定预设数量的测量深度，在每一个测量深度分别执行上述任一种地层波走时的确定方法，确定对应测量深度的地层纵波走时；
根据待测量井的深度范围内的全部测量深度对应的地层纵波走时，确定待测量井的地层波走时曲线。
本公开还提供一种电子装置，包括存储器和处理器，所述存储器中存储有用于进行地层波走时确定的计算机程序，所述处理器被设置为读取并运行所述用于进行地层波走时确定的计算机程序以执行上述任一种地层波走时的确定方法。
在阅读并理解了附图和详细描述后，可以明白其他方面。
附图说明
图1为本公开实施例中随钻声波测井仪进行测量的示意图；
图2为本公开实施例中一种地层波走时的确定方法的流程图；
图3为本公开实施例中滤波后的波形向前向后叠加处理示例；
图4为本公开实施例中滤波后的波形与波形叠加后的对比示例；
图5为本公开实施例中第4个接收器的滤波后的波形与波形叠加后对比示例；
图6为本公开实施例中第4个接收器的叠加后的波形与理论波形对比示例；
图7为本公开实施例中原始波形与地层波增强处理后的对比示例；
图8为本公开实施例中n个接收器波形及地层波走时曲线示例；
图9为本公开另一实施例中一种地层波走时的确定方法的流程图；
图10为本公开另一实施例中一种地层波走时的确定方法的流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图及具体实施例对本发明作进一步的详细描述。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互任意组合。
下述步骤编号不限定特定的执行顺序，根据具体实施例部分步骤能够调整其执行顺序。
实施例一
本公开实施例提供一种地层波走时的确定方法，利用随钻声波测井仪进行测量。本公开实施例中采用的随钻声波测井仪进行测量的场景如图1所示。随钻声波测井仪的探测部分深入待测量井中，探测部分至少包括：声波发射器、n个接收器(r1-rn)。
一种地层波走时的确定方法，其流程如图2所示，包括，
步骤0，获取一个深度测量到的随钻声波单极数据；所述随钻声波单极数据包括随钻声波测井仪的n个接收器在所述深度测量得到的n个原始随钻声波单极波形；
步骤1，根据随钻声波测井仪的钻铤波隔声阻带，设置钻铤波带通滤波区间，对n个原始随钻声波单极波形进行带通滤波，根据滤波后的n个波形确定一阵列滤波波形；
步骤2，对所述阵列滤波波形，采用时间-慢度相关分析法，确定地层波的纵波慢度和横波慢度；
步骤3，以所述地层波的纵波慢度，分别对n个接收器的波形进行向前和向后叠加处理；这里叠加处理的目的是压制钻铤波，增强地层波信号；n个接收波形分别叠加处理后确定一实测阵列波形；
步骤4，计算理论阵列波形；在一定时间窗内从滤波后的波形中截取预设周期的地层纵波波形；根据地层纵波慢度、横波慢度、地层密度及井眼半径计算对应的理论波形，根据n个理论波形确定一理论阵列波形；一些示例性实施例中，预设周期为2或3个周期；一些示例性实施例中，也可以截取其他预设周期数的地层纵波波形；
步骤5，计算所述时间窗内理论阵列波形与实测阵列波形的均方根误差；
步骤6，判断所述均方根误差是否满足收敛条件；如果满足，则执行步骤7；如果不满足，则修改计算参数，返回步骤4，重新计算所述理论阵列波形；所说计算参数包括：地层密度或/和井眼半径。一些示例性实施例中，根据预设的修改步长或规则来修改所述计算参数。
步骤7，根据所述理论阵列波形提取地层纵波走时。
步骤8，判断是否还有下一个深度的随钻声波单极数据，如果有，则返回步骤0，获取下一个深度的随钻声波单极数据，再次执行步骤1-7；如果没有，则执行步骤9；
步骤9，获取待测量井的全部深度的地层纵波走时，根据各深度和对应的地层纵波走时，确定待测量井的地层波走时曲线。
一些示例性实施例中，步骤1包括：对原始随钻单极波形进行有限冲击响应fir带通滤波，获取滤波后的波形。不同尺寸的钻铤，钻铤波的隔声阻带有很大差异。以外径为6.75in的钻铤为例，钻铤波的隔声阻带在为8khz～15khz附近，因此对应设置的带通滤波的范围为8khz～15khz。根据随钻声波测井仪的钻铤相应设置，不限于本示例所例举方式。
一些示例性实施例中，步骤2包括：对滤波波形分别进行时间-慢度相关处理，得到地层的纵波慢度dtc和横波慢度dts(硬地层的横波慢度)。
一些示例性实施例中，步骤3包括：
步骤2中虽然采用了带通滤波处理，但是钻铤波信号在整个频率区间内都存在，所以地层波中仍然有钻铤波信号。为了压制钻铤波，增强地层波信号，对n个接收器的波形在时间上以地层纵波慢度dtc采用了向前和向后叠加处理。
对n个接收器的滤波后波形分别进行如下处理：
以第l个接收器为例，包括：对于第1～l个接收器的波形采用的处理公式为：
而对于第l 1～n个接收器的波形采用的处理公式为：
因此，对于第l个接收器进行向前和向后叠加后的波形为：
上式中zl表示第l个接收器距离声源的距离。wn(t)表示第n个接收器的滤波后的波形，
根据1-n个接收器叠加后的波形确定所述实测阵列波形。
一些示例性实施例中，步骤4包括：
对步骤2中获取的滤波后的波形，截取地层纵波信号的2或3个周期作为声源的激励信号，给定源距离tr、井径、流体声速、地层密度及步骤2中获得的地层纵波慢度、横波慢度，采用实轴积分法计算理论阵列波形s(t)；所述源距离tr表示发射换能器(发射器)中心距离第一个接收换能器(第一个接收器r1)中心的距离。
一些示例性实施例中，还可以用一些数值计算的方法如：有限元或有限差分方法计算论阵列波形s(t)。本公开所提供的方案不限定特定的计算方法，本领域技术人员可以等效替换选择其计算方法。
一些示例性实施例中，步骤5包括：
在所述时间窗内，计算理论阵列波形与实测阵列波形的均方根误差，计算公式为：
sl(t)表示计算得到的第l个接收器对应的理论波形。
一些示例性实施例中，步骤6包括：
步骤5中的均方根误差如果小于预设的误差阈值，则满足收敛条件，进行步骤7，如果不满足收敛条件，则需要改变井径或地层密度，重新进行步骤4，直到均方根误差满足收敛条件，再进行步骤7。
一些示例性实施例中，步骤7包括：
计算步骤6理论计算的阵列波形的走时。
所述理论阵列波形中波形的基线平整，纵波波至点起跳明显，可用常规的能量比法或阈值法等从理论阵列波形中提取当前测量深度的地层波走时。
一些示例性实施例中，步骤9包括：
针对待测量井的全部深度位置，依次执行步骤0-7，得到全部深度位置各自对应的地层波走时；即，对于待测量井的全部深度范围内预设数量个不同深度的测量点，分别执行上述步骤0-7，得到对应的各深度的地层波走时。根据全部深度位置和对应的地层波走时绘制相应的地层波走时曲线。
实施例效果
图3为滤波后的波形(对原始波形进行滤波后的波形)向前和向后叠加示意图，图中以第l个接收器为例，第1～n个接收器所接收到的信号分别以地层纵波慢度传播不同距离后与第l个接收器信号叠加。
图4为滤波后的波形与波形叠加后的对比图，图中实线为滤波后的波形，虚线为向前向后移动后的波形，图中r1～r8个接收器信号分别以地层纵波慢度67μs/ft，按照不同的传播距离，在时间轴上移动后的结果。
图5为第4个接收器滤波后的波形与波形叠加后对比，其中实线为滤波后的波形信号，虚线为波形叠加后的信号，由图可见，滤波后的波形信号中由于钻铤波的干扰，纵波的波至点很难确定，而波形叠加后的信号，地层纵波信号明显幅度增大，更加容易确定地层纵波的波至点。图中矩形窗内为地层纵波信号。
图6为第4个接收器的叠加后的波形与理论波形对比，其中实线为波形叠加后的信号，虚线为理论计算信号，由图可见，在地层纵波的时间窗内，理论计算波形与实测波形吻合很好，且理论波形的波至点清晰，波形基线平整。
图7为采用本公开所述方法计算的原始波形与地层波增强处理后的对比成果曲线。图中第一道为第一个接收器的原始随钻单极变密度曲线，第二道为第一个接收器fir带通滤波后的变密度曲线，第三道为采用本公开方法计算的第一个接收器地层波增强后的变密度曲线，第四道为第一个接收器地层波增强前后的对比波形曲线，虚线为未增强曲线，而实线为增强后的曲线。由对比可见，采用本公开所提供的方法，地层波信号明显幅度增加，信噪比提高。
图8为n个接收器的滤波后的变密度曲线，同时根据本专利计算的理论波形得到的n个接收器地层纵波的走时曲线。由图可见，走时曲线与地层纵波的波至吻合很好。
可以看到，本公开为了克服现有技术存在的缺陷而提出了一种精确计算随钻声波地层波走时的方法。该方法利用了阵列声波相似性，对不同的接收器信号利用地层纵波慢度向前或向后叠加其他接收器信号，从而达到压制钻铤波增强地层波信号的目的。同时根据声波在井孔中的声传播模型，数值模拟阵列接收信号，与实测的阵列信号比较，调整相关参数，使理论波形与实测波形的均方根误差最小，从而精确获取地层波的走时信息。解决了在随钻测井环境中由于环境噪声及钻铤波噪声干扰导致地层波波至点很难提取的问题。
实施例二
本公开实施例还提供一种地层波走时的确定方法，应用于随钻声波测井仪，其流程如图9所示，包括：
步骤901，获取随钻声波测井仪在设定的测量深度测量到的多个接收器对应的多个原始随钻声波单极波形，并根据所述多个原始随钻声波单极波形，确定阵列滤波波形；
步骤902，根据所述阵列滤波波形，确定地层波的纵波慢度和横波慢度；
步骤903，根据所述地层波的纵波慢度，对所述阵列滤波波形进行叠加处理后，形成实测阵列波形；
步骤904，根据所述阵列滤波波形、所述纵波慢度、所述横波慢度和设定的计算参数，计算理论阵列波形；其中，所述计算参数包括：声源距离、流体声速、地层密度和井眼半径；
根据所述实测阵列波形和所述理论阵列波形，判断是否满足收敛条件，如果满足，则执行步骤905，根据所述理论阵列波形提取地层纵波走时；如果不满足，则执行步骤906，修改所述计算参数中的地层密度和/或井眼半径，重新计算所述理论阵列波形，并判断是否满足收敛条件。
一些示例性实施例中，步骤901中所述根据所述多个原始随钻声波单极波形，确定阵列滤波波形，包括：
根据所述随钻声波测井仪的钻铤波隔声阻带，设置钻铤波带通滤波区间；
根据所述钻铤波带通滤波区间，对所述多个原始随钻声波单极波形分别进行有限冲击响应fir带通滤波；
根据滤波后的多个波形，确定所述阵列滤波波形。
一些示例性实施例中，不同尺寸的钻铤，钻铤波的隔声阻带有很大差异。以外径为6.75in的钻铤为例，钻铤波的隔声阻带在为8khz～15khz附近，因此对应设置的带通滤波的范围为8khz～15khz。根据随钻声波测井仪的钻铤相应设置，不限于本示例所例举方式。
一些示例性实施例中，步骤902中所述根据所述阵列滤波波形，确定地层波的纵波慢度和横波慢度，包括：
根据所述阵列滤波波形，采用时间-慢度相关分析，确定所述地层波的纵波慢度和所述地层波的横波慢度。
一些示例性实施例中，步骤903中所述根据所述地层波的纵波慢度，对所述阵列滤波波形进行叠加处理后，形成实测阵列波形，包括：
对所述阵列滤波波形中的每一个滤波后波形，分别在时间上根据所述地层波的纵波慢度采用向前和向后叠加处理后，形成所述实测阵列波形。
一些示例性实施例中，所述多个原始随钻声波单极波形包括n个原始随钻声波单极波形，n为大于1的整数；
所述对所述阵列滤波波形中的每一个滤波后波形，分别在时间上根据所述地层波的纵波慢度采用向前和向后叠加处理后，形成所述实测阵列波形，包括：
对第l个滤波后波形，根据所述地层波的纵波慢度采用向前和向后叠加处理后得到的波形如下：
其中，zl表示第l个接收器距离声源的距离。wn(t)表示第n个接收器的滤波后的波形，
根据n个叠加处理后得到的波形，确定所述实测阵列波形。
一些示例性实施例中，对n个接收器的滤波后波形分别进行如下处理：
以第l个接收器为例，包括：对于第1～l个接收器的波形采用的处理公式为：
而对于第l 1～n个接收器的波形采用的处理公式为：
因此，对于第l个接收器进行向前和向后叠加后的波形为：
上式中zl表示第l个接收器距离声源的距离。wn(t)表示第n个接收器的滤波后的波形，
根据1-n个接收器叠加后的波形确定所述实测阵列波形。
一些示例性实施例中，步骤904中所述根据所述阵列滤波波形、所述纵波慢度、所述横波慢度和设定的计算参数，计算理论阵列波形，包括：
从所述阵列滤波波形中截取预设周期的地层纵波信号，根据所述纵波慢度、所述横波慢度、所述设定的计算参数和所截取的预设周期的地层纵波信号，采用实轴积分法计算所述理论阵列接收波形。
一些示例性实施例中，所述根据所述实测阵列波形和所述理论阵列波形，判断是否满足收敛条件，包括：
计算所述实测阵列波形和所述理论阵列波形的均方根误差；
如果所述均方根误差小于预设的误差阈值，则确定满足收敛条件；如果大于或等于预设的误差阈值，则确定不满足收敛条件。
一些示例性实施例中，根据以下公式计算所述实测阵列波形和所述理论阵列波形的均方根误差：
sl(t)表示计算得到的第l个接收器对应的理论波形。
一些示例性实施例中，步骤905中所述根据所说理论阵列波形提取地层纵波走时，包括：
根据所述理论阵列波形，采用能量比法或阈值法提取所述地层纵波走时。
实施例三
本公开实施例还提供一种地层波走时的确定方法，应用于随钻声波测井仪，其流程如图10所示，包括：
步骤101，在待测量井的深度范围内，选定预设数量的测量深度，在每一个测量深度分别执行上述地层波走时的确定方法，确定对应测量深度的地层纵波走时；
步骤102，根据待测量井的深度范围内的全部测量深度对应的地层纵波走时，确定待测量井的地层波走时曲线。
一些示例性实施例中，预设数量的测量深度(深度位置、深度点)在所述待测量井的深度范围内均匀分布；或者，按照预设的可变深度间隔分布。
实施例四
本公开实施例还提供一种地层波走时的确定装置110，应用于随钻声波测井仪，包括：
滤波模块1101，设置为获取随钻声波测井仪在设定的测量深度测量到的多个接收器对应的多个原始随钻声波单极波形，并根据所述多个原始随钻声波单极波形，确定阵列滤波波形；
慢度确定模块1102，设置为根据所述阵列滤波波形，确定地层波的纵波慢度和横波慢度；
叠加模块1103，设置为根据所述地层波的纵波慢度，对所述阵列滤波波形进行叠加处理后，形成实测阵列波形；
理论波形确定模块1104，设置为根据所述阵列滤波波形、所述纵波慢度、所述横波慢度和设定的计算参数，计算理论阵列波形；其中，所述计算参数包括：声源距离、流体声速、地层密度和井眼半径；
地层波走时确定模块1105，设置为根据所述实测阵列波形和所述理论阵列波形，判断是否满足收敛条件，如果满足，则根据所述理论阵列波形提取地层纵波走时；如果不满足，则修改所述计算参数中的地层密度和/或井眼半径，返回所述理论波形确定模块1104重新计算所述理论阵列波形后，再判断是否满足收敛条件。
一些示例性实施例中，所述确定装置还包括，曲线绘制模块1106；
所述曲线绘制模块1106设置为，根据待测量井的深度范围内全部预设数量的测量深度的地层纵波走时，确定待测量井的地层波走时曲线。
本公开实施例还提供一种电子装置，包括存储器和处理器，其特征在于，所述存储器中存储有用于进行地层波走时确定的计算机程序，所述处理器被设置为读取并运行所述用于进行地层波走时确定的计算机程序以执行上述任一种地层波走时的确定方法。
本公开实施例还提供一种存储介质，所述存储介质中存储有计算机程序，其中，所述计算机程序被设置为运行时执行上述任一种地层波走时的确定方法。
本领域普通技术人员可以理解，上文中所公开方法中的全部或某些步骤、系统、装置中的功能模块/单元可以被实施为软件、固件、硬件及其适当的组合。在硬件实施方式中，在以上描述中提及的功能模块/单元之间的划分不一定对应于物理组件的划分；例如，一个物理组件可以具有多个功能，或者一个功能或步骤可以由若干物理组件合作执行。某些组件或所有组件可以被实施为由处理器，如数字信号处理器或微处理器执行的软件，或者被实施为硬件，或者被实施为集成电路，如专用集成电路。这样的软件可以分布在计算机可读介质上，计算机可读介质可以包括计算机存储介质(或非暂时性介质)和通信介质(或暂时性介质)。如本领域普通技术人员公知的，术语计算机存储介质包括在用于存储信息(诸如计算机可读指令、数据结构、程序模块或其他数据)的任何方法或技术中实施的易失性和非易失性、可移除和不可移除介质。计算机存储介质包括但不限于ram、rom、eeprom、闪存或其他存储器技术、cd-rom、数字多功能盘(dvd)或其他光盘存储、磁盒、磁带、磁盘存储或其他磁存储装置、或者可以用于存储期望的信息并且可以被计算机访问的任何其他的介质。此外，本领域普通技术人员公知的是，通信介质通常包含计算机可读指令、数据结构、程序模块或者诸如载波或其他传输机制之类的调制数据信号中的其他数据，并且可包括任何信息递送介质。
1.一种地层波走时的确定方法，应用于随钻声波测井仪，其特征在于，包括，
获取随钻声波测井仪在设定的测量深度测量到的多个接收器对应的多个原始随钻声波单极波形，并根据所述多个原始随钻声波单极波形，确定阵列滤波波形；
根据所述阵列滤波波形，确定地层波的纵波慢度和横波慢度；
根据所述地层波的纵波慢度，对所述阵列滤波波形进行叠加处理后，形成实测阵列波形；
根据所述阵列滤波波形、所述纵波慢度、所述横波慢度和设定的计算参数，计算理论阵列波形；其中，所述计算参数包括：声源距离、流体声速、地层密度和井眼半径；
根据所述实测阵列波形和所述理论阵列波形，判断是否满足收敛条件，如果满足，则根据所述理论阵列波形提取地层纵波走时；如果不满足，则修改所述计算参数中的地层密度和/或井眼半径，重新计算所述理论阵列波形，并判断是否满足收敛条件。
2.根据权利要求1所述的方法，其特征在于，
所述根据所述多个原始随钻声波单极波形，确定阵列滤波波形，包括：
根据所述随钻声波测井仪的钻铤波隔声阻带，设置钻铤波带通滤波区间；
根据所述钻铤波带通滤波区间，对所述多个原始随钻声波单极波形分别进行有限冲击响应fir带通滤波；
根据滤波后的多个波形，确定所述阵列滤波波形。
3.根据权利要求1或2所述的方法，其特征在于，
所述根据所述阵列滤波波形，确定地层波的纵波慢度和横波慢度，包括：
根据所述阵列滤波波形，采用时间-慢度相关分析，确定所述地层波的纵波慢度和所述地层波的横波慢度。
4.根据权利要求1或2所述的方法，其特征在于，
所述根据所述地层波的纵波慢度，对所述阵列滤波波形进行叠加处理后，形成实测阵列波形，包括：
对所述阵列滤波波形中的每一个滤波后波形，分别在时间上根据所述地层波的纵波慢度采用向前和向后叠加处理后，形成所述实测阵列波形。
5.根据权利要求4所述的方法，其特征在于，
所述多个原始随钻声波单极波形包括n个原始随钻声波单极波形，n为大于1的整数；
所述对所述阵列滤波波形中的每一个滤波后波形，分别在时间上根据所述地层波的纵波慢度采用向前和向后叠加处理后，形成所述实测阵列波形，包括：
对第l个滤波后波形，根据所述地层波的纵波慢度采用向前和向后叠加处理后得到的波形如下：
其中，zl表示第l个接收器距离声源的距离，wn(t)表示第n个接收器的滤波后的波形，
根据n个叠加处理后得到的波形，确定所述实测阵列波形。
6.根据权利要求1或2所述的方法，其特征在于，
所述根据所述阵列滤波波形、所述纵波慢度、所述横波慢度和设定的计算参数，计算理论阵列波形，包括：
从所述阵列滤波波形中截取预设周期的地层纵波信号，根据所述纵波慢度、所述横波慢度、所述设定的计算参数和所截取的预设周期的地层纵波信号，采用实轴积分法计算所述理论阵列接收波形。
7.根据权利要求1所述的方法，其特征在于，
所述根据所述实测阵列波形和所述理论阵列波形，判断是否满足收敛条件，包括：
计算所述实测阵列波形和所述理论阵列波形的均方根误差；
如果所述均方根误差小于预设的误差阈值，则确定满足收敛条件；如果大于或等于预设的误差阈值，则确定不满足收敛条件。
8.根据权利要求1所述的方法，其特征在于，
所述根据所说理论阵列波形提取地层纵波走时，包括：
根据所述理论阵列波形，采用能量比法或阈值法提取所述地层纵波走时。
9.一种地层波走时的确定方法，应用于随钻声波测井仪，其特征在于，包括，
在待测量井的深度范围内，选定预设数量的测量深度，在每一个测量深度分别执行权利要求1至8任一项中所述的方法，确定对应测量深度的地层纵波走时；
根据待测量井的深度范围内的全部测量深度对应的地层纵波走时，确定待测量井的地层波走时曲线。
10.一种电子装置，包括存储器和处理器，其特征在于，所述存储器中存储有用于进行地层波走时确定的计算机程序，所述处理器被设置为读取并运行所述用于进行地层波走时确定的计算机程序以执行所述权利要求1至9任一项中所述的方法。
技术总结