CN101334380B - 油沙中油含量的声学测定方法 - Google Patents
油沙中油含量的声学测定方法 Download PDFInfo
- Publication number
- CN101334380B CN101334380B CN2008100206761A CN200810020676A CN101334380B CN 101334380 B CN101334380 B CN 101334380B CN 2008100206761 A CN2008100206761 A CN 2008100206761A CN 200810020676 A CN200810020676 A CN 200810020676A CN 101334380 B CN101334380 B CN 101334380B
- Authority
- CN
- China
- Prior art keywords
- amplitude
- oily
- sand
- oil
- saturation degree
- 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
- 238000000034 method Methods 0.000 title claims abstract description 14
- 239000003027 oil sand Substances 0.000 title claims abstract description 11
- 239000003102 growth factor Substances 0.000 claims abstract description 19
- 238000005259 measurement Methods 0.000 claims abstract description 12
- 239000004576 sand Substances 0.000 claims description 31
- 239000000523 sample Substances 0.000 claims description 23
- 229910052751 metal Inorganic materials 0.000 claims description 8
- 239000002184 metal Substances 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 5
- 238000007789 sealing Methods 0.000 claims description 4
- 230000003321 amplification Effects 0.000 claims description 3
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 3
- 230000001902 propagating effect Effects 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 3
- 239000003921 oil Substances 0.000 description 50
- 239000011521 glass Substances 0.000 description 9
- 238000012360 testing method Methods 0.000 description 4
- 239000004411 aluminium Substances 0.000 description 3
- 229910052782 aluminium Inorganic materials 0.000 description 3
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 3
- 230000009022 nonlinear effect Effects 0.000 description 3
- 239000003208 petroleum Substances 0.000 description 3
- 238000003556 assay Methods 0.000 description 2
- 238000013142 basic testing Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 229920006395 saturated elastomer Polymers 0.000 description 2
- 238000010183 spectrum analysis Methods 0.000 description 2
- 229940099259 vaseline Drugs 0.000 description 2
- 206010070834 Sensitisation Diseases 0.000 description 1
- DTOSIQBPPRVQHS-PDBXOOCHSA-N alpha-linolenic acid Chemical compound CC\C=C/C\C=C/C\C=C/CCCCCCCC(O)=O DTOSIQBPPRVQHS-PDBXOOCHSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 239000002480 mineral oil Substances 0.000 description 1
- 235000010446 mineral oil Nutrition 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000008313 sensitization Effects 0.000 description 1
- 239000002893 slag Substances 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Abstract
油沙中油含量的声学测定方法:通过观察声波在油沙中传播时产生的非线性耗散现象,发现表征该效应的声学参量——相对增长系数G和非线性耗散指数α与油沙的油饱和度s之间存在着依赖关系。根据对相对增长系数G与油饱和度s之间的对应关系以及非线性耗散指数α与油饱和度s之间的对应关系的测定,可以通过测得声波在油沙样本中传播时的相对增长系数G和非线性耗散指数α来反推得知该样本的油饱和度,从而获得该油沙的含油量。
Description
一、技术领域
本发明涉及油沙中油含量的测定方法,提出了利用声学方法对油沙中的油含量进行测定的一种新方法。
二、背景技术
随着世界经济对石油需求的不断增加,常规石油资源已不能满足石油需求的快速增长,人们纷纷把目光转向非常规石油资源,其中之一就是油沙。油沙属非金属矿产,同石油、煤炭一样是一种矿藏资源,是提炼石油类产品的优选材料。北美地区分布最广,加拿大是目前世界上最大的地面“油沙”开采国。我国油沙资源也十分丰富,分布非常广阔,在准葛尔、柴达木、鄂尔多斯、松辽、四川、二连、塔里木、吐哈等盆地都发现了油沙的分布。我国油沙点多面广,且含油率高,勘探前景十分喜人。非常规油气资源已经成为全球勘探的热点,与之相配套的评价方法、勘探技术也正在形成之中。因此,在当今“能源就是发展,能源就是经济,能源就是国力,能源就是生命”的国际大环境中,开展油沙中的矿物油含量测定的研究对油沙资源的勘探、评价和开发都具有重要的战略指导意义。
声波在沙中的传播伴随着被称为“自我增亮”的非线性效应,即大振幅声波要比小振幅声波衰减得慢,也就是说,大振幅声波的衰减系数要比小振幅声波的衰减系数来得小。这种效应在两种极端情况即干沙和完全油饱和的沙中,都表现得很明显,然而,其程度却在数量上有所不同。这意味着某些参量同沙中的油含量有着一定的依赖关系。因此,对沙中的油含量与“自我增亮”效应关系进行详细研究,并进而确定相关敏感参量是很有意义的。一旦参量被确定,它可以用来解决反问题,例如,沙中的油含量可以通过对这个参量值的测量而得到。目前未有油沙中的油含量的声学测定方法的公开。
三、发明内容
本发明的目的是:提出一种新的利用声学方法测定油沙中油含量的方法,利用声波在油沙中传播时产生的非线性耗散现象,找到一个与油含量密切相关的参量来表征油含量,从而通过该声学参量的测定反推得到油沙中的油含量。
本发明的技术方案:通过获得表征声学效应的声参量——相对增长系数G和非线性耗散指数α与油沙的油饱和度s之间存在着的依赖关系,从而获得油沙中的油含量;具体方案是:在系统的线性测量范围内,通过测量声波在不同油饱和度的沙子中传播时接收信号幅值与发射信号幅值之间的关系,测定相对增长系数G与油饱和度s之间的关系曲线以及非线性耗散指数α与油饱和度s之间的关系曲线;从而根据实际油沙样品的G值和α值,在G~s和α~s曲线上可以查询得到该样本的油饱和度,即可获得该油沙的含油量。
试验装置如图1所示:管内注满样品,用金属塞子密封两端,一端接发射换能器,另一端接压电加速度计,示波器一通道接经过前置放大的接收信号,另一通道接发射信号。发射信号选几十至几百KHz(一般取决于发射换能器的中心频率)的正弦猝发信号。该信号经功率放大器放大后接入超声换能器。换能器探头通过耦合剂(凡士林)与下端金属塞紧密耦合,向管中发射声波。通过固定在上端金属塞上的压电加速度计,将传播到玻璃管另一端的声压信号转换成电信号。这一电信号经由前级放大器放大后送往数字示波器,并通过PCI总线IEEE488接口卡采集到计算机中,完成波形采集和数据处理。管中样品是可以控制油含量的沙子样品。
本发明的测试是在发射-接收系统本身是线性的,即声波经过中空玻璃管时输出信号的幅度正比于输入信号的条件下进行的,可保证测量观察到的非线性效应完全是由油沙样品所引起的。发射信号幅度的测量:将加速度计紧贴发射换能器,测量不同正弦信号输入幅值UT(0~300mV)下的入射声波幅度U0,找出系统的线性测量范围,例如输入幅值为0~300mV时,入射声波幅度U0与输入幅值UT成正比,则系统的线性测量范围为0~300mV。以一定输入幅度(例如8mV)时的入射声波幅度作为参考小振幅入射声波幅度U0 *。
通过测量和分析声波在不同油饱和度的沙子中传播时,接收信号幅值与发射信号幅值之间的关系,就可以寻找并确定与油含量相关的敏感参量。发射信号幅值U0可以直接在信号源显示屏上读出,接收信号幅值UL可以通过对接收波形进行频谱分析来获得。输入正弦信号幅值依次从1mV到10mV,然后是每10mV递增,至300mV止。
油饱和度s是根据玻璃管中实际含油体积V与最大含油量即完全油饱和时的含油体积V0的比值得到的,s=V/V0。在一定的油饱和度(s=70%)下,管中的油就很难再被抽出,因此试验结果中的s在0.7到1之间变化。
根据试验结果,我们发现声波通过不同油饱和度的油沙时,相对增长系数不同。在小振幅或大振幅输入幅度下,输出与输入为线性关系,而在中等振幅范围内,UL比U0增长得更快。相对增长系数G定义为接收信号增长率与发射信号增长率之比。图2中油饱和度s=82%时,输入从X1到X2增长率a=X2/X1=10,而输出从Y1到Y2增长率b=Y2/Y1=18,相对增长系数G=b/a=1.8。不同油饱和度下,相对增长系数不同。由图3相对增长系数G与油饱和度s的依赖关系,我们可以根据油沙的相对增长系数反推得到油饱和度,从而测定油沙的含油率。
进一步的数据处理发现,如果进行坐标变换,令M=UL/UL *,N=U0/U0 *,其中U0和UL分别为不同振幅声波入射幅度和接收幅度,U0 *和UL *分别为参考小振幅声波入射幅度和接收幅度,则ln[ln(M/N)]与lnN之间存在着线性关系,这与理论预计是一致的。通过对声波在油沙中传播的非线性波动方程的求解,可以推导出ln[ln(M/N)]=αlnN+C,其中α称为非线性耗散指数,C是与材料性质有关的一个物理参数。对ln[ln(M/N)]和lnN的关系曲线进行直线拟合,得到的斜率即为非线性耗散指数α。α与油饱和度s之间存在着一一对应关系(见图4),因此我们也可以根据油沙的非线性耗散指数反推得到油饱和度。
本发明的特点和有益效果是:利用声学方法对油沙中的油含量进行测定。通过观察声波在油沙中传播时产生的非线性耗散现象,发现表征该效应的声学参量——相对增长系数G和非线性耗散指数α与油沙的油饱和度s之间存在着依赖关系,从而可以通过声学检测的方法来获得油沙中的油含量。该检测方法装置简单,操作方便,所需样品量少。此发明也可以推广到其他多孔媒质,例如土壤、粘土、矿渣、岩石等的水渗透度或油饱和度检测。
四、附图说明
图1是本发明实验装置框图
图2s=100%和s=82%时输入-输出信号之间的关系
图3相对增长系数G与油饱和度s的依赖关系
图4非线性耗散指数α与油饱和度s的依赖关系
图1中圆柱形玻璃管1、金属塞2、超声发射换能器3、加速度计4、橡皮软管5。
五、具体实施方式
1.系统搭试:玻璃(样品)管内注满样品,用金属塞子密封两端,下端接发射换能器,上端接压电加速度计,示波器一通道接经过前置放大的接收信号,另一通道接发射信号。发射信号中心频率39KHz(取决于发射换能器的中心频率),脉冲持续时间30Cycles,脉冲重复频率50Hz的正弦猝发信号。该信号经功率放大器放大后接入超声换能器。换能器探头通过耦合剂(凡士林)与下端金属铝塞紧密耦合,向玻璃管中发射声波。通过固定在上端金属铝塞上的压电加速度计,将传播到玻璃管另一端的声压信号转换成电信号。这一电信号经由前级放大器放大后送往数字示波器,并通过PCI总线IEEE488接口卡采集到计算机中,完成波形采集和数据处理。实验装置框图如图1所示。玻璃管内直径30.5mm,外直径33.5mm,长度37.5cm。最初,管中装满完全油饱和的沙子,其上下端均用金属铝塞密封,沙粒的平均直径在0.2mm左右。在靠近上下端处各有一个内直径2mm的旁支管开口,用软管与其相连。软管另一端连接医用注射器,用以控制沙中的油含量。
2.基本测试:验证发射-接收系统本身是线性的,即声波经过中空玻璃管时输出信号的幅度正比于输入信号。此基本测试保证了以下试验观察到的非线性效应完全是由油沙样品所引起的。
3.发射信号幅度的测量:将加速度计紧贴发射换能器,测量不同正弦信号输入幅值UT(0~300mV)下的入射声波幅度U0,找出系统的线性测量范围,例如输入幅值为0~300mV时,入射声波幅度U0与输入幅值UT成正比,则系统的线性测量范围为0~300mV。以一定输入幅度(例如8mV)时的入射声波幅度作为参考小振幅入射声波幅度U0 *。
4.接收信号幅度的测量:玻璃管内注满完全油饱和的沙子。测量不同入射幅度U0的声波经过管长为L的沙子后的接收信号,对信号进行频谱分析,得到接收信号幅度UL。以一定输入幅度(例如8mV)时的接收信号幅度作为UL *。并进行以下数据处理:
(1)确定该油饱和度下的相对增长系数G。画出声波在沙子中传播时接收幅度UL与入射幅度U0的关系曲线。在小振幅和大振幅入射声波时,接收幅度与入射幅度成线性关系,在中间振幅时,接收幅度比入射幅度增长得快。在曲线上选取中等入射幅度时的两点:A(X1,Y1)和B(X2,Y2),计算该油饱和度下中等幅度声波的相对增长系数G=(Y2/Y1)/(X2/X1)。
(2)确定非线性耗散指数α。由测量得到的参考小振幅声波入射幅度U0 *和接收幅度UL *,以及各中等振幅声波入射幅度U0和接收幅度UL,进行变量代换,得到M=UL/UL *和N=U0/U0 *,画出ln[ln(M/N)]和lnN的关系曲线。在中等振幅条件下,两者近似为线性关系,这与理论预期相吻合。经过直线拟合,获得该直线的斜率,即为非线性耗散指数α。
5.改变油饱和度,重复上述测量,确定不同油饱和度下的相对增长系数G和非线性耗散指数α。画出G~s的关系曲线和α~s的关系曲线。油饱和度的确定:利用医用注射器抽取玻璃管中的油,首先在样品准备时记录管中注满完全油饱和的沙子时油的初始体积V0,然后精确测量从管中抽出的油的体积V1,油饱和度s=(V0-V1)/V0×100%。
6.测量实际油沙样品的G值和α值,在G~s和α~s曲线上查询得到该样本的油饱和度,即获得该油沙的含油量。
Claims (2)
1.油沙中油含量的声学测定方法:其特征是通过观察声波在油沙中传播时产生的非线性耗散现象,发现表征该效应的声学参量——相对增长系数G和非线性耗散指数α与油沙的油饱和度s之间存在着依赖关系;在系统的线性测量范围内,通过分析声波在不同油饱和度的沙子中传播时接收信号幅值与发射信号幅值之间的关系,测定相对增长系数G与油饱和度s之间的关系曲线以及非线性耗散指数α与油饱和度s之间的关系曲线,相对增长系数G定义为接收信号增长率与发射信号增长率之比,非线性耗散指数α由测量得到的参考小振幅声波入射幅度U0*和接收幅度UL*,以及不同振幅声波入射幅度U0和接收幅度UL,进行变量代换,得到M=UL/UL*和N=U0/U0*,对ln[ln(M/N)]和lnN的关系曲线进行直线拟合,获得该直线的斜率,即为非线性耗散指数α,α与油饱和度s之间存在着一一对应关系,从而根据实际油沙样品的G值和α值,在G~s和α~s曲线上可以查询得到该样本的油饱和度,即可获得该油沙的含油量。
2.由权利要求1所述的油沙中油含量的声学测定方法:其特征是具体步骤是:在样品管内注满油沙样品,用金属塞子密封两端,一端接发射换能器,另一端接压电加速度计,固定在上端金属塞上的压电加速度计将接收到的声压信号转换成电信号;这一电信号经由前置放大器放大后送往数字示波器,并通过PCI总线IEEE488接口卡采集到计算机中,完成波形采集和数据处理;示波器一通道接经过前置放大的接收信号,另一通道接发射信号;发射信号选几十至几百KHz的正弦猝发信号,该信号经功率放大器放大后接入超声换能器;换能器探头通过耦合剂与下端金属塞紧密耦合,向样品管中发射声波,管中样品是可控制油含量的沙子样品。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100206761A CN101334380B (zh) | 2008-02-21 | 2008-02-21 | 油沙中油含量的声学测定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100206761A CN101334380B (zh) | 2008-02-21 | 2008-02-21 | 油沙中油含量的声学测定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101334380A CN101334380A (zh) | 2008-12-31 |
CN101334380B true CN101334380B (zh) | 2010-08-04 |
Family
ID=40197117
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008100206761A Expired - Fee Related CN101334380B (zh) | 2008-02-21 | 2008-02-21 | 油沙中油含量的声学测定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101334380B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104181080B (zh) * | 2014-08-28 | 2016-08-24 | 湖南中科富邦科技有限责任公司 | 一种污油水含油实时监测方法及其监测装置 |
CA2935712A1 (en) | 2016-07-07 | 2018-01-07 | Sacre-Davey Innovations Inc. | System and method for real time on-stream analysis of oil sands composition |
CN108267506A (zh) * | 2016-12-30 | 2018-07-10 | 中国石油天然气股份有限公司 | 一种重质油稳定性表征的装置及方法 |
CN107179258B (zh) * | 2017-07-26 | 2020-04-03 | 江西离子型稀土工程技术研究有限公司 | 稀土金属碳含量快速检测装置 |
-
2008
- 2008-02-21 CN CN2008100206761A patent/CN101334380B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN101334380A (zh) | 2008-12-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yang et al. | Ultrasonic P-wave propagation through water-filled rock joint: an experimental investigation | |
CN101942987B (zh) | 一种简易声速测井仪检验、标定及刻度方法 | |
CN102944608A (zh) | 波纹管孔道注浆密实度超声检查的装置及方法 | |
CN102890117A (zh) | 一种长锚索锚固质量检测方法及装置 | |
CN104101647A (zh) | 一种模拟储层条件的岩石超声波速度测试系统及方法 | |
CN102788738A (zh) | 多相流体密度和浓度超声阵列检测方法 | |
CN101334380B (zh) | 油沙中油含量的声学测定方法 | |
CN104251883A (zh) | 一种非接触式检测岩石声波速度的方法 | |
CN112557499B (zh) | 一种基于超声波的节理对应力波透反射规律影响的实验方法 | |
CN105136864B (zh) | 可现场测试地表下不同深度土含水率与干密度的探测器 | |
CN102540245A (zh) | 一种地雷谐振强度测量装置及测量方法 | |
CN115793035B (zh) | 一种基于声波主动激励的地下非金属管线探测系统及方法 | |
CN204330664U (zh) | 管道内泥浆浓度电容式测量装置 | |
Cui et al. | Numerical simulation of attenuation and group velocity of guided ultrasonic wave in grouted rock bolts | |
US5412989A (en) | Acoustic tracing of buried conduits | |
CN102435298A (zh) | 地表声阻抗率测量装置及方法 | |
RU2199005C1 (ru) | Способ диагностики состояния межтрубного пространства нефтяных добывающих скважин и устройство для его осуществления | |
Weber et al. | Design of a measurement assembly to study in situ rock damage driven by freezing | |
CN112197720A (zh) | 一种基于超声导波声场调控技术的隐蔽性柱状钢制结构长度测量方法及装置 | |
CN103033258A (zh) | 远距离高指向性低频声波发射与测量装置及方法 | |
CN201554458U (zh) | 一种数字声波和变密度综合测井仪 | |
CN102042003A (zh) | 一种数字声波和变密度综合测井仪 | |
CN108195725A (zh) | 基于shpb的颗粒物质一维应力波衰减效应测定方法 | |
XU et al. | Quality factor characteristics of six metamorphic kinds of coal in China | |
CN213398327U (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: 20100804 Termination date: 20130221 |