CN112051328B - 一种气体声弛豫吸收系数的测量方法 - Google Patents

一种气体声弛豫吸收系数的测量方法 Download PDF

Info

Publication number
CN112051328B
CN112051328B CN202010907679.8A CN202010907679A CN112051328B CN 112051328 B CN112051328 B CN 112051328B CN 202010907679 A CN202010907679 A CN 202010907679A CN 112051328 B CN112051328 B CN 112051328B
Authority
CN
China
Prior art keywords
ultrasonic
acoustic
relaxation
calculating
gas
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
Application number
CN202010907679.8A
Other languages
English (en)
Other versions
CN112051328A (zh
Inventor
张克声
李岩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guizhou Jiaotong College
Original Assignee
Guizhou Jiaotong College
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Guizhou Jiaotong College filed Critical Guizhou Jiaotong College
Priority to CN202010907679.8A priority Critical patent/CN112051328B/zh
Publication of CN112051328A publication Critical patent/CN112051328A/zh
Application granted granted Critical
Publication of CN112051328B publication Critical patent/CN112051328B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating 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/02Analysing fluids
    • G01N29/024Analysing fluids by measuring propagation velocity or propagation time of acoustic waves
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Algebra (AREA)
  • Mathematical Analysis (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Analytical Chemistry (AREA)
  • Chemical & Material Sciences (AREA)
  • Computational Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Mathematical Optimization (AREA)
  • Health & Medical Sciences (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

本发明公开了一种气体声弛豫吸收系数的测量方法,测量超声波频率分别为f1、f2和f3的超声波在待测气体中的声速值c(f1)、c(f2)和c(f3);计算弛豫强度ε;计算弛豫频率fm;计算吸收最大值μm;计算无量纲声弛豫吸收系数μ(f);计算声速频散谱c(f);计算声弛豫吸收系数α(f),本发明以三个频点上易于实时精确测量的声速测量值计算得到随频率变化的声弛豫吸收系数,其计算出的声弛豫吸收结果无需进行其他原因补偿,从而便可解决当前基于声弛豫现象的气体传感技术难以低成本地实时精确测量声弛豫吸收的技术瓶颈。

Description

一种气体声弛豫吸收系数的测量方法
技术领域
本发明属于超声波气体传感检测领域,尤其是涉及一种气体声弛豫吸收系数的测量方法。适用可激发气体(非单原子分子气体,例如:氮气、氧气、氯气甲烷、二氧化碳、六氟化硫)中的声弛豫吸收的测量。
背景技术
可激发(非单原子分子气体)气体中,声传播引起气体分子间发生非弹性碰撞,使得分子内外自由度温度跟随声压波动而变化;分子振动模式因较大的能级间隔造成振动温度的变化滞后于平动自由度温度的波动变化,需要一段弛豫时间才能完成回到热平衡态而出现弛豫过程.弛豫过程使得滞留在振动模式中的声激发能在声膨胀时进行热弛豫而发生损耗,从而造成随频率变化的声速频散和声弛豫吸收[1].半个多世纪以来,声弛豫吸收谱的测量一直是研究诸如分子振动弛豫能量转移速率等微观特性的一种重要的方法[2][3],同时它在气体定性、定量检测(例如:应用于天然气等级检测、六氟化硫泄露检测等)和气体腔体压强检测方面也具有广泛的应用[4][5][6][7][8].Ejakov等[9]利用四组不同频率的声换能器,通过多次改变发送和接收声换能器的距离和变化气体腔体压强,测量了包含氮气、氧气、氢气、二氧化碳和甲烷的多组气体的声吸收谱。随后,Petculescu等[10]设计了基于声吸收系数的气体传感器实验原型,利用一对固定距离的215KHz声换能器测量了二氧化碳、甲烷和空气混合气体的声吸收谱。但是,该传感器实验原型的实时性差(气体腔体压强需要从0.2atm依次变化至32atm),且单点声吸收测量值与理论值的最大误差可达44%。
为了避免传统测量方法需要不断改变气体腔体压强的问题,Petculescu和Lueptow(PL)[3]提出通过两频点上声速频散和声弛豫吸收系数的测量值合成气体主弛豫过程有效热容,进而重建整条声弛豫吸收谱的重建算法.但是PL算法需要获知气体平衡态压强和密度,当气体成分未知或成分比例未知时,对于气体密度的测量需求无疑会极大地增加测量装置的制作难度.针对PL算法的这一问题,张克声等[11]中提出了一个仅需测量两频点声吸收和声速值,既无需改变气体腔体压强又无需测量气体压强和密度,就可以获得整个频率范围声弛豫吸收谱的特征点重建算法.然而,该算法和PL算法一样,都需要精确测量两个频点上的声弛豫吸收系数。
低成本地实时精确测量声弛豫吸收仍是一项充满困难的工作,究其原因有二:(1)当发送和接收声换能器距离固定时,因它们具有不同的放大电路(例如:具有不同的放大系数和元器件参数),使得硬件电路上难以直接实现对声压幅度衰减的精确测量(此时声弛豫吸收系数实际测量的相对误差可高达44%左右[10]);如果为提升声弛豫吸收系数的测量精度和覆盖声弛豫吸收谱的频率范围,则需要数十次地改变收发声换能器之间的距离(通过接收信号幅度随距离变化的关系曲线求斜率而获得声弛豫吸收系数)和较大范围地改变气体腔体压强(例如:0.2-32atm)[22][23],进而导致测试系统实时性很差、结构复杂和制作成本高昂[9][10]。(2)在可激发气体中超声波换能器直接测量的是声衰减系数而非声弛豫吸收系数。声衍射损失、声散射损失和声吸收损失均可引起声衰减,且声吸收是声弛豫吸收和声经典吸收之和[1]。对于大多数可激发气体,虽然在低频(f/P≤106Hz/atm)时可以不考虑声经典吸收对声吸收的影响,但在高频(f/P>106Hz/atm)时声经典吸收则会取代声弛豫吸收占主导地位而不能被忽略[1]。因此,如果要精确获得声弛豫吸收系数,则需要从测量得到的声衰减中声经典吸收和补偿其他原因造成的衰减。然而,计算声经典吸收和其他原因造成的衰减,不仅需要预先知道气体成分构成和相关的气体分子参量(但这往往与气体检测的目的相背),且计算和补偿过程复杂[9][10]
所以,不论从测试系统的硬件电路设计要求,还是从测量结果构成的理论分析来看,都说明低成本地实时精确测量声弛豫吸收的工作具有很大的挑战性,难以满足实际气体检测中对传感器的精确度、实时性和制作成本的要求。
反观可激发气体中声速频散的测量(不同频点上的声速测量即为声速频散的测量),可在发送和接收声换能器距离固定时,通过声传播的飞渡时间(TOF)精确获得,其相对误差可小于0.01%[12];测量的实时性高、硬件模块成本低,与声弛豫吸收的测量技术相比较,其测量难度大为降低(例如:已广泛应用的超声测距仪就是通过精确测量TOF获得距离信息);另外,由于声速频散可认为只与分子弛豫过程相关[1][22],具有无需进行其他原因补偿的优点。
发明内容
针对目前现有技术中存在的上述问题,本发明提供一种气体声弛豫吸收系数的测量方法。
本发明的上述目的通过以下技术方案实现:
一种气体声弛豫吸收系数的测量方法,包括以下步骤:
步骤1、超声波频率分别为f1、f2和f3的超声波在待测气体中的声速值c(f1)、c(f2)和c(f3);
步骤2、根据下式计算弛豫强度ε:
其中,A=[c2(f3)-c2(f1)]f1 2
步骤3、根据下式计算弛豫频率fm
步骤4、根据下式计算吸收最大值μm
步骤5、根据下式计算无量纲声弛豫吸收系数μ(f):
其中,f为声频率;
步骤6、根据下式计算声速频散谱c(f):
步骤7、根据下式计算声弛豫吸收系数α(f):
α(f)=μ(f)f/c(f)。
声速值c(f1)、c(f2)和c(f3)通过以下步骤获得:
在待测气体中放置超超声波频率分别为f1、f2和f3的三组超声探头,每组超声波探头均包括超声波发射探头和超声波接收探头,通过单片机记录超声波发射探头和超声波接收探头之间的距离x,通过单片机记录超声波发射探头发送超声波的时间与超声波接收探头接收超声波的时间之间的时间差t,根据距离x和时间差t既可以获得超声波频率分别为f1、f2和f3的超声波在待测气体中的声速值c(f1)、c(f2)和c(f3)。
本发明相对于现有技术,具有以下有益效果:
本发明以三个频点上易于实时精确测量的声速测量值计算得到随频率变化的声弛豫吸收系数,其计算出的声弛豫吸收结果无需进行其他原因补偿,从而便可解决当前基于声弛豫现象的气体传感技术难以低成本地实时精确测量声弛豫吸收的技术瓶颈,进而可大幅提升传感器的测量精度、降低传感器的设计复杂度和成本,为探测气体分子微观特性提供一种更为简易和成本低廉的声学方法。
附图说明
图1为实施例2获得的声速频散谱(实线)示意图,其中圆圈:f1=40kHz声速测量值;菱形:f2=215kHz声速测量值;星型:f3=1.4MHz声速测量值;
图2为实施例2获得的声弛豫吸收谱(实线)与实测数据(点)的比对示意图,实测数据来源:J.Acoust.Soc.Am.,2006,120(4):1779-1782;
图3为实施例2获得的依赖于声频率的声弛豫吸收系数示意图;
图4为实施例1的流程示意图。
具体实施方式
下面结合实施示例对本发明作进一步的详细描述,但此处所描述的实施示例仅用于说明和解释本发明,并不用于限定本发明。
实施例1:
可在待测气体中放置测量设备,测量设备包括超声波频率分别为f1、f2和f3的三组超声探头、气压表和热电偶。每组超声波探头均包括超声波发射探头和超声波接收探头。
热电偶用于测量待测气体的温度值;
气压表用于测量待测气体的压强值;
三对超声探头分别用于测量超声波频率分别为f1、f2和f3的超声波在待测气体中的声速值,即获得声测量值,用于合成声弛豫吸收谱。
超声波频率f1、f2和f3气可根据气体种类和所测压强范围进行具体设置。例如,对于含有CH4、CO2、CL2等强弛豫可激发单一或混合气体,如果气体环境压强为1atm(1个标准大气压)时,超声波频率f1和f2可在10kHz-500kHz之间选择,超声波频率f3大于1.2M Hz;如果对于气体压强为P atm时,将上述气体压强为1atm条件下的各个超声波频率f1、f2和f3乘以P值即可。
超声波声速c可采用超声波发射探头和超声波接收探头距离一定时,测量声波传播时间的时差法进行测量,例如:每组超声波探头的超声波发射探头和超声波接收探头的距离均固定为x时,可通过一个单片机控制的计时器记录超声波在超声波发射探头和超声波接收探头之间的传播时间t,利用公式c=x/t来计算超声波频率分别为f1、f2和f3的超声波在待测气体中的声速值c(f1)、c(f2)和c(f3)。
对于气体分子单一振动弛豫过程所引起的声弛豫吸收谱,即依赖于声频率f的无量纲声弛豫吸收系数μ(f)=α(f)λ(f),其中α(f)是依赖于声频率的声弛豫吸收系数,单位是Np/m,f为声频率;λ(f)是声频率f对应的声波波长,且有λ(f)=c(f)/f,可表示为:
式(1)说明μ(f)可由弛豫频率fm和吸收最大值μm所确定。而依赖于声频率的声速频散谱c(f)为:
这里c(f3)表示超声波频率为f3的超声波在待测气体中的声速值,且满足c(f3)≈c(c是当f→∞时的声速高频极限值);ε是弛豫强度,且有:
假设待测气体为可激发气体,超声波频率分别为f1、f2和f3的超声波在待测气体中的声速值c(f1)、c(f2)和c(f3)。将声速值c(f1)、c(f2)和c(f3)代入(2)式,可得:
由(4)和(5)式可得:
其中,A=[c2(f3)-c2(f1)]f1 2,
利用(6)式,可求解得到弛豫强度ε为:
将(7)式代入(4)式,可计算得到弛豫频率fm为:
再将(7)式代入(3)式便可计算得到吸收最大值μm。再将计算获得的弛豫频率fm和吸收最大值μm代入(1)式就可获得无量纲声弛豫吸收系数μ(f);将ε、μm和c(f3)代入(2)式就可获得声速频散谱c(f)。
由于依赖于频率的声弛豫吸收系数α(f)可表示为:
α(f)=μ(f)/λ(f)=μ(f)f/c(f) (9)
由(9)式知,可通过声弛豫吸收谱μ(f)和声速频散谱c(f)进一步获得声弛豫吸收系数α(f)。
实施例2:
在本实施例中,利用实施例1的方法:当待测气体为98%CH4-2%air,在三个频点f1=40kHz,f2=215kHz,f3=1.4,测量得到的声速分别为c(f1)=447.446m/s,c(f2)=450.626m/s和c(f3)=451.831m/s;热电偶测量测量的环境温度T=300K;气压表测量的压强P=1atm。
利用(7)式计算得到的弛豫强度ε=0.0213,再将利用(8)式计算得到弛豫频率fm=1.241×105Hz,以及利用(3)式计算得到吸收最大值μm=0.0339。然后,基于计算获得的fm,μm和ε:利用(2)式便可重建声速频散谱c(f)(见图1);利用(1)式可重建声弛豫吸收谱μ(f)(见图2);利用(9)式重建依赖于频率的声弛豫吸收系数(见图3)。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。

Claims (1)

1.一种气体声弛豫吸收系数的测量方法,其特征在于,包括以下步骤:
步骤1、超声波频率分别为f1、f2和f3的超声波在待测气体中的声速值c(f1)、c(f2)和c(f3);
步骤2、根据下式计算弛豫强度ε:
其中,
步骤3、根据下式计算弛豫频率fm
步骤4、根据下式计算吸收最大值μm
步骤5、根据下式计算无量纲声弛豫吸收系数μ(f):
其中,f为声频率;
步骤6、根据下式计算声速频散谱c(f):
步骤7、根据下式计算声弛豫吸收系数α(f):
α(f)=μ(f)f/c(f),
声速值c(f1)、c(f2)和c(f3)通过以下步骤获得:
在待测气体中放置超声波频率分别为f1、f2和f3的三组超声探头,每组超声波探头均包括超声波发射探头和超声波接收探头,通过单片机记录超声波发射探头和超声波接收探头之间的距离x,通过单片机记录超声波发射探头发送超声波的时间与超声波接收探头接收超声波的时间之间的时间差t,根据距离x和时间差t即可以获得超声波频率分别为f1、f2和f3的超声波在待测气体中的声速值c(f1)、c(f2)和c(f3)。
CN202010907679.8A 2020-09-02 2020-09-02 一种气体声弛豫吸收系数的测量方法 Active CN112051328B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010907679.8A CN112051328B (zh) 2020-09-02 2020-09-02 一种气体声弛豫吸收系数的测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010907679.8A CN112051328B (zh) 2020-09-02 2020-09-02 一种气体声弛豫吸收系数的测量方法

Publications (2)

Publication Number Publication Date
CN112051328A CN112051328A (zh) 2020-12-08
CN112051328B true CN112051328B (zh) 2024-06-07

Family

ID=73606613

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010907679.8A Active CN112051328B (zh) 2020-09-02 2020-09-02 一种气体声弛豫吸收系数的测量方法

Country Status (1)

Country Link
CN (1) CN112051328B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112903811B (zh) * 2021-01-22 2023-12-01 湖北第二师范学院 一种用于检测混合气体的浓度比的声学传感方法
CN114002317B (zh) * 2021-10-22 2024-06-11 贵州交通职业技术学院 一种双原子气体的分子振动频率的声学测量方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4452082A (en) * 1981-04-30 1984-06-05 Fujitsu Limited Ultrasonic measuring method
CN106501357A (zh) * 2016-12-06 2017-03-15 华中科技大学 一种测量混合气体分子多模式弛豫过程的方法
CN107907590A (zh) * 2017-11-14 2018-04-13 中国计量大学 一种基于超声原理的气体浓度检测方法及装置
CN110161115A (zh) * 2019-05-23 2019-08-23 华中科技大学 一种基于声速频散强度变化率的气体成分探测方法和系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120118042A1 (en) * 2010-06-10 2012-05-17 Gillis Keith A Photoacoustic Spectrometer with Calculable Cell Constant for Quantitative Absorption Measurements of Pure Gases, Gaseous Mixtures, and Aerosols

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4452082A (en) * 1981-04-30 1984-06-05 Fujitsu Limited Ultrasonic measuring method
CN106501357A (zh) * 2016-12-06 2017-03-15 华中科技大学 一种测量混合气体分子多模式弛豫过程的方法
CN107907590A (zh) * 2017-11-14 2018-04-13 中国计量大学 一种基于超声原理的气体浓度检测方法及装置
CN110161115A (zh) * 2019-05-23 2019-08-23 华中科技大学 一种基于声速频散强度变化率的气体成分探测方法和系统

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
A simple measurement method of molecular relaxation in a gas by reconstructing acoustic velocity dispersion;Ming Zhu等;Measurement Science and Technology;第4页右栏-第5页 *
Algorithm for capturing primary relaxation processes in excitable gases by two-frequency acoustic measurements;Ke-Sheng Zhang等;MEASUREMENT SCIENCE AND TECHNOLOGY;第2页右栏-第7页左栏 *
Locating the inflection point of frequency-dependent velocity dispersion by acoustic relaxation to identify gas mixtures;Xiangqun Zhang等;Measurement Science and Technology;全文 *
声弛豫过程影响下的气体平衡态热容合成方法及其在气体检测中的应用;张克声;张向群;邵芳;;声学学报(第03期);全文 *
多原子分子气体中声波弛豫衰减谱的重建算法;鄢舒;王殊;;物理学报(第07期);全文 *

Also Published As

Publication number Publication date
CN112051328A (zh) 2020-12-08

Similar Documents

Publication Publication Date Title
Bond et al. Absorption of ultrasonic waves in air at high frequencies (10–20 MHz)
JP3494992B2 (ja) ガスの音速の測定
CN112051328B (zh) 一种气体声弛豫吸收系数的测量方法
CN114459649B (zh) 一种基于压电换能器阵列的无基线数据平面应力场在线监测方法、系统、设备和介质
Gedanitz et al. An apparatus for the determination of speeds of sound in fluids
CN110455919B (zh) 一种利用非线性效应评价固固界面接触特性的方法
CN114487109A (zh) 一种基于单模态多频率信号融合的无基线数据应力在线监测方法、系统、设备和介质
Jeong et al. Optimization and validation of dual element ultrasound transducers for improved pulse-echo measurements of material nonlinearity
Ishaku et al. Temperature effects on photoacoustic carbon dioxide sensor developed using mid-IR LED
US5672828A (en) Strength determination of sheet materials by utrasonic testing
CN108548631B (zh) 一种可激发气体压力容器压强测量方法
CN112903811B (zh) 一种用于检测混合气体的浓度比的声学传感方法
Revel et al. The development of a shock-tube based characterization technique for air-coupled ultrasonic probes
CN111397721A (zh) 一种基于水面边界测振技术的同振式矢量水听器绝对校准方法与系统
Li et al. Determining the responsivity of air-coupled piezoelectric transducers using a comparative method: Theory and experiments
Lu et al. A phase-comparison method for measurement of the acoustic nonlinearity parameter
Gan et al. Air-coupled ultrasonic tomographic imaging of high-temperature flames
RU2319956C2 (ru) Способ ультразвукового контроля вязкости по муни полимеров
Zhang et al. Synthesizing the frequency-dependent sound relaxational absorption by three-frequency sound speeds in excitable gases
Zorębski et al. Ultrasonic absorption measurements by means of a megahertz–range measuring set
JP2007309850A5 (zh)
Ma et al. Stress measurement of aero-engine thin-walled catheter based on ultrasonic guided wave
Zhan et al. Spontaneous and highly accurate ultrasonic temperature measurement system for air conditioner in automobiles
CN114002317B (zh) 一种双原子气体的分子振动频率的声学测量方法
Choi et al. Ultrasonic propagation in various gases at elevated pressures

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
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20240506

Address after: 551400 No.25 yunzhan Road, Qingzhen City, Guiyang City, Guizhou Province

Applicant after: GUIZHOU JIAOTONG College

Country or region after: China

Address before: School of Artificial Intelligence and Electrical Engineering, 7th Floor, Administrative Building, Guizhou Institute of Technology, No.1 Caiguan Road, Yunyan District, Guiyang City, Guizhou Province, 550003

Applicant before: Zhang Kesheng

Country or region before: China

Applicant before: Li Yan

GR01 Patent grant