CN107389794B - 测量岩石衰减系数的方法与系统 - Google Patents
测量岩石衰减系数的方法与系统 Download PDFInfo
- Publication number
- CN107389794B CN107389794B CN201710691134.6A CN201710691134A CN107389794B CN 107389794 B CN107389794 B CN 107389794B CN 201710691134 A CN201710691134 A CN 201710691134A CN 107389794 B CN107389794 B CN 107389794B
- Authority
- CN
- China
- Prior art keywords
- sample
- rock
- quality factor
- resonance
- measuring
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/041—Analysing solids on the surface of the material, e.g. using Lamb, Rayleigh or shear waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/01—Indexing codes associated with the measuring variable
- G01N2291/014—Resonance or resonant frequency
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/023—Solids
- G01N2291/0232—Glass, ceramics, concrete or stone
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/10—Number of transducers
- G01N2291/102—Number of transducers one emitter, one receiver
Abstract
本申请实施例提供了一种测量岩石衰减系数的方法与系统,该系统包括敞口容器,用于容纳用作传声介质的液体;震源,用于产生驱使系统发生共振的声波信号;浸没于液体中的共振腔,用于输出系统在共振时产生的共振声压场信号;共振声压场信号包括共振腔内未加载任何样品、单独加载第一样品、单独加载第二样品、及单独加载岩石样品下的共振声压场信号;置于共振腔内的水听器,用于采集共振声压场信号;锁相放大器,用于向震源输出震源激励信号,并将震源激励信号及水听器采集的共振声压场信号进行锁相处理后提供给上位机;上位机,用于根据锁相放大器提供的信号获得岩石样品的品质因子。本申请实施例可以提高测量岩石衰减系数的实用性和测量效率。
Description
技术领域
本申请涉及储层岩石物理性质研究技术领域,尤其是涉及一种测量岩石衰减系数的方法与系统。
背景技术
岩石物理性质的实验室测量是地球物理应用研究中的一项基础性研究工作。由于地下岩石的复杂性,一般认为,岩石的岩石物理性质是随频率而变化的,在实验室中准确测定现场应用频率段内的岩石物理性质对利用地震勘探和声波测井资料来认识地下岩石的性质,研究地震相和测井相随频率的变化特征等都具有重要的意义。而岩石中地震波衰减变化及其产生机理是其中被广泛关注研究的对象之一。
目前,实验室条件下主要存在着三种测量岩石衰减系数的方法:1)脉冲透射法,2)共振棒法,3)应力应变法。脉冲透射方法使用相同尺寸的岩石样品及参考样品(通常为铝),基于激发脉冲透过两样品的振幅之比求取岩石样品的衰减系数;其测量频段为几百千赫至兆赫,与地震勘探频段差异较大,二者可能存在着不同的衰减机理,测量结果不可直接用于解释地震数据。共振棒方法下,圆柱形或平行六面体形样品受迫运动进入或穿过共振状态,利用频谱扫描方法获取共振峰值fmax及半功率带宽Δf,然后据此确定品质因子;其测量频段为千赫,接近声波测井频段;当测量频率较低时,需较长长度、耐用且均匀的样品,然而这样的样品不容易获得,且样品的准备耗时较长。在应力应变方法下,周期性加力使物体产生形变,形变变化速率低于样品的自然共振频率,样品的衰减由施加的应力及相应应变所得相位差异计算;测量频段更低,为几赫至几千赫,包含地震勘探频段,在较宽的低频带范围内测量出相对准确且连续的衰减参数,但准备岩石样品耗时较长。
因此,目前亟需一种新的测量岩石衰减系数得技术方案,以提高测量岩石衰减系数的实用性和测量效率。
发明内容
本申请实施例的目的在于提供一种测量岩石衰减系数的方法与系统,以提高测量岩石衰减系数的实用性和测量效率。
为达到上述目的,一方面,本申请实施例提供了一种测量岩石衰减系数的系统,包括:
敞口容器,用于容纳用作传声介质的液体;
震源,用于产生驱使所述系统发生共振的声波信号;
浸没于所述液体中的共振腔,用于输出所述系统在共振时产生的共振声压场信号;所述共振声压场信号包括所述共振腔内未加载任何样品、单独加载第一样品、单独加载第二样品、及单独加载岩石样品下的共振声压场信号;
置于所述共振腔内的水听器,用于采集所述共振声压场信号;
锁相放大器,用于向所述震源输出震源激励信号,并将所述震源激励信号及所述水听器采集的共振声压场信号进行锁相处理后提供给上位机;
上位机,用于根据所述锁相放大器提供的信号获得所述岩石样品的品质因子。
本申请实施例的测量岩石衰减系数的系统,所述系统还包括样品移动测量装置,所述样品移动测量装置包括:
悬挂样品的传动机构,用于带动样品相对于所述共振腔上下移动;
与所述传动机构相连的步进电机,用于驱动所述传动机构;
位置控制器,用于在所述上位机的控制下向所述步进电机输出控制信号。
本申请实施例的测量岩石衰减系数的系统,所述共振腔的上端及下端敞口。
本申请实施例的测量岩石衰减系数的系统,所述系统还包括:
功率放大器,用于对所述锁相放大器输出的震源激励信号进行功率放大,并向所述震源输出功率放大后的震源激励信号;
预放大器,用于对所述水听器采集的共振声压场信号进行前置放大,并向所述锁相放大器输出前置放大后的共振声压场信号。
本申请实施例的测量岩石衰减系数的系统,所述共振腔包括圆筒形铝筒。
本申请实施例的测量岩石衰减系数的系统,所述敞口容器包括圆柱形玻璃桶。
本申请实施例的测量岩石衰减系数的系统,所述液体包括硅油。
本申请实施例的测量岩石衰减系数的系统,所述震源包括压电陶瓷片,其将所述震源激励信号转化为机械振动并以声波信号形式发出。
本申请实施例的测量岩石衰减系数的系统,所述传动机构包括滑轮组。
本申请实施例的测量岩石衰减系数的系统,所述震源位于所述敞口容器底部中心。
另一方面,本申请实施例还提供了一种利用上述测量岩石衰减系数的系统测量岩石衰减系数的方法,包括以下步骤:
确定岩石样品、第一样品及第二样品的体积;所述第一样品及所述第二样品的弹性性质已知;
确定所述共振腔内未加载任何样品、单独加载所述第一样品及单独加载所述第二样品下所述系统的共振频率;
根据对应条件下所述系统的共振频率,确定对应条件下所述系统的品质因子;
根据对应条件下所述系统的品质因子及所述第一样品及第二样品的体积,确定所述液体的品质因子及校正参数;
确定所述岩石样品的压缩系数,并确定所述共振腔内单独加载所述岩石样品下所述系统的品质因子;
根据所述液体的品质因子、所述校正参数、所述岩石样品的体积、所述共振腔内未加载任何样品下所述系统的品质因子、及所述共振腔内单独加载所述岩石样品下所述系统的品质因子,确定所述岩石样品的品质因子。
本申请实施例的测量岩石衰减系数的方法,所述根据对应条件下所述系统的共振频率,确定对应条件下所述系统的品质因子,包括:
将对应条件下所述系统的共振频率除以对应条件下的半功率带宽,得到对应条件下所述系统的品质因子。
本申请实施例的测量岩石衰减系数的方法,根据以下公式确定所述液体的品质因子及校正参数:
其中,为第一样品的品质因子,为液体的品质因子,C为校正参数,VAl为第一样品的体积,QS-Al为共振腔内单独加载第一样品下系统的品质因子,Q0为共振腔内未加载任何样品下系统的品质因子,κR0为液体的压缩系数,κAl为第一样品的压缩系数,为第二样品的品质因子,VLu为第二样品的体积,QS-Lu为共振腔内单独加载第二样品下系统的品质因子,κLu为第二样品的压缩系数。
本申请实施例的测量岩石衰减系数的方法,根据以下公式确定所述岩石样品的压缩系数:
其中,ωs为共振腔内单独加载岩石样品下系统的共振角频率,ω0为共振腔内未加载任何样品下系统的共振角频率,C为校正参数,Vs为岩石样品的体积,κs为岩石样品的压缩系数,κ0为液体的压缩系数。
本申请实施例的测量岩石衰减系数的方法,所述确定所述共振腔内单独加载所述岩石样品下所述系统的品质因子,包括:
确定所述共振腔内单独加载所述岩石样品下所述系统的共振频率;
将所述共振腔内单独加载所述岩石样品下所述系统的共振频率,除以所述共振腔内单独加载所述岩石样品下所述系统的半功率带宽,得到对应条件下所述系统的品质因子。
本申请实施例的测量岩石衰减系数的方法,根据以下公式确定所述岩石样品的品质因子:
其中,为所述岩石样品的品质因子,为液体的品质因子,C为校正参数,Vs为岩石样品的体积,QS为共振腔内单独加载岩石样品下系统的品质因子,Q0为共振腔内未加载任何样品下系统的品质因子,κR0为液体的压缩系数,κRs为岩石样品的压缩系数。
由以上本申请实施例提供的技术方案可见,本申请实施的系统包括敞口容器,用于容纳用作传声介质的液体;震源,用于产生驱使系统发生共振的声波信号;浸没于液体中的共振腔,用于输出系统在共振时产生的共振声压场信号;共振声压场信号包括共振腔内未加载任何样品、单独加载第一样品、单独加载第二样品、及单独加载岩石样品下的共振声压场信号;置于共振腔内的水听器,用于采集共振声压场信号;锁相放大器,用于向震源输出震源激励信号,并将震源激励信号及水听器采集的共振声压场信号进行锁相处理后提供给上位机;上位机,用于根据锁相放大器提供的信号获得岩石样品的品质因子,从而实现了岩石衰减系数的测量。由于本申请实施例对于岩石样品没有特定要求。因而本申请实施例提高了测量岩石衰减系数的实用性和测量效率。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1为本申请一实施例中共振腔内的压力场分布示意图;
图2为本申请一实施例中测量岩石衰减系数的系统的结构示意图;
图3为本申请一实施例中测量岩石衰减系数的方法的流程图;
图4为本申请一实施例在共振腔内未加载任何样品(即共振腔空腔)、单独加载标准铝样品、及单独加载标准玻璃样品下的共振曲线示意图;
图5为本申请一实施例中在共振腔内未加载任何样品(即共振腔空腔)、及单独加载待测量岩石样品下的共振曲线示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
参考图1所示,本申请实施例的测量岩石衰减系数的系统可以包括敞口容器1、共振腔2、震源3、水听器4、锁相放大器和上位机等。
在本申请一些实施例中,敞口容器1用于容纳用作传声介质的液体,因此只要其能够容纳用作传声介质的液体(例如硅油等),且其上端敞口以便于放入共振腔2、震源3和水听器4等,敞口容器1的具体形状和尺寸没有特定的要求。例如在本申请一示例性实施例中,敞口容器1可以为如图1所示的圆柱形有机玻璃桶,该圆柱形有机玻璃桶的直径可以为0.4米,高可以为1米,圆柱形有机玻璃桶内的硅油页面高度可以为0.9米。
在本申请一些实施例中,震源3用于产生驱使系统发生共振的声波信号。即震源3可以为将震源激励信号转化为机械振动并以声波信号形式发出的部件或装置等。在本申请一示例性实施例中,所述震源3可安装于所述敞口容器1底部中心。所述震源3例如可以为压电陶瓷片等。
在本申请一些实施例中,共振腔2浸没于所述液体中,以用于输出所述系统在共振时产生的共振声压场信号。其中,所述共振声压场信号可以包括所述共振腔2内未加载任何样品、单独加载第一样品、单独加载第二样品、及单独加载岩石样品5下的共振声压场信号。在本申请一示例性实施例中,所述共振腔2可以为上端及下端敞口的结构;例如可以为如图1所示的圆筒形铝筒,该圆筒形铝筒的壁厚0.01米,高0.61米。其中,第一样品和第二样品为弹性性质已知的标准样品(例如下文提及的标准率铝样品Al和标准有机玻璃样品Lu)。所述标准样品,即所述第一样品和第二样品的形状、大小和体积可均与岩石样品5的形状、大小和体积相同,这样可以有利于获得更为精确的岩石衰减系数。在本申请一示例性实施例中,共振腔2可如图1所示的那样悬挂于敞口容器1中。
在本申请一些实施例中,水听器4可以置于所述共振腔2内,以用于采集所述系统在共振时产生的共振声压场信号。需要说明的是,水听器4并不构成对本申请的限制,任何合适的可用于在液体中采集所述系统在共振时产生的共振声压场信号的装置或部件,都可应用于本申请。
锁相放大器,用于向所述震源3输出震源激励信号,并将所述震源激励信号及所述水听器4采集的共振声压场信号进行锁相处理后提供给上位机。在本申请一示例性实施例中,所述锁相放大器例如可以是从低频到高频每隔一定频率发射一个正弦波信号,以用作震源激励信号。相应的,震源3在对应正弦波信号的作用下可发出对应的声波信号。如此,从低频到高频总可以确定系统共振频率。需要说明的是,即便发射的是单频率的正弦波信号,比如发射1000Hz正弦波后系统发生共振,但在整个系统发生共振时波是有很多频率成分的,当然主要成分是1000Hz的,但是混杂着无数别的频率成分。因而,水听器4采集到的共振声压场信号里除了1000Hz正弦波外,也是含有众多其他频率成分的波。此时,利用锁相放大器则可以锁住希望监测驱使系统发生共振的频率信号,其他的频率成分则滤除掉。
上位机,用于根据所述锁相放大器提供的信号获得所述岩石样品5的品质因子。本申请实施例中,所述岩石样品5的品质因子即为所述岩石样品5的衰减系数,下文方法步骤中将具体说明上位机的处理过程。
为了进一步提高测量效率,所述系统还可以包括样品移动测量装置,所述样品移动测量装置可包括传动机构6、步进电机和位置控制器。
在本申请一些实施例中,传动机构6可悬挂样品(例如上文提到的标准率铝样品Al、标准有机玻璃样品Lu和岩石样品),用于带动样品相对于所述共振腔2上下移动。在本申请一示例性实施例中,传动机构6可以为如图1所示的滑轮组。当然,在本申请其他一些实施例中,根据需要也可以采用其他适合的传动机构。步进电机与所述传动机构6相连,以用于驱动所述传动机构6。在本申请其他一些实施例中,根据需要,所述步进电机也可以采用别的诸如其他可控电机等动力源。位置控制器,用于在所述上位机的控制下向所述步进电机输出控制信号。如图1所示,在本申请一示例性实施例中,所述上位机和位置控制器之间例如可以通过通用接口总线(General-Purpose Inter face Bus,简称GPIB)进行通信。
此外,在本申请其他一些实施例中,还可以利用上述样品移动测量装置获得不同位置下的岩石样品5的衰减系数,从而通过对不同位置下的岩石样品5的衰减系数进行取平均值计算,可以获得更为精确的岩石样品5的衰减系数。
如图1所示,在本申请一些实施例中,所述系统还可以包括功率放大器和预放大器。其中,功率放大器可用于对所述锁相放大器输出的震源激励信号进行功率放大,并向所述震源输出功率放大后的震源激励信号。预放大器(也称为前置放大器)可用于对所述水听器4采集的共振声压场信号进行前置放大,并向所述锁相放大器输出前置放大后的共振声压场信号。
如图2所示,当共振发生时,将产生沿着共振腔2轴向空间谐变的声压场。共振腔2中轴向一阶共振频率由如下公式决定:
其中,c0是充满共振腔2的硅油的声速,L是共振腔2的长度。这样,共振腔2中沿着轴向方向,声压曲线呈正弦分部,声速正比于声压的空间导数。相应地,在共振腔2中有两个特别的位置:共振腔2的中心与两端,分别称之为声压节点与速度节点。对于一阶共振频率,声压节点处声速为最大,而速度节点处声压为最大。在如图1所示的共振腔2中,待测的岩石样品5引入到充满硅油的共振腔2中产生声扰动。假设p1与p2分别为岩石样品5引入共振腔2前、后共振腔2中声压的分布函数,相应地,由于岩石样品5产生的声扰动,共振腔2的共振频率从ω1移动到ω2,用声波方程描述共振频率的一阶变换如下:
(1)在共振腔2的边界上有▽p=0;
(2)在共振腔2的两个开口端有p=0。
使用这样的两个边界条件可以解得方程(2)与方程(3)的声压解。在方程(2)的两端乘以p2,在方程(3)的两端乘以p1并对整个共振腔2进行体积积分(VC),可以得到
进一步操作方程(6)以及方程(7),最后可得到扰动方程(Perturbationequation)
在这些方程中,k1=ω1/c1是波数;VS与VC分别是共振腔2与岩石样品5的体积;通常而言方程(8)中的系数A与B可以表示成<p>2与<ρ1c1v>2,其物理意义分别为平均声压与共振腔2中平均粒子振动速度。可分别使用ω0与ωs(或f0与fs)表示岩石样品5引入共振腔2前后系统的共振频率,κ0与κs分别表示硅油与岩石样品5的压缩系数,而ρ0与ρs分别表示共振腔2中的硅油与岩石样品5的密度。这样,可以重写方程(8)为如下的方程
或
方程(9)即为扰动方程。
可以注意到,方程(9)说明了岩石样品5与共振腔2中硅油之间的压缩系数的对比((κs-κ0)/κ0)以及密度的对比((ρs-ρ0)/ρ0)直接贡献了样品引入共振腔2前后共振频率的偏移。通常而言,一块岩石样品5的密度是容易测得的,因此,岩石样品5的压缩系数是该系统测得的主要参数。如果一块岩石样品5放在共振腔2的速度节点或声压反节点上,在这个节点上声速为0,这样,方程(9)的第二项便可被消去了。方程(9)便可被简化为如下的表达式:
进一步获得方程(11)的紧凑形式:
其中,为校正参数。
发生共振时,共振曲线如图5所示。两条曲线分别对应空腔(即共振腔2内未加载人样样品)与引入样本(即引入待测的岩石样品5)后共振腔2的共振曲线。此时,系统的品质因子Q=fmax/Δf是包括样品在内的整个共振腔2系统的品质因子,其中,fmax为共振峰值;而Δf为半功率带宽(共振振幅峰值f的处),为了获得样品本身的衰减特性,本申请提出复角频率及复模量概念,以研究共振腔2的品质因子Q与待测的岩石样品5的品质因子Q*间的关系,最终求取待测的岩石样品5的品质因子Q*。
研究表明,共振腔2品质因子Q与复角频率ω的关系式,如下式所示:
ω=ωR+iωI (13)
在发生共振时,在速度节点(对应共振中心点)处的共振角频率与压缩系数间的关系,如上述公式(12)所示,即
(i.e.δω=CVsδκ)
其中,ω0和ωs分别是空腔与有样品(即加载了待测的岩石样品5)时的共振角频率,κ0和κs分别是腔内硅油与岩石样品5的压缩系数,Vs为岩石样品5的体积,C为校正参数,δω为共振角频率相对变化量,δκ为压缩系数相对变化量。为了得到共振腔2的品质因子Q与岩石样品5的品质因子Q*间的关系,首先需将复角频率和复压缩系数引入公式(12),并考虑ωI/ωR<<1以及ωRs≈ωR0:
其中,QS=ωRs/(2ωIs)=fs/Ws是当共振腔2内加载了一个岩石样品5时测量得到共振时系统的品质因子,而Q0=ωR0/(2ωI0)=f0/W0则是共振腔2内为空腔时测量得到共振时系统的品质因子。
将体积模量用复数形式表示,则有K=KR+iKI,此时用体积模量定义的样品的品质因子为:Q*=KR/KI。由于压缩系数是体积模量的倒数,因此压缩系数可以写成κ=κR-iκI,以其定义的品质因子写成Q*=-κR/κI。进而,公式(12)中等号右边的参数δκ可以表示为:
由于公式(12)等号左右的虚部应该相等,因此结合(12)、(15)和(16),可以得到:
根据品质因子在压缩系数下的定义Q*=-κR/κI,共振腔2内流体的品质因子定义为将其带入表达式(18)得到:
于是,岩石样品5压缩系数的虚部κIS为:
其中,是腔内流体的品质因子,κR0是腔内流体的复压缩系数实部(即硅油的复压缩系数实部),κRs是可由公式(12)中得到的样品复压缩系数实部。
在公式(21)中,岩石样品5的体积为VS、硅油的压缩系数κR0为已知参数;通过公式(12)可反演计算出岩石样品5的压缩系数κRs;参数Q0和QS可以通过在共振中心点位置处的扫频曲线,按照公式Q=fmax/Δf计算得到。因此,未知参数包括:硅油自身的品质因子校正参数C和待求的岩石样品5的品质因子为了求得首先需要得到硅油的品质因子和校正参数C。
其中,和分别为标准率铝样品Al和标准有机玻璃样品Lu自身的品质因子;Q0、QS-Al和QS-Lu分别为空腔、单独加载了标准率铝样品Al和单独加载了标准有机玻璃样品Lu时系统的品质因子(空腔、单独加载了标准率铝样品Al和单独加载了标准有机玻璃样品Lu时系统的共振曲线如图4所示),可以根据品质因子定义式Q=fmax/Δf计算得到;κAl和κLu为标准铝样品Al和标准有机玻璃样品Lu的压缩系数,也可以通过查表获得。因此,只要获得和即可联立公式(22)和公式(23),求解出硅油的品质因子和校正参数C。在实际操作中,还可以用谱比法获得和由于标准率铝样品Al以及标准有机玻璃Lu的品质因子无频散(即不随频率变化),因此可采用谱比法获得的和为kHz(本测量系统的测试频段)数量级的品质因子。
据此,参考图3所示,利用本申请上述实施例的测量岩石衰减系数的系统测量岩石衰减系数的方法,可以包括以下步骤:
S101、确定岩石样品、第一样品及第二样品的体积;所述第一样品及所述第二样品的弹性性质已知。
S102、确定所述共振腔内未加载任何样品、单独加载所述第一样品及单独加载所述第二样品下所述系统的共振频率。
S103、根据对应条件下所述系统的共振频率,确定对应条件下所述系统的品质因子。具体的,可将对应条件下所述系统的共振频率除以对应条件下的半功率带宽,得到对应条件下所述系统的品质因子。
S104、根据对应条件下所述系统的品质因子及所述第一样品及第二样品的体积,确定所述液体的品质因子及校正参数。具体实现可参见上述公式(22)和(23)。
S105、确定所述岩石样品的压缩系数,并确定所述共振腔内单独加载所述岩石样品下所述系统的品质因子。其中,确定所述岩石样品的压缩系数的具体实现可参见上述公式(12)。
S106、根据所述液体的品质因子、所述校正参数、所述岩石样品5的体积、所述共振腔2内未加载任何样品下所述系统的品质因子、及所述共振腔2内单独加载所述岩石样品5下所述系统的品质因子,确定所述岩石样品5的品质因子。具体实现可参见上述公式(21)。
虽然上文描述的过程流程包括以特定顺序出现的多个操作,但是,应当清楚了解,这些过程可以包括更多或更少的操作,这些操作可以顺序执行或并行执行(例如使用并行处理器或多线程环境)。
还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的方法或者系统不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种的方法或者系统所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的的方法或者系统中还存在另外的相同要素。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于方法实施例而言,由于其基本相似于系统实施例,所以描述的比较简单,相关之处参见系统实施例的部分说明即可。
以上所述仅为本申请的实施例而已,并不用于限制本申请。对于本领域技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本申请的权利要求范围之内。
Claims (14)
1.一种利用测量岩石衰减系数的系统测量岩石衰减系数的方法,其特征在于,所述系统包括:
敞口容器,用于容纳用作传声介质的液体;
震源,用于产生驱使所述系统发生共振的声波信号;
浸没于所述液体中的共振腔,用于输出所述系统在共振时产生的共振声压场信号;所述共振声压场信号包括所述共振腔内未加载任何样品、单独加载第一样品、单独加载第二样品、及单独加载岩石样品下的共振声压场信号;
置于所述共振腔内的水听器,用于采集所述共振声压场信号;
锁相放大器,用于向所述震源输出震源激励信号,并将所述震源激励信号及所述水听器采集的共振声压场信号进行锁相处理后提供给上位机;
上位机,用于根据所述锁相放大器提供的信号获得所述岩石样品的品质因子;
所述方法包括以下步骤:
确定岩石样品、第一样品及第二样品的体积;所述第一样品及所述第二样品的弹性性质已知;
确定所述共振腔内未加载任何样品、单独加载所述第一样品及单独加载所述第二样品下所述系统的共振频率;
根据对应条件下所述系统的共振频率,确定对应条件下所述系统的品质因子;
根据对应条件下所述系统的品质因子及所述第一样品及第二样品的体积,确定所述液体的品质因子及校正参数;
确定所述岩石样品的压缩系数,并确定所述共振腔内单独加载所述岩石样品下所述系统的品质因子;
根据所述液体的品质因子、所述校正参数、所述岩石样品的体积、所述共振腔内未加载任何样品下所述系统的品质因子、及所述共振腔内单独加载所述岩石样品下所述系统的品质因子,确定所述岩石样品的品质因子;其中,
根据以下公式确定所述岩石样品的品质因子:
2.如权利要求1所述的测量岩石衰减系数的方法,其特征在于,所述根据对应条件下所述系统的共振频率,确定对应条件下所述系统的品质因子,包括:
将对应条件下所述系统的共振频率除以对应条件下的半功率带宽,得到对应条件下所述系统的品质因子。
5.如权利要求1所述的测量岩石衰减系数的方法,其特征在于,所述确定所述共振腔内单独加载所述岩石样品下所述系统的品质因子,包括:
确定所述共振腔内单独加载所述岩石样品下所述系统的共振频率;
将所述共振腔内单独加载所述岩石样品下所述系统的共振频率,除以所述共振腔内单独加载所述岩石样品下所述系统的半功率带宽,得到对应条件下所述系统的品质因子。
6.如权利要求1所述的测量岩石衰减系数的方法,其特征在于,所述系统还包括样品移动测量装置,所述样品移动测量装置包括:
悬挂样品的传动机构,用于带动样品相对于所述共振腔上下移动;
与所述传动机构相连的步进电机,用于驱动所述传动机构;
位置控制器,用于在所述上位机的控制下向所述步进电机输出控制信号。
7.如权利要求1或6所述的测量岩石衰减系数的方法,其特征在于,所述共振腔的上端及下端敞口。
8.如权利要求1所述的测量岩石衰减系数的方法,其特征在于,所述系统还包括:
功率放大器,用于对所述锁相放大器输出的震源激励信号进行功率放大,并向所述震源输出功率放大后的震源激励信号;
预放大器,用于对所述水听器采集的共振声压场信号进行前置放大,并向所述锁相放大器输出前置放大后的共振声压场信号。
9.如权利要求7所述的测量岩石衰减系数的方法,其特征在于,所述共振腔包括圆筒形铝筒。
10.如权利要求1所述的测量岩石衰减系数的方法,其特征在于,所述敞口容器包括圆柱形玻璃桶。
11.如权利要求1所述的测量岩石衰减系数的方法,其特征在于,所述液体包括硅油。
12.如权利要求1所述的测量岩石衰减系数的方法,其特征在于,所述震源包括压电陶瓷片,其将所述震源激励信号转化为机械振动并以声波信号形式发出。
13.如权利要求6所述的测量岩石衰减系数的方法,其特征在于,所述传动机构包括滑轮组。
14.如权利要求1所述的测量岩石衰减系数的方法,其特征在于,所述震源位于所述敞口容器底部中心。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710691134.6A CN107389794B (zh) | 2017-08-14 | 2017-08-14 | 测量岩石衰减系数的方法与系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710691134.6A CN107389794B (zh) | 2017-08-14 | 2017-08-14 | 测量岩石衰减系数的方法与系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107389794A CN107389794A (zh) | 2017-11-24 |
CN107389794B true CN107389794B (zh) | 2020-01-10 |
Family
ID=60355689
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710691134.6A Active CN107389794B (zh) | 2017-08-14 | 2017-08-14 | 测量岩石衰减系数的方法与系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107389794B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108287118B (zh) * | 2017-12-26 | 2020-09-08 | 中国石油天然气股份有限公司 | 一种岩石低频弹性模量及衰减系数测量方法及系统 |
CN111637962B (zh) * | 2020-06-05 | 2021-04-20 | 无锡鸣石峻致医疗科技有限公司 | 一种剪切波衰减系数测量方法与系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003184062A (ja) * | 2001-12-18 | 2003-07-03 | Kiso Jiban Consultants Kk | 計測ボーリングを応用したサンプリング条件の最適化システム |
CN102095804A (zh) * | 2010-12-01 | 2011-06-15 | 中国石油天然气集团公司 | 一种模拟测试低频岩石模量的方法和实验装置 |
CN102778695A (zh) * | 2012-07-23 | 2012-11-14 | 中国科学院地质与地球物理研究所 | 一种获取煤岩品质因子的控制方法 |
CN104570084A (zh) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | 跨尺度地震岩石物理衰减模型及预测衰减和频散的方法 |
CN106291714A (zh) * | 2016-08-31 | 2017-01-04 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 获得地层吸收衰减品质因子的方法及其地层划分方法 |
CN106814397A (zh) * | 2016-12-21 | 2017-06-09 | 长江大学 | 一种多参数联合反演计算岩石散射衰减的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9442204B2 (en) * | 2012-08-06 | 2016-09-13 | Exxonmobil Upstream Research Company | Seismic inversion for formation properties and attenuation effects |
-
2017
- 2017-08-14 CN CN201710691134.6A patent/CN107389794B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003184062A (ja) * | 2001-12-18 | 2003-07-03 | Kiso Jiban Consultants Kk | 計測ボーリングを応用したサンプリング条件の最適化システム |
CN102095804A (zh) * | 2010-12-01 | 2011-06-15 | 中国石油天然气集团公司 | 一种模拟测试低频岩石模量的方法和实验装置 |
CN102778695A (zh) * | 2012-07-23 | 2012-11-14 | 中国科学院地质与地球物理研究所 | 一种获取煤岩品质因子的控制方法 |
CN104570084A (zh) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | 跨尺度地震岩石物理衰减模型及预测衰减和频散的方法 |
CN106291714A (zh) * | 2016-08-31 | 2017-01-04 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 获得地层吸收衰减品质因子的方法及其地层划分方法 |
CN106814397A (zh) * | 2016-12-21 | 2017-06-09 | 长江大学 | 一种多参数联合反演计算岩石散射衰减的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107389794A (zh) | 2017-11-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Brocanelli et al. | Measurement of low-strain material damping and wave velocity with bender elements in the frequency domain | |
US8166801B2 (en) | Non-invasive fluid density and viscosity measurement | |
AU2013394872B2 (en) | Method and device for the concurrent determination of fluid density and viscosity in-situ | |
US9857278B2 (en) | Apparatus for and a method of characterising mechanical properties of a sample | |
EP2267486B1 (en) | Control system for marine vibrators and seismic acquisition system using such control system | |
US9008977B2 (en) | Determining fluid density | |
NO20060885L (no) | Anordning og fremgangsmate for kompensering av en Coriolis-maler | |
JP2017525968A (ja) | 気泡特性の音響推定用の方法及びデバイス | |
JP2011525619A (ja) | 弾性的および散逸的非線形性および粘弾性を測定するための局所的な非接触音響デバイス | |
CN107389794B (zh) | 测量岩石衰减系数的方法与系统 | |
CN103217358B (zh) | 一种同时获取储层岩石低频弹性性质与密度的方法及系统 | |
CN103267804A (zh) | 一种储层岩石低频岩石物理模量测量系统 | |
Lott et al. | From local to global measurements of nonclassical nonlinear elastic effects in geomaterials | |
CN102095804A (zh) | 一种模拟测试低频岩石模量的方法和实验装置 | |
Cascante et al. | A new mathematical model for resonant-column measurements including eddy-current effects | |
RU2411500C1 (ru) | Способ измерения параметров вязкоупругих жидких сред и устройство для его реализации | |
JP2004012149A (ja) | 液体物性測定装置 | |
Mert et al. | A new method to measure viscosity and intrinsic sound velocity of liquids using impedance tube principles at sonic frequencies | |
Yin et al. | Numerical and experimental investigation of a low-frequency measurement technique: differential acoustic resonance spectroscopy | |
Camacho-Tauta et al. | Frequency domain method in bender element testing–experimental observations | |
Wisniewiski | Second generation shock tube calibration system | |
Dong et al. | Numerical simulation and analysis for low-frequency rock physics measurements | |
Henann et al. | Small-amplitude acoustics in bulk granular media | |
Gan et al. | Modulation Method of Nonlinear Acoustical Imaging | |
Malarkodi et al. | Design considerations and establishment of a low frequency hydrophone calibration setup using the principle of vibrating water column |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |