CN104237383A - 一种岩石超声波速度及品质因子同时测量方法及系统 - Google Patents
一种岩石超声波速度及品质因子同时测量方法及系统 Download PDFInfo
- Publication number
- CN104237383A CN104237383A CN201410446685.2A CN201410446685A CN104237383A CN 104237383 A CN104237383 A CN 104237383A CN 201410446685 A CN201410446685 A CN 201410446685A CN 104237383 A CN104237383 A CN 104237383A
- Authority
- CN
- China
- Prior art keywords
- signal
- measured
- rock sample
- quality factor
- rock
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明公开了一种岩石超声波速度及品质因子同时测量方法及系统,其中,该方法包括:选取参考样品,利用超声波脉冲透射测量系统获取参考样品的标准信号,拾取标准起跳时间;选取待测岩石样品,获取待测岩石样品的测量信号,确定测量信号的初始起跳时间;按照时差对测量信号进行时移,获取时移信号;设定待测岩石样品的品质因子初始值及剩余相对时移初始值,将其作用到标准信号进行正演,得到正演信号;在有效时窗内对比时移信号及正演信号,建立目标函数;从品质因子初始值及剩余相对时移初始值开始,利用优化算法搜索目标函数的最小值,获取待测岩石样品的品质因子最优值及剩余相对时移最优值;计算获得待测岩石样品的走时及超声波速度。
Description
技术领域
本发明涉及基于超声波脉冲透射技术的速度及衰减测量领域,尤其涉及一种岩石超声波速度及品质因子同时测量方法及系统。
背景技术
现今,随着地震及测井技术的发展,为了提高储层预测及流体检测的准确性,人们迫切需要更详细地了解储层流体及岩石的弹性波性质,包括速度、衰减、各向异性等方面。
超声波脉冲透射技术(Birch,1960)是目前研究岩石的波传播特征比较常用的实验测量手段。采集的岩芯通常被加工成一定形态和尺寸的样品(例如圆柱体),该样品的孔隙度、渗透率、密度甚至组构特征均可以事先通过其他实验方法获取得到,所研究的岩石样品非常清楚。此外,界面及射线路径可通过人为设定,有些实验装置可实现对温度、压力、饱和流体类型及饱和状态的准确控制(Timur,1968;Gregory,1976;Wang and Nur,1988等)。因此,通过该方法可以建立弹性参数与目标储层参数之间比较可靠的关系,具有实际应用价值。
对于速度测量而言,一般通过拾取信号起跳时间,由下式计算速度V:
其中,T0和T1分别为探头对接信号和岩石样品信号的起跳时间。样品长度L通常可以比较准确地测量得到。
而在衰减介质中,随着传播距离增加,信号的幅度及形态会发生变换,起跳时间拾取存在较大不确定性,导致速度的测量误差可能会超过2%,使得实验结果准确度较低。
另外,地下储层岩石普遍表现出黏弹性特征,尤其是流体饱和岩石,其衰减机制主要归咎于流体与岩石骨架的相互作用。弹性波衰减通常用品质因子Q来进行描述,Q值越小衰减越严重。对于给定的信号源,发射信号的频带宽度相对有限,一般可以认为Q值不随频率发生变化。谱比法是一种非常流行方法来估计Q值。我们可以比较两个不同传播距离时信号的振幅谱,然后从根据斜率来确定Q值(Bourbiē等,1987):
其中,G为包括几何扩散、透射、反射等影响的系数。
等(1979)对上式进行修改,通过对比参考样品和岩石样品的振幅谱成功获得岩石样品的Q值。然而,谱比法的成功应用需要建立在不受干扰波影响的基础上。如果信号范围比较宽,那么样品的边缘反射和多次波都会污染信号的尾部,从而影响计算Q值的准确性。
速度测量中信号起跳时间的拾取受衰减影响,而Q值计算方法中速度又是重要的输入参数,因此很难独自对速度和衰减进行准确测量。很显然,需要发展出一种新的方法,能够准确地从信号中同时提取速度和Q值。
发明内容
由于现有技术中,速度测量中信号起跳时间的拾取受衰减影响,而Q值计算方法中速度又是重要的输入参数,很难独自对速度和衰减进行准确测量,本发明提出了一种同时获取岩石超声波速度及品质因子的方法及系统,解决了现有超声波实验技术不能准确获取岩石速度及衰减的问题。
为达到上述目的,本发明提出了一种岩石超声波速度及品质因子同时测量方法,包括:步骤1,选取一参考样品,利用超声波脉冲透射测量系统获取所述参考样品的标准信号,并拾取所述标准信号的标准起跳时间;步骤2,选取一待测岩石样品,利用所述超声波脉冲透射测量系统获取所述待测岩石样品的测量信号,并与所述标准信号进行相关性分析,根据所述标准起跳时间确定所述测量信号的初始起跳时间;步骤3,按照所述标准起跳时间与所述初始起跳时间的时差对所述测量信号进行时移,获取时移信号;步骤4,设定所述待测岩石样品的品质因子初始值及剩余相对时移初始值,将其作用到所述标准信号进行正演,得到正演信号;步骤5,设定有效时窗,在所述有效时窗内对比所述时移信号及正演信号,建立目标函数;步骤6,从所述品质因子初始值及剩余相对时移初始值开始,利用优化算法搜索所述目标函数的最小值,根据所述目标函数的最小值获取所述待测岩石样品的品质因子最优值及剩余相对时移最优值;步骤7,根据所述初始起跳时间以及剩余相对时移最优值计算获得所述待测岩石样品的走时,并根据所述待测岩石样品的走时计算获得待测岩石样品的超声波速度。
为达到上述目的,本发明还提出了一种岩石超声波速度及品质因子同时测量系统,包括:标准信号获取模块,用于根据一参考样品,通过超声波脉冲透射测量系统获取所述参考样品的标准信号,并拾取所述标准信号的标准起跳时间;测量信号获取模块,用于根据一待测岩石样品,通过所述超声波脉冲透射测量系统获取所述待测岩石样品的测量信号,并与所述标准信号进行相关性分析,根据所述标准起跳时间确定所述测量信号的初始起跳时间;时移模块,用于按照所述标准起跳时间与所述初始起跳时间的时差对所述测量信号进行时移,获取时移信号;正演模块,用于设定所述待测岩石样品的品质因子初始值及剩余相对时移初始值,将其作用到所述标准信号进行正演,得到正演信号;目标函数建立模块,用于设定有效时窗,在所述有效时窗内对比所述时移信号及正演信号,建立目标函数;品质因子最优值及剩余相对时移最优值获取模块,用于从所述品质因子初始值及剩余相对时移初始值开始,利用优化算法搜索所述目标函数的最小值,根据所述目标函数的最小值获取所述待测岩石样品的品质因子最优值及剩余相对时移最优值;超声波速度计算模块,根据所述初始起跳时间以及剩余相对时移最优值计算获得所述待测岩石样品的走时,并根据所述待测岩石样品的走时计算获得待测岩石样品的超声波速度。
本发明提出的岩石超声波速度及品质因子同时测量方法及系统,与现有技术对这两个量单独进行测量的方法不同,本发明在时间域有效信号时窗范围内进行分析,能够同时测量得到岩石样品的超声波速度及品质因子,克服了两者之间的相互作用,还排除了后续干扰波的影响,可有效地提高测量精度,同时大大提高了实验测量的效率。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,并不构成对本发明的限定。在附图中:
图1为本发明一实施例的岩石超声波速度及品质因子同时测量方法流程图。
图2为本发明一实施例岩石超声波速度及品质因子同时测量的数据关系示意图。
图3为本发明一实施例的岩石超声波速度及品质因子同时测量方法中步骤4的详细流程图。
图4为本发明一具体实施例的参考样品及待测岩石样品的超声波信号示意图。
图5为本发明一具体实施例时移后的信号及最优正演理论信号示意图。
图6为本发明一实施例的岩石超声波速度及品质因子同时测量系统结构示意图。
图7为本发明一实施例的正演模块的具体结构示意图。
具体实施方式
以下配合图式及本发明的较佳实施例,进一步阐述本发明为达成预定发明目的所采取的技术手段。
图1为本发明一实施例的岩石超声波速度及品质因子同时测量方法流程图。图2为本发明一实施例岩石超声波速度及品质因子同时测量的数据关系示意图。结合图1及图2所示,该方法包括以下步骤:
步骤1,选取一参考样品,利用超声波脉冲透射测量系统获取参考样品的标准信号S1(t),并拾取标准信号的标准起跳时间T1;
在本实施例中,参考样品的长度L1、密度ρ1、超声波速度V1、品质因子Q1都是已知量。
步骤2,选取一待测岩石样品,利用超声波脉冲透射测量系统获取待测岩石样品的测量信号S2(t),并与标准信号S1(t)进行相关性分析,根据标准起跳时间T1确定测量信号的初始起跳时间T2;
在步骤2中,待测岩石样品的长度L2、密度ρ2为已知量,而其超声波速度V2和品质因子Q2是需要利用本发明后面步骤计算获得的量。
在本实施例中,超声波脉冲透射测量系统利用的超声波信号包括P波、SV波、SH波。
步骤3,按照标准起跳时间T1与初始起跳时间T2的时差ΔT对测量信号S2(t)进行时移,获取时移信号S3(t);其中,ΔT=T2-T1。
步骤4,设定待测岩石样品的品质因子初始值Q0及剩余相对时移初始值Δτ0,将其作用到标准信号S1(t)进行正演,得到正演信号S4(t)。
在步骤4中,如图3所示,正演信号的获取包括以下具体步骤:
步骤41,对标准信号S1(t)进行傅里叶变换得到其相应的标准频谱S1(f);
步骤42,利用以下公式计算正演信号的频谱S4(f):
其中,S4(f)为正演信号的频谱;
G为待测岩石样品相对于参考样品的透射校正系数,利用公式(2)求得;
f为超声波脉冲透射测量系统的测量频率,为已知量;
T0为超声波测量系统的延迟走时,为已知量;
T1为标准起跳时间;
T2为初始起跳时间;
Δτ0为待测岩石样品的剩余相对时移初始值;
Q0为待测岩石样品的品质因子初始值;
Q1为参考样品的品质因子,为已知量;
S1(f)为标准频谱。
另外,待测岩石样品相对于参考样品的透射校正系数G是利用以下公式所得:
其中,ρ1为参考样品的密度,为已知量;
ρ2为待测岩石样品的密度,为已知量;
Z0为超声波脉冲透射测量系统的探头阻抗,为已知量;
V1参考样品的超声波速度,为已知量;
L2为待测岩石样品的长度,为已知量;
步骤43,对正演信号的频谱S4(f)进行傅里叶反变换得到正演信号S4(t)。
步骤5,设定有效时窗W1,在有效时窗W1内对比时移信号S3(t)及正演信号S4(t),建立目标函数Y。
在步骤5中,建立目标函数利用的公式如下:
其中,Y为目标函数;
n为有效时窗W1内的采样点个数;
S3(ti)为时间域时移信号中每个采样点对应的幅值;
S4(ti)为时间域正演信号中每个采样点对应的幅值。
步骤6,从品质因子初始值Q0及剩余相对时移初始值Δτ0开始,利用优化算法搜索目标函数的最小值Ymin,根据目标函数的最小值Ymin获取待测岩石样品的品质因子最优值及剩余相对时移最优值;
在步骤6中,搜索到目标函数的最小值Ymin后,即可通过计算机软件计算出待测岩石样品的品质因子最优值Q2及剩余相对时移最优值Δτ2。
在本实施例中,步骤6利用的优化算法包括:模式搜索法、Rosenbrock方法、单纯形搜索法、Powell方法、最速下降法、牛顿法、共轭梯度法、拟牛顿法、信赖域方法、最小二次法、可行方向法、遗传算法、蚂蚁算法、模拟退火法以及神经网络法等等方法的至少其中之一。
步骤7,根据初始起跳时间T2以及剩余相对时移最优值Δτ2计算获得待测岩石样品的走时并根据待测岩石样品的走时计算获得待测岩石样品的超声波速度V2。
在步骤7中,计算待测岩石样品的走时利用的公式如下:
其中,为待测岩石样品的走时;
Δτ2为步骤6得到的待测岩石样品的剩余相对时移最优值。
在得到待测岩石样品的走时后,即可计算待测岩石样品的超声波速度V2,利用的公式如下:
其中,V2为待测岩石样品的超声波速度;
为公式(4)得到的待测岩石样品的走时数据;
L2为待测岩石样品的长度,为已知量。
为了对上述岩石超声波速度及品质因子同时测量方法进行更为清楚的解释,下面结合一个具体的实施例来进行说明,然而值得注意的是该实施例仅是为了更好地说明本发明,并不构成对本发明不当的限定。
首先,结合步骤1、2,采用休斯敦大学研制的多功能超声波测量系统对一块铝块样品和一块观雾山组石英砂岩样品(蒸馏水饱和)进行测量,得到铝块的标准信号S1(t)和砂岩的测量信号S2(t),如图4所示,图4为本具体实施例的铝块样品(实线)及砂岩样品(虚线)的超声波信号示意图,并拾取信号的起跳点T1和T2。
铝块样品和砂岩样品的基本参数如表1所示:
表1:铝块和砂岩样品的基本参数
其中,由于铝块几乎没有衰减,因此其品质因子可以认为是无穷大。表1中,砂岩样品的纵波速度和纵波品质因子为待测量。
结合步骤3,按照时移量ΔT=3.75μs,对测量信号S2(t)进行向前时移得到时移信号S3(t),如图5中所示的实线;
结合步骤4,设定砂岩样品的品质因子初始值Q0=50,剩余相对时移初始值Δτ0=0μs,从标准信号S1(t)正演得到砂岩样品的正演信号S4(t),正演信号S4(t)如图5中所示的虚线。
结合步骤5,确定有效时窗W1(至少包含S1(t)的前三个半周期),在该时窗内对比时移信号S3(t)及正演信号S4(t),建立目标函数Y。
结合步骤6,从Q0和Δτ0开始,利用优化算法搜索寻找目标函数的最小值Ymin,利用计算机软件计算得到砂岩样品的纵波品质因子Q2=40.95,剩余相对时移最优值Δτ2=0.05μs。
结合步骤7,该超声波测量系统中探头对接时信号的起跳点T0=11.52μs,通过计算得到砂岩样品的走时
最后,计算得到石英砂岩样品的纵波速度
基于同一发明构思,本发明实施例中还提供了一种岩石超声波速度及品质因子同时测量系统,如下面的实施例所述。由于岩石超声波速度和品质因子的获取系统解决问题的原理与岩石超声波速度和品质因子的获取方法相似,因此该系统的实施可以参见前述方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
图6为本发明一实施例的岩石超声波速度及品质因子同时测量系统结构示意图。如图6所示,该系统包括:
标准信号获取模块11,用于根据一参考样品,通过超声波脉冲透射测量系统获取参考样品的标准信号,并拾取标准信号的标准起跳时间;
测量信号获取模块12,用于根据一待测岩石样品,通过超声波脉冲透射测量系统获取待测岩石样品的测量信号,并与标准信号进行相关性分析,根据标准起跳时间确定测量信号的初始起跳时间;
在本实施例中,超声波脉冲透射测量系统利用的超声波信号包括P波、SV波、SH波。
时移模块13,用于按照标准起跳时间与初始起跳时间的时差对测量信号进行时移,获取时移信号;
正演模块14,用于设定待测岩石样品的品质因子初始值及剩余相对时移初始值,将其作用到标准信号进行正演,得到正演信号;
目标函数建立模块15,用于设定有效时窗,在有效时窗内对比时移信号及正演信号,建立目标函数;
品质因子最优值及剩余相对时移最优值获取模块16,用于从品质因子初始值及剩余相对时移初始值开始,利用优化算法搜索目标函数的最小值,根据目标函数的最小值获取待测岩石样品的品质因子最优值及剩余相对时移最优值;
超声波速度计算模块17,根据初始起跳时间以及剩余相对时移最优值计算获得待测岩石样品的走时,并根据待测岩石样品的走时计算获得待测岩石样品的超声波速度。
在本实施例中,图7为本发明一实施例的正演模块的具体结构示意图。如图7所示,正演模块14包括:
标准频谱获取单元141,用于对标准信号进行傅里叶变换得到其相应的标准频谱;
正演信号频谱获取单元142,其中利用以下公式计算获得正演信号的频谱:
其中,S4(f)为正演信号的频谱;
G为待测岩石样品相对于参考样品的透射校正系数,利用公式(2)求得;
f为超声波脉冲透射测量系统的测量频率,为已知量;
T0为超声波测量系统的延迟走时,为已知量;
T1为标准起跳时间;
T2为初始起跳时间;
Δτ0为待测岩石样品的剩余相对时移初始值;
Q0为待测岩石样品的品质因子初始值;
Q1为参考样品的品质因子,为已知量;
S1(f)为标准频谱;
待测岩石样品相对于参考样品的透射校正系数是利用以下公式所得:
其中,ρ1为参考样品的密度,为已知量;
ρ2为待测岩石样品的密度,为已知量;
Z0为超声波脉冲透射测量系统的探头阻抗,为已知量;
V1参考样品的超声波速度,为已知量;
L2为待测岩石样品的长度,为已知量;
正演信号获取单元143,对正演信号的频谱进行傅里叶反变换得到正演信号。
在本实施例中,目标函数建立模块15中,建立目标函数利用的公式如下:
其中,Y为目标函数;
n为有效时窗内的采样点个数;
S3(ti)为时间域时移信号中每个采样点对应的幅值;
S4(ti)为时间域正演信号中每个采样点对应的幅值。
在本实施例中,品质因子最优值及剩余相对时移最优值获取模块16中利用的优化算法包括:模式搜索法、Rosenbrock方法、单纯形搜索法、Powell方法、最速下降法、牛顿法、共轭梯度法、拟牛顿法、信赖域方法、最小二次法、可行方向法、遗传算法、蚂蚁算法、模拟退火法以及神经网络法的至少其中之一。
在本实施例中,超声波速度计算模块17中,计算待测岩石样品的走时利用的公式如下:
其中,为待测岩石样品的走时;
Δτ2为步骤6得到的待测岩石样品的剩余相对时移最优值。
计算待测岩石样品的超声波速度利用的公式如下:
其中,V2为待测岩石样品的超声波速度;
为公式(4)得到的待测岩石样品的走时数据;
L2为待测岩石样品的长度,为已知量。
本发明提出的岩石超声波速度及品质因子同时测量方法及系统,与现有技术对这两个量单独进行测量的方法不同,本发明在时间域有效信号时窗范围内进行分析,能够同时测量得到岩石样品的超声波速度及品质因子,克服了两者之间的相互作用,还排除了后续干扰波的影响,可有效地提高测量精度,同时大大提高了实验测量的效率。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (14)
1.一种岩石超声波速度及品质因子同时测量方法,其特征在于,包括:
步骤1,选取一参考样品,利用超声波脉冲透射测量系统获取所述参考样品的标准信号,并拾取所述标准信号的标准起跳时间;
步骤2,选取一待测岩石样品,利用所述超声波脉冲透射测量系统获取所述待测岩石样品的测量信号,并与所述标准信号进行相关性分析,根据所述标准起跳时间确定所述测量信号的初始起跳时间;
步骤3,按照所述标准起跳时间与所述初始起跳时间的时差对所述测量信号进行时移,获取时移信号;
步骤4,设定所述待测岩石样品的品质因子初始值及剩余相对时移初始值,将其作用到所述标准信号进行正演,得到正演信号;
步骤5,设定有效时窗,在所述有效时窗内对比所述时移信号及正演信号,建立目标函数;
步骤6,从所述品质因子初始值及剩余相对时移初始值开始,利用优化算法搜索所述目标函数的最小值,根据所述目标函数的最小值获取所述待测岩石样品的品质因子最优值及剩余相对时移最优值;
步骤7,根据所述初始起跳时间以及剩余相对时移最优值计算获得所述待测岩石样品的走时,并根据所述待测岩石样品的走时计算获得待测岩石样品的超声波速度。
2.根据权利要求1所述的岩石超声波速度及品质因子同时测量方法,其特征在于,所述超声波脉冲透射测量系统利用的超声波信号包括P波、SV波、SH波。
3.根据权利要求1所述的岩石超声波速度及品质因子同时测量方法,其特征在于,所述步骤4,设定所述待测岩石样品的品质因子初始值及剩余相对时移初始值,将其作用到所述标准信号进行正演,得到正演信号包括以下步骤:
步骤41,对所述标准信号进行傅里叶变换得到其相应的标准频谱;
步骤42,利用以下公式计算所述正演信号的频谱:
其中,S4(f)为所述正演信号的频谱;
G为所述待测岩石样品相对于所述参考样品的透射校正系数,利用公式(2)求得;
f为所述超声波脉冲透射测量系统的测量频率,为已知量;
T0为所述超声波测量系统的延迟走时,为已知量;
T1为所述标准起跳时间;
T2为所述初始起跳时间;
Δτ0为所述待测岩石样品的剩余相对时移初始值;
Q0为所述待测岩石样品的品质因子初始值;
Q1为所述参考样品的品质因子,为已知量;
S1(f)为所述标准频谱;
所述待测岩石样品相对于参考样品的透射校正系数是利用以下公式所得:
其中,ρ1为所述参考样品的密度,为已知量;
ρ2为所述待测岩石样品的密度,为已知量;
Z0为所述超声波脉冲透射测量系统的探头阻抗,为已知量;
V1所述参考样品的超声波速度,为已知量;
L2为所述待测岩石样品的长度,为已知量;
步骤43,对所述正演信号的频谱进行傅里叶反变换得到所述正演信号。
4.根据权利要求3所述的岩石超声波速度及品质因子同时测量方法,其特征在于,所述步骤5,建立目标函数利用的公式如下:
其中,Y为所述目标函数;
n为所述有效时窗内的采样点个数;
S3(ti)为时间域所述时移信号中每个采样点对应的幅值;
S4(ti)为时间域所述正演信号中每个采样点对应的幅值。
5.根据权利要求4所述的岩石超声波速度及品质因子同时测量方法,其特征在于,所述步骤6利用的优化算法包括:模式搜索法、Rosenbrock方法、单纯形搜索法、Powell方法、最速下降法、牛顿法、共轭梯度法、拟牛顿法、信赖域方法、最小二次法、可行方向法、遗传算法、蚂蚁算法、模拟退火法以及神经网络法的至少其中之一。
6.根据权利要求5所述的岩石超声波速度及品质因子同时测量方法,其特征在于,所述步骤7,计算所述待测岩石样品的走时利用的公式如下:
其中,为所述待测岩石样品的走时;
Δτ2为所述步骤6得到的所述待测岩石样品的剩余相对时移最优值。
7.根据权利要求6所述的岩石超声波速度及品质因子同时测量方法,其特征在于,所述步骤7,计算所述待测岩石样品的超声波速度利用的公式如下:
其中,V2为所述待测岩石样品的超声波速度;
为公式(4)得到的待测岩石样品的走时数据;
L2为所述待测岩石样品的长度,为已知量。
8.一种岩石超声波速度及品质因子同时测量系统,其特征在于,包括:
标准信号获取模块,用于根据一参考样品,通过超声波脉冲透射测量系统获取所述参考样品的标准信号,并拾取所述标准信号的标准起跳时间;
测量信号获取模块,用于根据一待测岩石样品,通过所述超声波脉冲透射测量系统获取所述待测岩石样品的测量信号,并与所述标准信号进行相关性分析,根据所述标准起跳时间确定所述测量信号的初始起跳时间;
时移模块,用于按照所述标准起跳时间与所述初始起跳时间的时差对所述测量信号进行时移,获取时移信号;
正演模块,用于设定所述待测岩石样品的品质因子初始值及剩余相对时移初始值,将其作用到所述标准信号进行正演,得到正演信号;
目标函数建立模块,用于设定有效时窗,在所述有效时窗内对比所述时移信号及正演信号,建立目标函数;
品质因子最优值及剩余相对时移最优值获取模块,用于从所述品质因子初始值及剩余相对时移初始值开始,利用优化算法搜索所述目标函数的最小值,根据所述目标函数的最小值获取所述待测岩石样品的品质因子最优值及剩余相对时移最优值;
超声波速度计算模块,根据所述初始起跳时间以及剩余相对时移最优值计算获得所述待测岩石样品的走时,并根据所述待测岩石样品的走时计算获得待测岩石样品的超声波速度。
9.根据权利要求8所述的岩石超声波速度及品质因子同时测量系统,其特征在于,所述超声波脉冲透射测量系统利用的超声波信号包括P波、SV波、SH波。
10.根据权利要求8所述的岩石超声波速度及品质因子同时测量系统,其特征在于,所述正演模块包括:
标准频谱获取单元,用于对所述标准信号进行傅里叶变换得到其相应的标准频谱;
正演信号频谱获取单元,其中利用以下公式计算获得所述正演信号的频谱:
其中,S4(f)为所述正演信号的频谱;
G为所述待测岩石样品相对于所述参考样品的透射校正系数,利用公式(2)求得;
f为所述超声波脉冲透射测量系统的测量频率,为已知量;
T0为所述超声波测量系统的延迟走时,为已知量;
T1为所述标准起跳时间;
T2为所述初始起跳时间;
Δτ0为所述待测岩石样品的剩余相对时移初始值;
Q0为所述待测岩石样品的品质因子初始值;
Q1为所述参考样品的品质因子,为已知量;
S1(f)为所述标准频谱;
所述待测岩石样品相对于参考样品的透射校正系数是利用以下公式所得:
其中,ρ1为所述参考样品的密度,为已知量;
ρ2为所述待测岩石样品的密度,为已知量;
Z0为所述超声波脉冲透射测量系统的探头阻抗,为已知量;
V1所述参考样品的超声波速度,为已知量;
L2为所述待测岩石样品的长度,为已知量;
正演信号获取单元,对所述正演信号的频谱进行傅里叶反变换得到所述正演信号。
11.根据权利要求10所述的岩石超声波速度及品质因子同时测量系统,其特征在于,所述目标函数建立模块中,建立目标函数利用的公式如下:
其中,Y为所述目标函数;
n为所述有效时窗内的采样点个数;
S3(ti)为时间域所述时移信号中每个采样点对应的幅值;
S4(ti)为时间域所述正演信号中每个采样点对应的幅值。
12.根据权利要求11所述的岩石超声波速度及品质因子同时测量系统,其特征在于,所述品质因子最优值及剩余相对时移最优值获取模块中利用的优化算法包括:模式搜索法、Rosenbrock方法、单纯形搜索法、Powell方法、最速下降法、牛顿法、共轭梯度法、拟牛顿法、信赖域方法、最小二次法、可行方向法、遗传算法、蚂蚁算法、模拟退火法以及神经网络法的至少其中之一。
13.根据权利要求12所述的岩石超声波速度及品质因子同时测量系统,其特征在于,所述超声波速度计算模块中,计算所述待测岩石样品的走时利用的公式如下:
其中,为所述待测岩石样品的走时;
Δτ2为所述步骤6得到的所述待测岩石样品的剩余相对时移最优值。
14.根据权利要求13所述的岩石超声波速度及品质因子同时测量系统,其特征在于,所述超声波速度计算模块中,计算所述待测岩石样品的超声波速度利用的公式如下:
其中,V2为所述待测岩石样品的超声波速度;
为公式(4)得到的待测岩石样品的走时数据;
L2为所述待测岩石样品的长度,为已知量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410446685.2A CN104237383B (zh) | 2014-09-03 | 2014-09-03 | 一种岩石超声波速度及品质因子同时测量方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410446685.2A CN104237383B (zh) | 2014-09-03 | 2014-09-03 | 一种岩石超声波速度及品质因子同时测量方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104237383A true CN104237383A (zh) | 2014-12-24 |
CN104237383B CN104237383B (zh) | 2016-08-17 |
Family
ID=52225878
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410446685.2A Active CN104237383B (zh) | 2014-09-03 | 2014-09-03 | 一种岩石超声波速度及品质因子同时测量方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104237383B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105572229A (zh) * | 2014-10-15 | 2016-05-11 | 空中客车运营简化股份公司 | 部件无损测试系统的校准的验证方法及组装件 |
CN110245462A (zh) * | 2019-07-09 | 2019-09-17 | 西南石油大学 | 一种碳酸盐岩孔隙度预测方法及预测模型构建方法 |
CN111530845A (zh) * | 2020-05-25 | 2020-08-14 | 重庆大学 | 一种基于超声波的均压电极手持式除垢装置及除垢方法 |
CN113884575A (zh) * | 2021-08-18 | 2022-01-04 | 山东省交通规划设计院集团有限公司 | 一种围岩分级方法及系统 |
CN114861563A (zh) * | 2022-04-27 | 2022-08-05 | 中国石油大学(华东) | 物理嵌入深度学习地层压力预测方法、装置、介质及设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102023311A (zh) * | 2010-08-10 | 2011-04-20 | 中国石油大学(华东) | 地层的品质因子谱及其求取方法 |
CN102778695A (zh) * | 2012-07-23 | 2012-11-14 | 中国科学院地质与地球物理研究所 | 一种获取煤岩品质因子的控制方法 |
CN103412336A (zh) * | 2013-07-22 | 2013-11-27 | 中国石油天然气股份有限公司 | 一种非均质油藏中岩石系统的纵波速度预测方法 |
-
2014
- 2014-09-03 CN CN201410446685.2A patent/CN104237383B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102023311A (zh) * | 2010-08-10 | 2011-04-20 | 中国石油大学(华东) | 地层的品质因子谱及其求取方法 |
CN102778695A (zh) * | 2012-07-23 | 2012-11-14 | 中国科学院地质与地球物理研究所 | 一种获取煤岩品质因子的控制方法 |
CN103412336A (zh) * | 2013-07-22 | 2013-11-27 | 中国石油天然气股份有限公司 | 一种非均质油藏中岩石系统的纵波速度预测方法 |
Non-Patent Citations (5)
Title |
---|
D.H.JOHNSTON ET AL.: "Attenuation of seismic waves in dry and saturated rocks:Ⅱ.Mechanisms", 《GEOPHYSICS》 * |
KHAZANEHDARI J. ET AL.: "Acoustic and petrophysical relationships in low-shale sandstone reservoir rocks", 《GEOPHYSICAL PROSPECTING》 * |
张建利等: "横波偏振方法在煤样超声测量实验中的应用", 《煤炭学报》 * |
李生杰等: "准噶尔盆地岩石品质因子与速度分析", 《内陆地震》 * |
魏建新等: "岩块超声横波速度高精度测试方法试验研究", 《中国石油大学学报(自然科学版)》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105572229A (zh) * | 2014-10-15 | 2016-05-11 | 空中客车运营简化股份公司 | 部件无损测试系统的校准的验证方法及组装件 |
CN110245462A (zh) * | 2019-07-09 | 2019-09-17 | 西南石油大学 | 一种碳酸盐岩孔隙度预测方法及预测模型构建方法 |
CN111530845A (zh) * | 2020-05-25 | 2020-08-14 | 重庆大学 | 一种基于超声波的均压电极手持式除垢装置及除垢方法 |
CN111530845B (zh) * | 2020-05-25 | 2022-05-31 | 重庆大学 | 一种基于超声波的均压电极手持式除垢装置及除垢方法 |
CN113884575A (zh) * | 2021-08-18 | 2022-01-04 | 山东省交通规划设计院集团有限公司 | 一种围岩分级方法及系统 |
CN113884575B (zh) * | 2021-08-18 | 2023-08-01 | 山东省交通规划设计院集团有限公司 | 一种围岩分级方法及系统 |
CN114861563A (zh) * | 2022-04-27 | 2022-08-05 | 中国石油大学(华东) | 物理嵌入深度学习地层压力预测方法、装置、介质及设备 |
CN114861563B (zh) * | 2022-04-27 | 2022-12-13 | 中国石油大学(华东) | 物理嵌入深度学习地层压力预测方法、装置、介质及设备 |
Also Published As
Publication number | Publication date |
---|---|
CN104237383B (zh) | 2016-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109425896B (zh) | 白云岩油气储层分布预测方法及装置 | |
CN102692645B (zh) | 利用纵波、转换波数据联合反演储层纵横波速度比的方法 | |
Tang et al. | Fast inversion of formation permeability from Stoneley wave logs using a simplified Biot-Rosenbaum model | |
CN104237383B (zh) | 一种岩石超声波速度及品质因子同时测量方法及系统 | |
NO341729B1 (no) | Bestemmelse av seismisk anisotropi omfattende analyse av signaldata reflektert fra en laggrense i undergrunnen | |
CN107092030B (zh) | 通过瑞雷波探测软黏土包的方法 | |
CN104570108A (zh) | 估算等效品质因子方法及用其估算地层品质因子的方法 | |
CN103852785B (zh) | 地层各向异性的评价方法 | |
CN103257363A (zh) | 一种探测地下裂缝性储层中裂缝倾角的方法 | |
CN107544087B (zh) | 一种测量近地表地层品质因子的方法及装置 | |
CN105242311B (zh) | 预估初至拾取的时窗的方法和装置 | |
CN104570103B (zh) | 一种低信噪比地震资料的井约束速度谱拾取方法 | |
Liu et al. | Detection of karst cavity beneath cast-in-place pile using the instantaneous phase difference of two receiver recordings | |
CN102830170A (zh) | 一种基于超声测试获取煤样横波信号的控制方法及装置 | |
CN106770668A (zh) | 一种用于单孔的基桩质量声波透射法检测方法 | |
US4933911A (en) | Method for determining seismic velocities | |
CN108693560A (zh) | 一种基于互相关道的散射波成像方法及系统 | |
Huang et al. | Fast-forward modeling of compressional arrival slowness logs in high-angle and horizontal wells | |
Holliger et al. | Attenuation of broad‐band (50–1500 Hz) seismic waves in granitic rocks near the Earth’surface | |
RU2375729C1 (ru) | Геофизический радиолокатор | |
Callerio et al. | Cross-hole tests at Zelazny most tailings pond, Poland-highlights and statistical interpretation of results | |
CN105445795A (zh) | 倾斜界面的速度分析方法及其系统 | |
CN112305595B (zh) | 基于折射波分析地质体结构的方法及存储介质 | |
CN104123440B (zh) | 一种在振幅不匹配时使时间域地震波形反演收敛的方法 | |
CN103576189B (zh) | 一种随钻地震直达波组搜索方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |