一种三维数字岩心纵横波速度的确定方法与流程

1.本发明涉及岩石物理研究中的数字岩石物理研究领域，具体涉及一种三维数字岩心纵横波速度的确定方法。


背景技术：

2.数字岩心技术在众多岩石物理研究人员的不懈努力下，近二十年已经在油气勘探和开发领域都获得了快速的发展。从数字岩心技术的研究内容来看，主要包括数字岩心建模和岩石物理性质数值计算两个方面。基于数字岩心的岩石物理性质数值计算方面，目前已经开展了岩石声学性质、电学性质、渗流性质和核磁共振的模拟和计算。在岩石声学性质计算方面，目前主要有声格子方法、有限元以及旋转交错网格有限差分方法两大类。
3.声格子模型在格子气模型基础上引入了变速粒子，用来模拟复杂介质中的波动现象。该模型揭示了粒子运动与波动之间的内在联系。与有限差分方法就不同模型的对比计算结果表明，声格子模型可以有效地克服常规计算方法在操作过程中对模型的平滑，更适合复杂介质模型的计算。但是由于该方法的计算量大，目前只应用于二维数值模拟研究。
4.基于三维数字岩心计算纵波速度和横波速度目前主要采用有限元方法。makarynska等人采用有限元的方法研究了部分饱和岩石的弹性性质，并把计算结果与低频gassmann-wood方程和高频gassmann-hill方程进行了比较，发现岩石均匀饱和情况下计算结果与低频gassmann-wood方程计算结果一致。旋转交错网格有限差分模拟可以避免普通离散对于高度非均质介质容易导致交错网格失稳的问题。使用此修改网格可以模拟介质中含裂缝、孔隙及没有明确边界条件时的弹性波传播特性。
5.从目前三维数字岩心纵横波速度计算方法来看，主要存在两个方面的问题。一是目前的方法计算量和存储量都非常大、耗时长，在特定计算资源条件下，只能计算较小尺度数字岩心的纵横波速度；二是计算结果只反映特定频率条件下的纵横波速度，无法满足跨频段纵横波速度变化研究。


技术实现要素：

6.本发明的目的在于提供一种三维数字岩心纵横波速度的确定方法，以解决现有技术在三维数字岩心纵横波速度计算中存在的困难。
7.为了实现以上目的，本发明采用以下技术方案：
8.一种三维数字岩心纵横波速度的确定方法，包括以下步骤：
9.为三维数字岩心中的每一种矿物组分和流体赋予纵波速度和横波速度；
10.设定三维数字岩心的底面为声波传播起始点；
11.从声波传播起始点体素开始，设定数字岩心中与所述声波传播起始点体素相邻的一个体素为中心体素，找出与所述中心体素相邻的体素；
12.确定声波由所述相邻的体素传播到所述中心体素的最短时间；
13.确定所述中心体素与所述声波传播起始点之间是否存在同一相中的直线路径；
14.如果存在所述同一相中的直线路径，将该直线路径对应的时间与所述声波从所述相邻的体素传播到所述中心体素的最短时间中的较短时间确定为所述中心体素的初至时间；如果不存在所述同一相中的直线路径，将所述声波从所述相邻的体素传播到所述中心体素的最短时间确定为中心体素的初至时间；
15.接着确定与所述中心体素相邻的下一个体素的初至时间，直到遍历所有体素；
16.取三维数字岩心顶面中初至时间的最小值作为纵波或者横波传播的初至时间，进而确定三维数字岩心的纵波或横波速度。
17.根据本发明的确定方法，优选地，所述为三维数字岩心中的每一种矿物组分和流体赋予纵波速度和横波速度的步骤中，根据三维数字岩心骨架矿物组分和孔隙流体类型，分别为每一种矿物组分和流体赋予纵波速度和横波速度。
18.根据本发明的确定方法，优选地，所述设定三维数字岩心的底面为声波传播起始点的步骤中，所述声波传播起始点根据实际需要设定为点源、线源或者面源。
19.根据本发明的确定方法，优选地，所述声波传播起始点为点源时，将三维数字岩心底面上该点源所在的体素，设置初至时间为0；
20.所述声波传播起始点为线源时，将三维数字岩心底面上该线源上所有体素，设置初至时间为0；
21.所述声波传播起始点为面源时，将三维数字岩心底面上该面源上所有体素，设置初至时间为0。
22.根据本发明的确定方法，优选地，与所述中心体素相邻的体素包括：与所述中心体素面面相邻、棱相邻、点相邻的所有体素。
23.根据本发明的确定方法，优选地，所述找出与所述中心体素相邻的体素的步骤包括：
24.设定所述中心体素的笛卡尔坐标为(x0，y0，z0)；
25.找出与所述中心体素面相邻的6个体素，其笛卡尔坐标分别为：(x0+1，y0，z0)、(x
0-1，y0，z0)、(x0，y0+1，z0)、(x0，y
0-1，z0)、(x0，y0，z0+1)、(x0，y0，z
0-1)；
26.找出与所述中心体素棱相邻的12个体素，其笛卡尔坐标分别为：(x0+1，y0+1，z0)、(x0+1，y
0-1，z0)、(x
0-1，y0+1，z0)、(x
0-1，y
0-1，z0)、(x0+1，y0，z0+1)、(x0+1，y0，z
0-1)、(x
0-1，y0，z0+1)、(x
0-1，y0，z
0-1)、(x0，y0+1，z0+1)、(x0，y0+1，z
0-1)、(x0，y
0-1，z0+1)、(x0，y
0-1，z
0-1)；
27.找出与中心体素点相邻的8个体素，其笛卡尔坐标分别为：(x0+1，y0+1，z0+1)、(x0+1，y0+1，z
0-1)、(x0+1，y
0-1，z0+1)、(x0+1，y
0-1，z
0-1)、(x
0-1，y0+1，z0+1)、(x
0-1，y0+1，z
0-1)、(x
0-1，y
0-1，z0+1)、(x
0-1，y
0-1，z
0-1)。
28.根据本发明的确定方法，优选地，所述确定声波由所述相邻的体素传播到所述中心体素的最短时间的步骤，根据费马原理进行。
29.根据本发明的确定方法，优选地，所述确定声波由所述相邻的体素传播到所述中心体素的最短时间的步骤包括：
30.对于与所述中心体素面相邻的体素，设定其中一个面相邻体素的纵波初至时间为t
pf
，则纵波由该体素传播到所述中心体素的时间t
pf0
为：
[0031][0032]
式中，a为三维数字岩心的分辨率；v
p0
为中心体素的纵波速度；v
pf
为该面相邻体素的纵波速度；
[0033]
对于与所述中心体素棱相邻的体素，设定其中一个棱相邻体素的纵波初至时间为t
pl
，则纵波由该体素传播到所述中心体素的时间t
pl0
为：
[0034][0035]
式中，v
pl
为该棱相邻体素的纵波速度；
[0036]
对于与所述中心体素点相邻的体素，设定其中一个点相邻体素的纵波初至时间为t
pp
，则纵波由该体素传播到中心体素的时间t
pp0
为：
[0037][0038]
式中，v
pp
为该点相邻体素的纵波速度；
[0039]
根据费马原理，取纵波由所有26个面相邻、棱相邻和点相邻体素传播到所述中心体素的最短时间作为纵波传播到所述中心体素的最短时间；
[0040]
将以上纵波速度改为横波速度，即得到横波传播所述中心体素的最短时间。
[0041]
根据声波在同一相里直线传播时间最短原则，本发明方法中首先确定中心体素与源之间是否存在同一相中的直线路径，如果存在，该直线路径的长度除以该相的速度，即为初至时间。
[0042]
根据本发明的确定方法，优选地，所述确定所述中心体素与所述声波传播起始点之间是否存在同一相中的直线路径的步骤包括：
[0043]
对于与所述声波传播起始点在同一底面中的体素来说，确定方法如下：
[0044]
设所述中心体素笛卡尔坐标为(x0，y0，0)，声波传播起始点体素笛卡尔坐标为(xi，yi，0)；
[0045]
以所述中心体素和所述声波传播起始点体素为直线上的两点，计算直线的斜率k和截距b，分别为：
[0046]
k＝(y
0-yi)/(x
0-xi)x0≠xi[0047]
b＝y
i-k
·
xi[0048]
根据所述斜率和截距，确定所述中心体素和所述声波传播起始点体素之间的直线经过的所有体素对应的笛卡尔坐标；这些体素与所述中心体素和所述声波传播起始点体素一起组成一个集合；
[0049]
判断所述集合中的所有体素是否为同一种矿物组分或流体类型，即是否为同一相；如果是同一相，则说明所述中心体素与所述声波传播起始点之间存在同一相中的直线路径；
[0050]
针对所述声波传播起始点为面源或线源，更换所述声波传播起始点体素，重复以上步骤，找出所述中心体素与所述声波传播起始点之间的其它同一相中的直线路径；从所有的同一相中直线路径，选出最短的直线路径作为所述中心体素与所述声波传播起始点之
间的同一相中的直线路径；
[0051]
对于与所述声波传播起始点不在同一底面中的体素来说，确定方法如下：
[0052]
设所述中心体素笛卡尔坐标为(x0，y0，z0)，声波传播起始点体素笛卡尔坐标为(xi，yi，0)；
[0053]
以所述中心体素和所述声波传播起始点体素为直线上的两点，确定三维空间的直线方程：
[0054]
(x-xi)/(x
0-xi)＝(y-yi)/(y
0-yi)＝z/z0；
[0055]
根据所述直线方程，确定所述中心体素和所述声波传播起始点体素之间的直线经过的所有体素对应的笛卡尔坐标；这些体素与所述中心体素和所述声波传播起始点体素一起组成一个集合；
[0056]
判断所述集合中的所有体素是否为同一种矿物组分或流体类型，即是否为同一相；如果是同一相，说明所述中心体素与所述声波传播起始点之间存在同一相中的直线路径；
[0057]
针对所述声波传播起始点为面源或线源，更换所述声波传播起始点体素，重复以上步骤，找出所述中心体素与所述声波传播起始点之间的其它同一相中的直线路径；从所有的同一相中直线路径，选出最短的直线路径作为所述中心体素与所述声波传播起始点之间的同一相中的直线路径。
[0058]
根据本发明的确定方法，优选地，所述进而确定三维数字岩心的纵波或横波速度的步骤包括：
[0059]
用三维数字岩心在z方向的长度，除以所述纵波或者横波传播的初至时间，得到三维数字岩心在z方向的纵波或横波速度。
[0060]
本发明提出的三维数字岩心纵横波速度确定方法计算量小、占用内存小，计算速度块，不但可以用于较大尺度三维数字岩心的纵横波速度计算，而且可以用于计算不同频率条件下的纵横波速度，为开展基于三维数字岩心的频散研究、岩石声学响应机理研究提供了计算手段。该方法还可以推广到标准柱塞岩样数字岩心和其它尺度数字岩心的纵横波速度计算。
附图说明
[0061]
图1为本发明一实施例之三维数字岩心纵横波速度的确定方法流程图。
[0062]
图2a-图2c为本发明一具体实施例之三维数字岩心底面不同类型源的设置示意图。
[0063]
图3a-图3c为本发明一具体实施例之与中心体素相邻的26个体素。
[0064]
图4为本发明一具体实施例之中心体素与声波发射源之间直线路径示意图。
[0065]
图5为本发明一具体实施例之三维数字岩心纵波初至时间场。
[0066]
图6为本发明一具体实施例之三维数字岩心横波初至时间场。
[0067]
图7a-图7b为本发明一具体实施例之三维打印岩心样品纵横波速度实验测试结果与本发明结果的交会图。其中，图7a为本发明实施例计算的纵波速度与实验测试的纵波速度的关系图；图7b为本发明实施例计算的横波速度与实验测试的横波速度的关系图。
具体实施方式
[0068]
为了更清楚地说明本发明，下面结合优选实施例对本发明做进一步的说明。本领域技术人员应当理解，下面所具体描述的内容是说明性的而非限制性的，不应以此限制本发明的保护范围。
[0069]
如图1所示，本发明在此提供一优选实施例，一种三维数字岩心纵横波速度的确定，具体包括以下步骤：
[0070]
s100、根据三维数字岩心骨架矿物组分和孔隙流体类型，分别为每一种矿物组分和流体赋予纵波速度和横波速度。
[0071]
s200、设定三维数字岩心的底面为声波传播起始点，根据实际需要，可以设定为点源、线源或者面源。
[0072]
为了采用费马原理计算三维数字岩心的纵横波速度，需要设定声波传播的起始点。根据实际需要，可以设置三维数字岩心不同的部位作为声波传播的起始点，本发明优选将三维数字岩心的底面作为声波传播的起始点。起始点设置方法如下：
[0073]
s201、如果设置声波传播起始点为点源，将三维数字岩心底面上该点所在的体素，设置初至时间为0。
[0074]
s202、如果设置声波传播起始点为线源，将三维数字岩心底面上该线上所有体素，设置初至时间为0。
[0075]
s203、如果设置声波传播起始点为面源，将三维数字岩心底面上该面上所有体素，设置初至时间为0。
[0076]
s300、从声波传播起始点体素开始，设定数字岩心中与声波起始点体素相邻的一个体素为中心体素，找出与中心体素相邻的体素。
[0077]
为了找出与中心体素相邻的体素，需要对相邻体素进行分类，一类是与中心体素面接触的体素，称为面相邻体素；一类是与中心体素棱相邻的体素，称为棱相邻体素；一类是与中心体素点接触的体素，称为点相邻体素。找出上述三类相邻体素的步骤如下：
[0078]
s301、设定中心体素的笛卡尔坐标为(x0，y0，z0)，首先找出与中心体素面相邻的6个体素，其笛卡尔坐标分别为：(x0+1，y0，z0)、(x
0-1，y0，z0)、(x0，y0+1，z0)、(x0，y
0-1，z0)、(x0，y0，z0+1)、(x0，y0，z
0-1)。
[0079]
s302、然后找出与中心体素棱相邻的12个体素，其笛卡尔坐标分别为：(x0+1，y0+1，z0)、(x0+1，y
0-1，z0)、(x
0-1，y0+1，z0)、(x
0-1，y
0-1，z0)、(x0+1，y0，z0+1)，(x0+1，y0，z
0-1)、(x
0-1，y0，z0+1)、(x
0-1，y0，z
0-1)、(x0，y0+1，z0+1)、(x0，y0+1，z
0-1)、(x0，y
0-1，z0+1)、(x0，y
0-1，z
0-1)。
[0080]
s303、最后找出与中心体素点相邻的8个体素，笛卡尔坐标分别为：(x0+1，y0+1，z0+1)、(x0+1，y0+1，z
0-1)、(x0+1，y
0-1，z0+1)、(x0+1，y
0-1，z
0-1)、(x
0-1，y0+1，z0+1)、(x
0-1，y0+1，z
0-1)、(x
0-1，y
0-1，z0+1)、(x
0-1，y
0-1，z
0-1)。
[0081]
s400、根据费马原理，确定声波由相邻的体素传播到中心体素的最短时间，设定该时间为t0。
[0082]
需要计算纵波或横波按照直线传播方式，由相邻体素传播到中心体素的时间，然后在所有26个相邻体素传播到中心体素的时间集合中，选出最短的传播时间，步骤如下：
[0083]
s401、对于与中心体素面相邻的体素，设定其中一个面相邻体素的纵波初至时间
为t
pf
，则纵波由该体素传播到中心体素的时间t
pf0
为：
[0084][0085]
上式中a为三维数字岩心的分辨率；v
p0
为中心体素的纵波速度；v
pf
为其中一个面相邻体素的纵波速度。其它面相邻体素传播到中心体素的时间的计算方法与上式相同。
[0086]
s402、对于与中心体素棱相邻的体素，设定其中一个棱相邻体素的纵波初至时间为t
pl
，则纵波由该体素传播到中心体素的时间t
pl0
为：
[0087][0088]
上式中v
pl
为其中一个棱相邻体素的纵波速度，其它参数的物理意义同上。其它棱相邻体素传播到中心体素的时间的计算方法与上式相同。
[0089]
s403、对于与中心体素点相邻的体素，设定其中一个点相邻体素的纵波初至时间为t
pp
，则纵波由该体素传播到中心体素的时间t
pp0
为：
[0090][0091]
上式中v
pp
为其中一个点相邻体素的纵波速度，其它参数的物理意义同上。其它点相邻体素传播到中心体素的时间的计算方法与上式相同。
[0092]
s404、根据费马原理，取所有26个面相邻、棱相邻和点相邻体素传播到中心体素的最短时间作为纵波传播到中心体素的时间，计为t0。
[0093]
s405、横波传播到中心体素的最短时间的计算方法与纵波相似，只需要将所有体素的纵波速度改为横波速度即可。
[0094]
s500、根据声波在同一相里直线传播时间最短原则，因此首先确定中心体素与源之间是否存在同一相中的直线路径。
[0095]
为了找出中心体素与源之间是否存在同一相中直线路径，需要分两种情况来进行处理，步骤如下：
[0096]
s501、对于与源在同一底面中的体素来说，确定方法如下：
[0097]
①
设中心体素笛卡尔坐标为(x0，y0，0)，面源(或线源、点源)体素笛卡尔坐标为(xi，yi，0)；
[0098]
②
以中心体素和面源体素为直线上的两点，计算直线的斜率k和截距b，分别为
[0099]
k＝(y
0-yi)/(x
0-xi)x0≠xi[0100]
b＝y
i-k
·
xi[0101]
③
根据斜率和截距，确定中心体素和面源体素之间的直线经过的所有体素对应的笛卡尔坐标。这些体素与中心体素和面源体素一起组成一个集合；
[0102]
④
判断集合中的所有体素是否为同一种矿物组分或流体类型，也就是说是否为同一相。如果是同一相，说明中心体素与源之间存在同一相中的直线路径；
[0103]
⑤
更换面源体素，重复步骤
①
到步骤
④
，找出中心体素与源之间的其它同一相中的直线路径；
[0104]
⑥
从所有的同一相中直线路径，选出最短的路径作为中心体素与源之间的同一相
中直线路径。
[0105]
s502、对于与源不在同一底面中的体素来说，确定方法如下：
[0106]
①
设中心体素笛卡尔坐标为(x0，y0，z0)，面源(或线源、点源)体素笛卡尔坐标为(xi，yi，0)。
[0107]
②
以中心体素和面源体素为直线上的两点，确定三维空间的直线方程：
[0108]
(x-xi)/(x
0-xi)＝(y-yi)/(y
0-yi)＝z/z0[0109]
③
根据直线方程，确定中心体素和面源体素之间的直线经过的所有体素对应的笛卡尔坐标。这些体素与中心体素和面源体素一起组成一个集合。
[0110]
④
判断集合中的所有体素是否为同一种矿物组分或流体类型，也就是说是否为同一相。如果是同一相，说明中心体素与源之间存在同一相中的直线路径。
[0111]
⑤
更换面源体素，重复步骤
①
到步骤
④
，找出中心体素与源之间的其它同一相中的直线路径。
[0112]
⑥
从所有的同一相中直线路径，选出最短的路径作为中心体素与源之间的同一相中直线路径。
[0113]
在本实施例的步骤s107中，在完成了单个体素纵波或横波初至时间的计算后，需要对三维数字岩心中除源之外所有体素计算纵波或横波初至时间，最后实现三维数字岩心纵横波速度的计算，步骤如下：
[0114]
s503、完成三维数字岩心中除源体素以外所有体素纵波和横波初至时间的计算。
[0115]
s504、选取三维数字岩心顶面体素中最小初至时间值作为纵波或者横波传播的初至时间。
[0116]
s505、用三维数字岩心在z方向的长度，除以纵波或者横波传播的初至时间，就可以得到三维数字岩心在z方向的纵波或横波的速度。
[0117]
s600、如果存在同一相中的直线路径，计算最短路径对应的时间，设定该时间为t1，比较时间t0和t1，将最小时间设定为中心体素的初至时间；如果不存在，将t0设定为中心体素的初至时间；
[0118]
s700、接着计算与中心体素相邻的下一个体素的时间，直到所有体素完成计算；取三维数字岩心顶面中最小初至时间值作为纵波或者横波传播的初至时间，计算三维数字岩心的纵波或横波速度。
[0119]
为了对上述三维数字岩心纵横波速度的确定方法更为清楚、直观的解释，下面结合一个具体的实施例来说明，然而值得注意的是该具体实施例仅是为了更好的说明本发明，并不构成对本发明的不当限定。
[0120]
结合步骤s100，三维数字岩心是数字化的岩石，对于岩石中的每一种矿物组分和流体类型，都用特定的数字与之对应。以下表1为针对致密砂岩储层岩石的矿物组分和孔隙中的流体类型，分别给定的纵波速度和横波速度数值。由于这些值是采用超声波测试的高频结果，因此，本次计算的三维数字岩心的纵横波速度也是高频条件下的值。
[0121]
表1针对致密砂岩储层岩石的矿物组分和孔隙中的流体类型分别给定的纵波速度和横波速度数值及其对应的数字化标签
[0122][0123]
结合步骤s200，三维数字岩心中每一个最小化的单元称为体素，与二维图像中的像素对应。体素是一个边长为分辨率的正方体，本实施例中选用的三维数字岩心的分辨率为3.83μm。图2a-图2c中每一个正方体就是一个体素，其中灰色的体素表示声波传播起始点，可以设置不同的起始点类型，例如图2a为点源，由一个单一体素组成。图2b为线源，由若干在同一直线上的体素组成。图2c为面源，由若干在同一面上的体素组成。本实施例选用面源，由9
×
9共81个体素组成。所有面源体素的纵波传播初至时间和横波传播初至时间均设置为0。
[0124]
结合步骤s300和步骤s400，对于三维数字岩心中每一个非边界体素，总可以找出与其相邻的三类体素：面相邻体素、棱相邻体素和点相邻体素。图3a-图3c分别是这三类相邻体素的示意图，其中相邻体素均用灰色表示，中心体素用白色表示。图3a为面相邻体素的示意图，图3a中与中心体素面相邻的体素共有6个。声波从任意一个面相邻体素中心传播到中心体素中心的距离为体素分辨率，在本实施例中为3.84μm。传播时间由步骤s401中公式计算。图3b为棱相邻体素的示意图，图3b中与中心体素棱相邻的体素共有12个。声波从任意一个棱相邻体素中心传播到中心体素中心的距离为体素分辨率乘以在本实施例中为5.43μm。传播时间由步骤s402中公式计算。图3c为点相邻体素的示意图，图3c中与中心体素点相邻的体素共有8个。声波从任意一个点相邻体素中心传播到中心体素中心的距离为体素分辨率乘以在本实施例中为6.65μm。传播时间由步骤s403中公式计算。待所有26个相邻体素传播到中心体素的时间计算完毕，从中选出最小值作为传播到中心体素的时间，计为t0。
[0125]
结合步骤s500和步骤s600，图4底面中白色体素为面源体素，单一的灰色体素p表示三维数字岩心中除面源体素的任意一个体素。由于岩心在数字化过程中，会形成离散化的规则正方体单元，在实际计算面源体素和任意非面源体素之间的传播时间时，所用的距离是折线形式，而非两点之间最短的直线。因此，需要判断面源体素与任意非面源体素之间的直线经过的体素是否为同一相，也就是需要判断图4中直线经过的体素是否同为石英(或粘土、长石、天然气)。如果存在经过同一相的直线，此时面源体素与任意体素之间距离采用两点之间的距离公式计算，然后除以该相的纵波速度或横波速度计算出传播时间，记为t1，一般情况下，t1≤t0。此时将t1作为从源出发传播到中心体素的最短时间。如果不存在经过同一相的直线，将t0作为从源出发传播到中心体素的最短时间。
[0126]
结合步骤s700，图5为本具体实施例最终获得的三维数字岩心纵波初至时间场。图6为本具体实施例最终获得的三维数字岩心横波初至时间场。图5、图6中颜色的深浅表示了初至时间的差异；颜色越深，表示初至时间越短；颜色越浅，表示初至时间越长。选取三维数字岩心顶面体素中最小初至时间值作为纵波或者横波传播的初至时间。用三维数字岩心在z方向的长度，除以纵波或者横波传播的初至时间，就可以得到三维数字岩心在z方向的纵
波或横波的速度。
[0127]
为了验证本发明提出方法的准确性，利用三维打印裂缝样品(打印材料为光敏树脂)的超声波实验室测量的纵横波速度，与本发明提出的方法得到的纵横波速度进行了对比，具体结果见图7a和图7b。图7a为计算的纵波速度与实验测试的纵波速度的关系图；图7b为计算的横波速度与实验测试的横波速度的关系图。图7a和图7b中可以看出，本发明提出方法计算的纵波速度和横波速度与实验室超声波测量的纵波速度和横波速度接近，说明了本发明提出方法的可靠性。
[0128]
本发明提出的三维数字岩心纵横波速度确定方法计算量小、占用内存小，计算速度块，不但可以用于较大尺度三维数字岩心的纵横波速度计算，而且可以用于计算不同频率条件下的纵横波速度，为开展基于三维数字岩心的频散研究、岩石声学响应机理研究提供了计算手段。该方法还可以推广到标准柱塞岩样数字岩心和其它尺度数字岩心的纵横波速度计算。
[0129]
显然，本发明的上述实施例仅仅是为清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定，对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动，这里无法对所有的实施方式予以穷举，凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。