CN102183585B - 一种岩心选样的方法 - Google Patents
一种岩心选样的方法 Download PDFInfo
- Publication number
- CN102183585B CN102183585B CN2011100563968A CN201110056396A CN102183585B CN 102183585 B CN102183585 B CN 102183585B CN 2011100563968 A CN2011100563968 A CN 2011100563968A CN 201110056396 A CN201110056396 A CN 201110056396A CN 102183585 B CN102183585 B CN 102183585B
- Authority
- CN
- China
- Prior art keywords
- rock core
- rock
- wave
- core
- cores
- 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.)
- Expired - Fee Related
Links
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明涉及一种岩心选样的方法,该方法依次包括以下步骤:对待筛选的岩心进行称重和测量;将每个岩心分别用不同频率的超声波测试,读取超声波透射过岩心所用时间、探头对接首波振幅和加岩心后的首波振幅;结合误差分析理论,得到每个岩心波阻抗的取值范围;对在不同频率下的波阻抗取值范围均出现交集的岩心,将其初步筛选出来;对初步筛选出来的岩心,利用声波衰减系数进行进一步筛选;对在不同频率下的衰减系数取值范围均出现交集的岩心,将其最终筛选出来。本发明原理可靠,结构简单,操作方便,能用较少的时间和成本优选出大量实验用性质相同的岩心,利用本发明可对岩心的质量、密度等物性和内部孔隙结构进行定量分析,在岩心分析测试中具有重大意义。
Description
技术领域
本发明涉及石油工业领域,利用不同频率超声波测试数据结合误差分析理论、岩心波阻抗和超声波衰减系数处理实现物性和孔隙结构相似岩心筛选的方法。
背景技术
在很多岩心分析测试中,由于同一块岩心在一次实验后就被破坏或者性质发生改变。因此,为了进行同步、对比实验,实验前必须准备多块性质相同(岩性、物性、结构、力学性质等都相同)的岩心。例如三轴力学实验、酸化影响及酸化效果评价实验、储层损害评价实验等。但是岩石作为一种特殊的非均质性极强的材料,其内部孔、洞、缝的发育未知,因此,如何挑选物性和结构相同或相似的岩心就成了岩心分析测试中亟待解决的基础问题。
目前行业内主要是利用地质上对岩石沉积后组系的划分,将同一组系或同一层系区块的岩心视为性质相同的岩心进行分析测试,这样做具有很明显的缺点:
(1)常规的选样方法只能保证岩心在岩性上大体一样,但岩心内部孔隙结构、机械力学性质等均无法评估,不能做定量分析,而其中任何一个因素都可能导致最终的实验结果缺乏对比性或得出错误认识。
(2)对均质性好的砂泥岩地层岩心,利用上述方法效果较好,但对于现在大量开发的低渗透致密砂岩和碳酸盐岩地层,由于其缝、洞发育,非均质极强,传统方法就显得过于粗糙。
为了解决上述岩心选样难的问题,确保岩心实验结果的科学性、正确性。本发明提出了综合利用不同频率超声波测试数据、误差分析理论、岩心波阻抗和超声波衰减系数进行处理,实现相同岩心筛选的方法。
发明内容
本发明的目的在于提供一种岩心选样的方法,该方法使用含多个不同频率超声波探头的超声波测试系统并结合误差分析理论、岩心波阻抗和超声波衰减系数处理来进行岩心选样,提高了岩心选样的效率及相似性。
为了达到上述技术目的,本发明提供以下技术方案。
所谓性质相同的岩心,主要是指质量、密度等物性和孔隙结构类似的岩心。岩石内部包含有大小不同的缝洞,而不同频率的超声波能反演岩石内部不同尺度的缝洞,裂缝对高频吸收强。表1是实际生产过程中所需检测的裂缝尺度范围,按照传统超声波测试精度,相应的探头频率能检测的裂缝尺度为λ/3(λ为波长),岩心的骨架速度在3600m/s~7000m/s范围内波动,因此本发明的测试系统中采用的探头频率为25KHz、50KHz、100KHz、250KHz、500KHz、1.25MHz、2.5MHz。
表1裂缝分类及规格
利用电子天平对实验用岩心进行称重,用游标卡尺测量实验用岩心的长度和直径,利用公式(1)(2)可以求出岩心密度:
ρ,V,M分别为岩心的密度(g·cm-3)、体积(cm3)和质量(g)D,L分别为岩心的直径(cm)和长度(cm)。
将上述长度和质量已经测试过的岩心进行多探头不同频率超声波测试,利用示波器读取超声波透射过岩心首波到达时间T(μs),并记录下首波波形以获取首波振幅用于衰减分析。利用公式(3)可得到超声波在岩心内传播的首波速度:
密度和速度这两个物理量分别表征岩心的物性和孔隙结构,在已知上述两个物理量的情况下可得到岩心的波阻抗,式(4)为波阻抗的定义式,式(5)为岩心波阻抗的计算式:
Z=ρ·V (4)
设实验中所用游标卡尺的精度为Δl、电子天平的精度为Δm、示波器读取时间的精度为Δt,由于仪器系统误差的存在致使质量、长度和时间的测量中均存在误差。上述参数(质量、直径和时间)的直接测量值加减所用测试仪器的精度将得到每个参数的一个取值范围。根据误差合成理论,质量、长度和时间这三个直接测量值的误差将传递给波阻抗这个间接量,在式(5)中带入质量的下限值和直径与时间的上限值将得到波阻抗的下限值Zmin,同理可求出波阻抗的上限值Zmax,由此得到波阻抗的一个取值范围Z=[Zmin,Zmax]。
根据上述可知每块岩心在不同频率超声波测试后都将得到一个波阻抗的取值范围,当一批岩心对不同频率下的波阻抗取值范围均出现交集时,表征这批岩心具有类似的物性和孔隙结构,以此就筛选出了性质基本相似的岩心。
但是由于单一超声波具有多解性即一个声波速度可能对应多个孔隙结构,因此我们引入声波衰减系数来进一步优选孔隙结构和强度特征类似的岩心,由式(6)可得到声波的衰减系数α:
AC,A分别为探头对接和探头间加岩心后透射波的首波振幅。
设实验用示波器读取振幅的精度为Δa,同样由于系统误差的存在致使首波振幅的读取存在误差,振幅和长度这两个直接测量值的误差将传递给衰减系数这个间接量。用振幅和长度的测量值分别加减Δa和Δl可得到首播振幅与岩心长度的一个取值范围,在(6)式中带入A0的上限值和A与L的下限值将得到衰减系数的上限值αmax,同理可求出衰减系数的下限值αmin,由此得到衰减系数的一个取值范围α=[αmin,αmax]。
不同频率指示了不同尺度范围的缝洞结构,同波阻抗分析方法类似,当一批岩心不同频率衰减系数取值范围均存在交集时,表征这批岩心具有高度相似的物性和孔隙结构。
一种岩心选样的方法,该方法依次包括以下步骤:
(1)利用电子天平对待筛选的岩心进行称重,用游标卡尺测量岩心直径和长度,得到每个岩心的质量M0、直径D0和长度L0;
(2)将每个岩心分别用探头频率为25KHz、50KHz、100KHz、250KHz、500KHz、1.25MHz、2.5MHz的超声波测试,用示波器读取不同频率超声波透射过岩心所用时间T0、探头对接首波振幅Ac0和加岩心后的首波振幅A0;
(3)设游标卡尺的精度为Δl、电子天平的精度为Δm、示波器读取时间的精度为Δt,将M0-Δm、D0+Δl、T0+Δt和M0+Δm、D0-Δl、T0-Δt分别带入下式:
可得到波阻抗的下限值Zmin和上限值Zmax,由此得到每个岩心波阻抗的取值范围Z=[Zmin,Zmax];
(4)对在不同频率下的波阻抗取值范围均出现交集的岩心,说明其具有相似的物性和孔隙结构,将其初步筛选出来;
(5)对于(4)中初步筛选出来的岩心,利用声波衰减系数进行进一步筛选,设示波器读取振幅的精度为Δa,将Ac0-Δa、A0+Δa、L0+Δl和Ac0+Δa、A0-Δa、L0-Δl分别带入下式:
将得到衰减系数的下限值αmin和上限值αmax,由此得到每个岩心的衰减系数的取值范围α=[αmin,αmax];
(6)对在不同频率下的衰减系数取值范围均出现交集的岩心,说明其具有高度相似的物性和孔隙结构,将其最终筛选出来。
与现有技术相比,本发明提供的测量方法原理可靠,结构简单,操作方便,实用性强,能用较少的时间和成本优选出大量实验用性质相同的岩心;利用本发明可对岩心的质量、密度等物性和内部孔隙结构进行定量分析,在岩心分析测试中具有重大意义。
附图说明
图1是本发明流程图。
具体实施方式
下面结合附图和实施例进一步说明本发明。
参看图1。
一种岩心选样的方法,该方法依次包括以下步骤:
(1)选择同一组系层段和区块的5块岩心,利用精度Δl=0.02mm的游标卡尺进行长度和直径测量,精度Δm=0.0001g的电子天平称重,用时间和振幅精度分别为Δt=0.4μs、Δa=0.0039的数字示波器进行读值,数据结果见表2;
(2)将(1)中处理过的岩心,利用多探头不同频率超声波测试系统进行超声波测试,将发射和接受探头分别置于岩心两端面,加适当围压让探头与岩心端面充分接触,为了较好的耦合可适当在端面涂抹凡斯林,开启脉冲发生器,开始测试;
(3)利用示波器获取不同频率超声波探头对接和透射过岩心时的时间与首波声波振幅,数据见表3、表4、表5;
(4)利用计算机数据处理系统计算表2、表3中的数据得到表6,即每个岩心不同频率下的波阻抗值范围,由表6可知岩心编号为2,3,4的三块岩心的波阻抗值对不同频率均有交集,即这三块岩心具有相似的物性和孔隙结构,将其初步筛选出来;
(5)对在(4)中筛选出来的岩心,通过衰减系数来对其做进一步筛选,由表2、表4、表5中的数据可计算得到表7,即编号为2,3,4的三块岩心在不同频率下的衰减系数范围;
(6)由表7可知岩心编号为2和3的两块岩心其衰减系数对不同频率均具有交集,即在本次岩样筛选中2号和3号岩心具有物性和孔隙结构的高度一致性,由此完成性质相同岩心的筛选工作。
表2岩心的基本物理参数记录
表3不同频率超声波透射过岩心首波时间记录(μs)
表4不同频率超声波探头对接首波振幅AC记录
表5不同频率声波透射过岩心首波振幅A记录
表6岩心对不同频率超声波的波阻抗Z范围(g·mm-2·μs-1)
表7不同频率超声波的岩心衰减系数α范围
Claims (1)
1.一种岩心选样的方法,该方法依次包括以下步骤:
(1)利用电子天平对待筛选的岩心进行称重,用游标卡尺测量岩心直径和长度,得到每个岩心的质量M0、直径D0和长度L0;
(2)将每个岩心分别用探头频率为25KHz、50KHz、100KHz、250KHz、500KHz、1.25MHz、2.5MHz的超声波测试,用示波器读取不同频率超声波透射过岩心所用时间T0、探头对接首波振幅Ac0和加岩心后的首波振幅A0;
(3)设游标卡尺的精度为Δl、电子天平的精度为Δm、示波器读取时间的精度为Δt,将M0-Δm、D0+Δl、T0+Δt和M0+Δm、D0-Δl、T0-Δt分别带入下式:
得到每个岩心波阻抗的取值范围Z=[Zmin,Zmax];
(4)对在不同频率下的波阻抗取值范围均出现交集的岩心,将其初步筛选出来;
(5)对于(4)中筛选出来的岩心,设示波器读取振幅的精度为Δa,将Ac0-Δa、A0+Δa、L0+Δl和Ac0+Δa、A0-Δa、L0-Δl分别带入下式:
得到每个岩心衰减系数的取值范围α=[αmin,αmax];
(6)对在不同频率下的衰减系数取值范围均出现交集的岩心,将其最终筛选出来。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011100563968A CN102183585B (zh) | 2011-03-09 | 2011-03-09 | 一种岩心选样的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011100563968A CN102183585B (zh) | 2011-03-09 | 2011-03-09 | 一种岩心选样的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102183585A CN102183585A (zh) | 2011-09-14 |
CN102183585B true CN102183585B (zh) | 2012-07-11 |
Family
ID=44569799
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011100563968A Expired - Fee Related CN102183585B (zh) | 2011-03-09 | 2011-03-09 | 一种岩心选样的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102183585B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103176055A (zh) * | 2013-02-07 | 2013-06-26 | 吉林大学 | 层状结构岩石垂直纹理方向介电性质测量样品的制备方法 |
CN105068120B (zh) * | 2015-07-16 | 2019-01-22 | 长江大学 | 致密砂岩裂缝的声波实验方法与识别方法 |
CN106442269B (zh) * | 2016-10-28 | 2019-01-25 | 西南石油大学 | 一种筛选室内物理模拟实验用非变量岩心的方法 |
CN106597565A (zh) * | 2016-11-05 | 2017-04-26 | 中国石油化工股份有限公司 | 一种测定岩心裂缝走向及倾向的方法 |
CN106769790B (zh) * | 2017-02-23 | 2023-10-31 | 西南石油大学 | 一种超声波作用下基于液体压力脉冲的页岩渗透率测试装置及方法 |
CN108593771A (zh) * | 2018-04-27 | 2018-09-28 | 西南石油大学 | 损伤强度计算方法和损伤强度计算装置 |
CN109596713A (zh) * | 2018-12-29 | 2019-04-09 | 西南石油大学 | 一种基于分形分维的煤岩超声波衰减系数评价方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5265461A (en) * | 1991-03-19 | 1993-11-30 | Exxon Production Research Company | Apparatuses and methods for measuring ultrasonic velocities in materials |
CN2638064Y (zh) * | 2003-07-09 | 2004-09-01 | 石油大学(华东) | 岩屑声波波速测量装置 |
CN201352212Y (zh) * | 2009-02-11 | 2009-11-25 | 中国石油天然气股份有限公司 | 岩心波速各向异性测试台 |
-
2011
- 2011-03-09 CN CN2011100563968A patent/CN102183585B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5265461A (en) * | 1991-03-19 | 1993-11-30 | Exxon Production Research Company | Apparatuses and methods for measuring ultrasonic velocities in materials |
CN2638064Y (zh) * | 2003-07-09 | 2004-09-01 | 石油大学(华东) | 岩屑声波波速测量装置 |
CN201352212Y (zh) * | 2009-02-11 | 2009-11-25 | 中国石油天然气股份有限公司 | 岩心波速各向异性测试台 |
Also Published As
Publication number | Publication date |
---|---|
CN102183585A (zh) | 2011-09-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102183585B (zh) | 一种岩心选样的方法 | |
CN108169330B (zh) | 基于非线性超声谐波法的混凝土构件轴向应力无损检测的装置和方法 | |
CN101140265B (zh) | 采用非线性共振频率的偏移测定混凝土含水量和裂纹的方法 | |
Hsiao et al. | Detecting flaws in concrete blocks using the impact-echo method | |
CN100495033C (zh) | 超声回弹综合检测山砂混凝土抗压强度方法 | |
CN104142195A (zh) | 基于超声波法的钢结构构件内部初始应力检测装置和方法 | |
CN102636568B (zh) | 一种检测混凝土内部缺陷的有限元超声成像方法 | |
CN103018337A (zh) | 一种基于面波的混凝土内部状况无损检测方法及系统 | |
CN106908177B (zh) | 一种测量各向异性材料平面应力的装置 | |
CN102954914A (zh) | 真三轴试验超声波和声发射测试系统及测试方法 | |
CN101571513A (zh) | 用于复合材料层合板质量评定的超声导波检测装置 | |
CN106802202B (zh) | 一种测量各向异性材料平面应力的方法 | |
CN107449664B (zh) | 应用超声散射波法测量单轴受压下混凝土绝对应力的方法 | |
Cui et al. | Assessing the effects of insufficient rebar and missing grout in grouted rock bolts using guided ultrasonic waves | |
CN112557499B (zh) | 一种基于超声波的节理对应力波透反射规律影响的实验方法 | |
CN102608212A (zh) | 基于声压反射系数功率谱测量薄层声阻抗和声衰减的方法 | |
CN201653999U (zh) | 一种智能超声横波纵波综合测试仪 | |
CN107315048A (zh) | 一种基于横波频谱分析的钢构件绝对应力识别方法 | |
Payan et al. | Ultrasonic methods | |
CN107271280A (zh) | 一种冲击回波‑回弹综合评定混凝土抗压强度方法 | |
Wiciak et al. | Sensor and dimensions effects in ultrasonic pulse velocity measurements in mortar specimens | |
CN103091342A (zh) | 一种对岩芯样品进行ct扫描分析处理的方法 | |
CN104897778A (zh) | 混凝土内部缺陷检测的断面网格射线平均声速法 | |
CN101819182B (zh) | 重构非均匀介质中缺陷形状的方法 | |
CN113533504A (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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120711 Termination date: 20130309 |