CN113639934A - 气体泄漏位置三维定位方法、定位系统 - Google Patents

气体泄漏位置三维定位方法、定位系统 Download PDF

Info

Publication number
CN113639934A
CN113639934A CN202111011720.4A CN202111011720A CN113639934A CN 113639934 A CN113639934 A CN 113639934A CN 202111011720 A CN202111011720 A CN 202111011720A CN 113639934 A CN113639934 A CN 113639934A
Authority
CN
China
Prior art keywords
sound field
information
sensor
leakage
dimensional
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
Application number
CN202111011720.4A
Other languages
English (en)
Other versions
CN113639934B (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.)
Zhengzhou University
Original Assignee
Zhengzhou University
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 Zhengzhou University filed Critical Zhengzhou University
Priority to CN202111011720.4A priority Critical patent/CN113639934B/zh
Publication of CN113639934A publication Critical patent/CN113639934A/zh
Application granted granted Critical
Publication of CN113639934B publication Critical patent/CN113639934B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M3/00—Investigating fluid-tightness of structures
    • G01M3/02—Investigating fluid-tightness of structures by using fluid or vacuum
    • G01M3/04—Investigating fluid-tightness of structures by using fluid or vacuum by detecting the presence of fluid at the leakage point
    • G01M3/24—Investigating fluid-tightness of structures by using fluid or vacuum by detecting the presence of fluid at the leakage point using infrasonic, sonic, or ultrasonic vibrations
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17D—PIPE-LINE SYSTEMS; PIPE-LINES
    • F17D5/00—Protection or supervision of installations
    • F17D5/005—Protection or supervision of installations of gas pipelines, e.g. alarm
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17D—PIPE-LINE SYSTEMS; PIPE-LINES
    • F17D5/00—Protection or supervision of installations
    • F17D5/02—Preventing, monitoring, or locating loss
    • F17D5/06—Preventing, monitoring, or locating loss using electric or acoustic means
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M3/00—Investigating fluid-tightness of structures
    • G01M3/02—Investigating fluid-tightness of structures by using fluid or vacuum
    • G01M3/04—Investigating fluid-tightness of structures by using fluid or vacuum by detecting the presence of fluid at the leakage point
    • G01M3/24—Investigating fluid-tightness of structures by using fluid or vacuum by detecting the presence of fluid at the leakage point using infrasonic, sonic, or ultrasonic vibrations
    • G01M3/243—Investigating fluid-tightness of structures by using fluid or vacuum by detecting the presence of fluid at the leakage point using infrasonic, sonic, or ultrasonic vibrations for pipes

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Examining Or Testing Airtightness (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

气体泄漏位置三维定位方法，包括：建立基于基尔霍夫衍射远场声全息理论的声场模型；在测量平面获取泄漏源声场信息；处理获取的声场信息，得到泄漏源声场的声全息数据；利用声场模型和声全息数据，确定泄漏源的三维位置信息。气体泄漏位置三维定位方法能够获得更准确的声场信息，进而能够引入深度信息实现泄漏源的三维检测效果展示；通过构建基于基尔霍夫衍射理论的远场声全息模型，能够实现较大距离泄漏源的声场定位；利用虚拟相控阵列技术，在气体泄漏检测与定位中降低了系统成本和复杂度，提高了泄漏源检测精度和定位精度，在气体泄漏位置的准确定位方面具有巨大实用价值和工业应用潜力。

Description

气体泄漏位置三维定位方法、定位系统
技术领域
本申请属于气体检测技术领域，具体涉及气体泄漏位置三维定位方法、定位系统。
背景技术
气体的运输和存储在人们的工业应用和日常生活中随处可见。但是，在气体管道或存储设备的使用过程中由于自身或其他人为因素的破坏等原因，在压缩气体的运输或存储过程中经常会产生泄漏。在许多工业应用场合，这些压缩气体通常具有易燃易爆或腐蚀毒害性等特征，一旦发生泄漏事故，往往会造成严重环境污染、巨大资源浪费和重大安全事故。能够快速准确的定位气体泄漏位置对于防止发生气体运输或储存事故至关重要。
目前，针对气体泄漏源的一维检测方法可以做到对管道及容器中泄漏孔的检测和定位。例如专利CN109813501A公开了一种测量气体管道泄漏位置的方法、装置及系统，利用较少数量的传感器保证数据采集的统一性，提高了测量精度，提高了对气体管道泄漏位置定位的有效性；然而在通常情况下，实际发生气体泄漏的环境非常复杂，由于受外界噪声干扰的影响以及泄漏声学信号自身反射波和应力波在管道或储存工具中传播的复杂机理，气体泄漏源的一维检测方式在识别和定位管道泄漏孔时其精度受到了很大的限制。
针对气体泄漏源的二维检测方法目前仍处于采用线阵或其他小型阵列的初步研究阶段，成像分辨率和定位精度低。另一方面，由于这类方法往往无法对气体泄漏源深度信息(即距离信息)进行有效估计，目前仍缺少适用于气体泄漏源的三维定位方法，无法实现泄漏源的三维定位。
发明内容
有鉴于此，一方面，一些实施例公开的技术方案是气体泄漏位置三维定位方法，该方法包括：
建立基于基尔霍夫衍射远场声全息理论的声场模型；
在测量平面获取泄漏源声场信息；
处理获取的声场信息，得到泄漏源声场的声全息数据；
利用声场模型和声全息数据，确定泄漏源的三维位置信息。
进一步，一些实施例公开的气体泄漏位置三维定位方法，在测量平面获取泄漏源声场信息具体包括：
利用超声波传感器依照预设阵列规则在测量平面形成阵元间距相等的虚拟相控传感器阵列，超声波传感器包括至少一个参考传感器和至少一个扫描传感器；
利用扫描传感器逐一获取虚拟相控传感器阵列每一个阵元位置的声场信息，利用参考传感器同步获取其所在位置的声场信息。
一些实施例公开的气体泄漏位置三维定位方法，处理获取的声场信息包括：
利用互功率谱算法处理扫描传感器和参考传感器获得的声场信息，消除扫描时间间隔，达到采集时间一致性，得到泄漏源测量平面的声全息数据。
一些实施例公开的气体泄漏位置三维定位方法，声场模型的表达式为：
Figure BDA0003238664430000021
其中，U(∈,η)表示泄漏源平面的声压分布，被定义为空间波函数h(x,y)，M、N为虚拟相控传感器阵列(M*N)的行数和列数；h*(x,y)为空间波函数h(x,y)的共轭函数，包含虚拟相控传感器阵列第m行、第n列位置的声全息信息；rmn是虚拟相控传感器阵列第m行、第n列位置到泄漏源之间的距离；ΔxΔy是虚拟相控传感器阵列第m行、第n列传感器的覆盖面积；C为全息常数；z0是泄漏源到测量平面的垂直距离；
Figure BDA0003238664430000022
为复变函数；k为波数；m≤M,n≤N。
一些实施例公开的气体泄漏位置三维定位方法，处理获取的声场信息具体包括：
对参考传感器的信息xr(t)和第i个扫描传感器的信息xi(t)进行傅里叶变换，得到如下表达式：
Figure BDA0003238664430000031
其中，
Figure BDA0003238664430000032
和
Figure BDA0003238664430000033
分别表示xr(t)和xi(t)的傅里叶变换，e-j2πft为复变函数；
进一步得到如下表达式：
Figure BDA0003238664430000034
xi(t)和xr(t)的互功率谱可以表达为：
Figure BDA0003238664430000035
其中，
Figure BDA0003238664430000036
是xr(t)的共轭谱，
Figure BDA0003238664430000037
为复变函数；
第i个扫描传感器采集的信号与参考传感器同时采集的参考信号之间的相位差
Figure BDA0003238664430000038
表达为：
Figure BDA0003238664430000039
其中，Im[R(f)]表示R(f)的虚部，Re[R(f)]表示R(f)的实部；
虚拟相控传感器阵列中第一个扫描传感器采集到的信号为X1(t)，则，第i个传感器采集到的信号Xi(t)表达为：
Figure BDA00032386644300000310
其中，
Figure BDA00032386644300000311
是X1(t)与同时接收的参考信号之间的相位差，
Figure BDA00032386644300000312
是Xi(t)与同时接收的参考信号之间的相位差，
Figure BDA00032386644300000313
是复变函数。
一些实施例公开的气体泄漏位置三维定位方法，确定泄漏源的三维位置信息包括：
利用声全息数据和声场模型进行泄漏源声场重建计算，将计算结果显示为泄漏源图像信息；
以设定重建距离间隔在不同距离上进行多次声场重建计算，显示多个泄漏源图像信息；
确定清晰度最高的泄漏源图像信息对应的泄漏源位置信息为泄漏源三维位置信息。
一些实施例公开的气体泄漏位置三维定位方法，泄漏源图像信息包括二维图像信息或三维图像信息。
一些实施例公开的气体泄漏位置三维定位方法，预设阵列规则包括阵元间距、阵元数量、阵列形状。
一些实施例公开的气体泄漏位置三维定位方法，确定泄漏源的三维位置信息包括全息显示泄漏源位置。
另一方面，一些实施例公开了气体泄漏位置三维定位系统，用于实现一些实施例公开的气体泄漏位置三维定位方法，该三维定位系统包括：
信息采集组件，包括至少一个扫描传感器和至少一个参考传感器，设置用于获取泄漏源在测量平面的声场信息；
线性移动平台，用于控制扫描传感器和参考传感器形成虚拟相控阵列；
数据处理组件，其中嵌入有声场模型，设置用于传输并处理声场信息，获取声全息数据，建立声场模型，利用建立的声场模型和处理得到的声全息数据，计算确定泄漏源的三维位置信息。
本申请实施例公开的气体泄漏位置三维定位方法，能够获得更准确的声场信息，进而能够引入深度信息实现泄漏源的三维检测效果展示；通过构建基于基尔霍夫衍射理论的远场声全息模型，能够实现较大距离泄漏源的声场定位；利用虚拟相控阵列技术，在气体泄漏检测与定位中降低了系统成本和复杂度，提高了泄漏源检测精度和定位精度，在气体泄漏位置的准确定位方面具有巨大实用价值和工业应用潜力。
附图说明
图1衍射声全息方法原理示意图
图2声场信息采集原理示意图
图3气体泄漏三维定位系统组成示意图
图4气体泄漏三维定位结果二维平面投影示意图
图5气体泄漏三维定位结果二维平面投影和三维投影示意图
图6气体泄漏三维定位检测效果图
附图标记
1 参考传感器 2 扫描传感器
3 信数据处理组件 4 线性移动平台
12 虚拟相控阵列 100 泄漏源
101 虚拟阵元 x1 横向间距
y1 纵向间距 102 空气压缩机
103 喷嘴 104 喷嘴声源全息图像
105 喷嘴声源全息图像放大图
具体实施方式
在这里专用的词“实施例”，作为“示例性”所说明的任何实施例不必解释为优于或好于其它实施例。本申请实施例中性能指标测试，除非特别说明，采用本领域常规试验方法。应理解，本申请中所述的术语仅仅是为描述特别的实施方式，并非用于限制本申请公开的内容。
除非另有说明，否则本文使用的技术和科学术语具有本申请所属技术领域的普通技术人员通常理解的相同含义；作为本申请中其它未特别注明的试验方法和技术手段均指本领域内普通技术人员通常采用的实验方法和技术手段。
本文所用的术语“基本”和“大约”用于描述小的波动。例如，它们可以是指小于或等于±5％，如小于或等于±2％，如小于或等于±1％，如小于或等于±0.5％，如小于或等于±0.2％，如小于或等于±0.1％，如小于或等于±0.05％。在本文中以范围格式表示或呈现的数值数据，仅为方便和简要起见使用，因此应灵活解释为不仅包括作为该范围的界限明确列举的数值，还包括该范围内包含的所有独立的数值或子范围。例如，“1～5％”的数值范围应被解释为不仅包括1％至5％的明确列举的值，还包括在所示范围内的独立值和子范围。因此，在这一数值范围中包括独立值，如2％、3.5％和4％，和子范围，如1％～3％、2％～4％和3％～5％等。这一原理同样适用于仅列举一个数值的范围。此外，无论该范围的宽度或所述特征如何，这样的解释都适用。本文中述及的传感器，一般指能够采集声音信号的传感器，如超声波传感器等。
在本文中，包括权利要求书中，连接词，如“包含”、“包括”、“带有”、“具有”、“含有”、“涉及”、“容纳”等被理解为是开放性的，即是指“包括但不限于”。连接词“由……构成”和“由……组成”是封闭连接词。
为了更好的说明本申请内容，在下文的具体实施例中给出了众多的具体细节。本领域技术人员应当理解，没有某些具体细节，本申请同样可以实施。在实施例中，对于本领域技术人员熟知的一些方法、手段、仪器、设备等未作详细描述，以便凸显本申请的主旨。
在不冲突的前提下，本申请实施例公开的技术特征可以任意组合，得到的技术方案属于本申请实施例公开的内容。
在一些实施方式中，气体泄漏位置三维定位方法包括：建立基于基尔霍夫衍射远场声全息理论的声场模型，在测量平面获取泄漏源声场信息，处理获取的声场信息，得到泄漏源声场的声全息数据，利用声场模型和声全息数据，确定泄漏源的三维位置信息。建立基于基尔霍夫远场全息理论的声场模型与获得泄漏源声场信息为独立进行的两个阶段，一般地，三维定位方法中不限定其先后顺序。除非在具体实施过程中表明了实际进行的次序，或者根据上下文的记载，确定了严格的先后顺序。
作为可选实施例，在测量平面获取泄漏源声场信息具体包括：利用超声波传感器依照预设阵列规则在测量平面形成阵元间距相等的虚拟相控传感器阵列，超声波传感器包括至少一个参考传感器和至少一个扫描传感器；利用扫描传感器逐一获取虚拟相控传感器阵列每一个阵元位置的声场信息，同时利用参考传感器同步获取其所在位置的声场信息。
作为可选实施例，处理获取的声场信息包括：利用互功率谱算法处理扫描传感器和参考传感器获得的声场信息，消除不同扫描传感器的扫描时间间隔得到采集时间一致性，得到泄漏源测量平面的声全息数据。
作为可选实施例，声场模型表达式为：
Figure BDA0003238664430000061
其中，U(∈,η)表示泄漏源平面的声压分布，被定义为空间波函数h(x,y)，M、N为虚拟相控传感器阵列(M*N)的行数和列数；h*(x,y)为全息函数，包含虚拟相控传感器阵列第m行、第n列位置的全息信息；rmn是虚拟相控传感器阵列第m行、第n列位置到泄漏源之间的距离；ΔxΔy是虚拟相控传感器阵列第m行、第n列传感器的覆盖面积；C为全息常数；z0是泄漏源到测量平面的垂直距离；
Figure BDA0003238664430000071
为复变函数；k为波数；m≤M,n≤N。
声场模型的建立
在气体泄漏源的声全息分析中，目标源辐射的声场中存在两个重要的平行面：声源面和测量面，它们是平面声全息方法中重要的两个平行平面。利用建立的远场声全息模型，通过投影算子能够关联两个空间分离平面的声场，即通过了解一个平面中的场，可以用数值估计出与其平行的另一个平面中的场。所以，通过在测量面测得的测量场可以计算得到声源面的辐射场，从而从声源面的辐射场识别泄漏位置。
以下示例性说明利用压力场传播的线性模型在二维平面测量的函数中构建声场的全息表达。
声压波P在线性不可压缩流体中的传播可以用亥姆霍茨方程表达：
Figure BDA0003238664430000072
其中，U(P)是观测点P(x,y,z)的复振幅，k是波数，并由下式表示：
Figure BDA0003238664430000073
其中，ω是信号角频率，λ是泄漏信号的波长，c是声波的传播速度。可以看出亥姆霍兹方程与时间因素无关，它描述了不随时间变化的稳定声场。
经典的无穷基尔霍夫衍射公式如下式：
Figure BDA0003238664430000074
其中，U是空间中P(x,y)点的声场分布，可以被定义为空间波函数h(x,y)；n是测量平面的法向量；r是半径向量r的模，eikr表示复变函数；
当h(x,y)给定时，基尔霍夫衍射积公式可以计算出被动空间中波传播方向任意点P处的场。利用这个特性，可以从已知平面的波场分布推导出下一个平面上的波场分布。例如，可以通过以下实施例实现。
如图1所示，θ是半径向量r和测量平面法向量n之间的夹角，∑为空间中测量平面位置处的无限平面，H是该无限平面上的有限测量部分，即测量平面；R为声源点S(∈,η)所在的平面，可以称为声源平面，∈,η分别表示x和y轴上的坐标。
对空间波函数h(x,y)进行共轭处理后可以得到h*(x,y)，从而可以得到包含测量信号共轭波信息的全息函数H(x,y)，共轭波继续传播，会聚在声源S(∈,η)的实像S1上(S1可以看作是虚拟声源)。共轭波的传播如图1(b)所示。
法向量n、半径向量r和z轴的方向向量z之间的关系由下式给出：
Figure BDA0003238664430000081
Figure BDA0003238664430000082
全息函数H(x,y)可以表达为：
H(x,y)＝U＝h*(x,y)ejkz ……(6)
其中h*(x,y)是h(x,y)的共轭表示，k是波数，ejkz表示复变函数。
将式(4)、(5)和(6)代入式(3)，可得到声源平面内的声压分布：
Figure BDA0003238664430000083
积分全息平面xOy是无穷的，但在实际测量中，无法测量平面内连续的声压点，因此，声压分布方程应当离散化处理。例如，在测量平面内设置虚拟相控传感器阵列，该阵列中的阵列单元传感器代表离散的声音信息测量点，离散化后的声压分布公式为：
Figure BDA0003238664430000084
其中，U(∈,η)表示泄漏源平面的声压分布，被定义为空间波函数h(x,y)，M、N为虚拟相控传感器阵列(M*N)的行数和列数；h*(x,y)为全息函数，包含虚拟相控传感器阵列第m行、第n列位置的全息信息；rmn是虚拟相控传感器阵列第m行、第n列位置到泄漏源之间的距离；ΔxΔy是虚拟相控传感器阵列第m行、第n列传感器的覆盖面积；C为全息常数；z0是泄漏源到测量平面的垂直距离；
Figure BDA0003238664430000091
为复变函数；k为波数；m≤M,n≤N。
对每个泄漏源点重复相同的过程，可以获得整个泄漏源平面的声压分布。
声场信息的采集与处理
构建虚拟相控阵列
如图2所示，在声场信息测量平面内设置虚拟相控阵列12，该相控阵列中包括M×N个阵元作为声场信息采集点，形成M*N阶虚拟相控阵列，其中所有阵元之间等间隔设置，M×N个采集点形成方形的虚拟相控阵列12；泄漏源100设置在三维直角坐标系空间中，在泄漏源100的右侧适当位置形成虚拟相控阵列12，在虚拟相控阵列12中设置一个位置固定的参考传感器1，同时设置一个位置可移动的扫描传感器2，参考传感器1和扫描传感器2与数据处理组件3连接；
利用扫描传感器2逐一地采集虚拟相控阵列12中各个采集点的声场信息，同时利用参考传感器1同步测量参考传感器设置位置的声场信息；其中，参考传感器1采集的信息表示为xr(t)，扫描传感器2在第i个采集点采集的信息表示为xi(t),分别表达为：
xr(t)＝α1s(t-τi-tsr)+nr(t)(i＝1,2,…n) ……(9)
xi(t)＝α2s(t-τi-tsi)+ni(t)(i＝1,2,…n) ……(10)
其中，s(t)是采集到的声场信号，由泄漏源产生；α1和α2是衰减因子；nr(t)和ni(t)分别是参考传感器和扫描传感器的噪声；τi是扫描传感器在第i采集点与第i-1个采集点之间的扫描时间间隔。
可以得到扫描传感器2在第i个采集点采集信息的时间与参考传感器1采集信号的时间差Δti为：
Δti＝tsr-tsi(i＝1,2,…,n) ……(11)
由上式可知，Δti的表达式中不包含扫描时间间隔τi，说明τi与扫描时间间隔无关。因此，采用互功率谱法可以消除虚拟相控阵列中扫描传感器逐一扫描的扫描时间间隔，得到与传统传感器阵列相同的型号模型。
声场信息处理
作为可选实施例，利用互功率谱法处理获取的声场信息具体包括：
对参考传感器的信息xr(t)和扫描传感器在第i个信息采集点的信息xi(t)进行傅里叶变换，得到如下表达式：
Figure BDA0003238664430000101
其中，
Figure BDA0003238664430000102
和
Figure BDA0003238664430000103
分别表示xr(t)和xi(t)的傅里叶变换，e-j2πft为复变函数；
进一步得到如下表达式：
Figure BDA0003238664430000104
xi(t)和xr(t)的互功率谱可以表达为：
Figure BDA0003238664430000105
其中，
Figure BDA0003238664430000106
是xr(t)的共轭谱，
Figure BDA0003238664430000107
为复变函数；
扫描传感器在第i个信息采集点采集的信号与参考传感器同时采集的参考信号之间的相位差
Figure BDA0003238664430000108
表达为：
Figure BDA0003238664430000109
其中，Im[R(f)]表示R(f)的虚部，Re[R(f)]表示R(f)的实部；
虚拟相控传感器阵列中扫描传感器在第一个信息采集点采集到的信号为X1(t)，则，早第i个信息采集点采集到的信号Xi(t)表达为：
Figure BDA00032386644300001010
其中，
Figure BDA00032386644300001011
是X1(t)与同时接收的参考信号之间的相位差，
Figure BDA00032386644300001012
是Xi(t)与同时接收的参考信号之间的相位差，
Figure BDA00032386644300001013
是复变函数。
通过对每个阵元重复相同的过程，可以获得类似于传统传感器阵列的全息数据。利用全息数据，结合声全息方法，利用声场模型计算整个泄漏源平面的声压分布。
一些实施例公开的气体泄漏位置三维定位方法，确定泄漏源的三维位置信息包括：
利用声全息数据和声场模型进行泄漏源重建计算，将计算结果显示为泄漏源图像信息；
以设定重建距离间隔多次进行重建计算，显示多个泄漏源图像信息；通常以第一次重建计算的距离为基础，以等间隔的距离，距离增加和距离减小的两个方向进行泄漏源重建，显示不同重建距离处的多个泄漏源图像信息；
确定清晰度最高的泄漏源图像信息对应的泄漏源位置信息为泄漏源三维位置信息。通常图像的清晰度最高，意味着成像过程的聚焦效果最好。
作为可选实施例，泄漏源图像信息包括二维图像信息或三维图像信息。
作为可选实施例，预设阵列规则包括阵元间距、阵元数量、阵列形状。一般地可以根据泄漏源的具体情况，如测量距离、测量空间等合理设置阵元数量、阵元间距，设置合理的阵列形状，以满足泄漏源的测量精度和测量效率等需求。
一些实施例公开了气体泄漏位置三维定位系统，用于实现气体泄漏位置三维定位方法，该三维定位系统包括：
信息采集组件，包括至少一个扫描传感器和至少一个参考传感器，设置用于获取泄漏源在测量平面的声场信息；一般地，设置一个参考传感器和一个扫描传感器即可实现虚拟相控阵列的设置和声场信息采集；也可以设置多个扫描传感器和参考传感器实现相控阵列的设置和声场信息采集，一般地，传感器的数量要远小于虚拟相控阵列的阵元数量，可以有效地发挥虚拟相控阵列的优势；
线性移动平台，用于控制所述扫描传感器和所述参考传感器形成虚拟相控阵列；一般地线性移动平台可以精确控制扫描传感器和参考传感器的位置和移动，依照设定的规则控制传感器在虚拟相控阵列中移动，采集每一个相控阵元位置的声场信息；
数据处理组件，设置用于传输并处理所述声场信息，获取声全息数据，建立声场模型，利用建立的声场模型和处理得到的声全息数据，计算确定泄漏源的三维位置信息。一般地，数据处理组件中嵌入或植入有已经建立的声场模型，例如声场模型可以植入上位机中；通常，数据处理组件包括数据采集部件，包括多通道数据采集卡，采集程序面板等，可以将传感器采集的信息传输至上位机；上位机用于操控传感器和线性移动平台形成虚拟相控阵列并进行声场信息采集，同时控制数据采集部件；进一步，上位机还根据建立的远场声全息模型处理数据，显示三维定位结果。
气体泄漏三维定位系统构建
如图3所示，设置一参考传感器1和一扫描传感器2，利用线性移动平台4控制扫描传感器2和参考传感器1，将参考传感器1和扫描传感器2设置在三维直角坐标系XYZ空间中，在监测平面XOY内设置虚拟相控阵列12，该虚拟相控阵列12包括多个等间隔设置的虚拟阵元101，虚拟阵元101横向间隔为x1，纵向间隔为y1；设置数据处理组件3，将其与扫描传感器2和参考传感器1连接，并与线性移动平台4连接；数据处理组件3通常包括模拟量输入通道、ADC芯片、PFGA芯片、PCI控制器、存储器等，还可以包括精密基准源、校准电路等；
线性移动平台4控制扫描传感器2在监测平面内移动，能够逐一在每个虚拟阵元101位置采集声场信息；参考传感器1设置在坐标系原点O位置，与扫描传感器2同步采集该位置的声场信息；
气体泄漏三维定位系统采集数据过程中，利用参考传感器1采集坐标系原点O位置的声场信息，扫描传感器2依次逐个采集虚拟阵元101位置的声场信息，扫描传感器2采集每一个虚拟阵元101位置的声场信息时，参考传感器1同步采集坐标系原点位置的声场信息；参考传感器1和扫描传感器2采集的声场信息输入到数据处理组件3，通过模拟量输入通道一和模拟量输入通道二，传输到数据处理组件3的信息存储部件存储并经由信息处理部件进行互功率谱计算，消除扫描时间间隔，获取全息数据；根据建立的声全息模型和获得的声全息数据，计算出泄漏源100声场的声压分布，获取泄漏源深度信息，实现三维定位。
具体地，以空气压缩机作为目标气源，将目标气源从喷嘴中喷出形成泄漏源，设置在三维直角坐标系XYZ空间，设定横向间距x1为3mm，纵向间距y1为3mm，设置形成包含60×60虚拟阵元的虚拟相控阵列，设置相邻传感器之间的间距小于声波型号的半波长，确保满足空间采样定理限制，使用全数字双通道记录器以采样率为1MHz的频率记录泄漏声场信号，测量平面位于距声源zh＝0.6m处，且测量平面为边长为18cm的正方形平面。测量平面完成泄漏声场数据测量后，通过后处理对泄漏源声场进行三维检测和成像。
利用距离泄漏源0.6m的测量平面处获得的声场信息和声场模型，构建多个距离处的声场，并进行三维成像，每次声场重建的距离间隔0.15m。结果如图4所示，图中，距离z表示与测量平面之间的距离，图中泄漏源图像的边界越清晰，说明泄露声场越强(实际操作中，泄漏源图像显示为红色，颜色越红表示泄漏声场越强)，据此可直观判断泄漏声场最强的的位置为泄漏源。
图4中每隔0.15m进行一次声场重建，可以看出，随着重建距离z的增加，重建图像逐渐收敛；当重建距离z超过实际距离0.6m，重建图像逐渐发散，分辨率逐渐降低，直至失去意义。可见，只有当重构距离z等于实际距离时，才能很好地重构泄漏源。相反，在其他距离上定位形成的图像显示出强烈的伪影，泄漏源未得到很好的重建。通过在不同的深度上反复进行泄漏声场重构，以此达到泄漏源深度信息的有效估计，从而实现对泄漏源的三维定位。可以确定，泄漏源的位置为z＝0.6m。
图5中，显示了泄漏源在0.3m，0.6m和1.2m距离处的定位结果，包括二维平面投影成像及对应的三维成像。
图6为气体泄漏三维定位检测效果图，直接展示了实际测试现场的定位效果。图6中，空气压缩机102设置与喷嘴103连接，形成气体泄漏源，测量平面完成泄漏声场数据测量后，对获得的数据进行处理，并进行三维检测和成像，在喷嘴气体泄漏源位置直接显示喷嘴声源全息图像104，喷嘴声源全息图像放大图105是将喷嘴声源全息图像104放大显示的效果。
定位系统可靠性试验
为了验证三维定位系统的稳定性和可靠性，对定位系统记性了多组重复实现，泄漏源位于0.6m处，并对测试结果进行了误差分析，结果如表1。
表1 三维定位系统结果及分析
Figure BDA0003238664430000141
在x和y方向定位的相对误差(％)定义为泄漏源估计位置与实际位置的差值与重建平面尺寸的比值。从表1中可看出该方法在x和y方向的定位误差约为5％。
为了分析该方法在z方向的定位精度，将测量平面前方声场定义为一个0.6m×0.6m×1.2m的长方体空间，并将其离散化为60×60×24网格(z方向上划分为24个网格)。从表1中可以看出z方向上定位误差大约2-3个网格。误差统计表明该定位系统的误差表现可以满足大多数实际工程需要，并且定位效果稳定、可靠。
一般地，随着气体泄漏源的距离增加，泄漏源最终成像的分辨率略有下降；通常这种情况的发生与检测距离增加造成的噪声、声反射等负面影响因素有关。
本申请实施例公开的气体泄漏位置三维定位方法，能够获得更准确的声场信息，进而能够引入深度信息实现泄漏源的三维检测效果展示；通过构建基于基尔霍夫衍射理论的声全息模型，能够实现较大距离泄漏源的声场定位；利用虚拟相控阵列技术，在气体泄漏检测与定位中降低了系统成本和复杂度，提高了泄漏源检测精度和定位精度，在气体泄漏位置的准确定位方面具有巨大实用价值和工业应用潜力。
本申请公开的技术方案和实施例中公开的技术细节，仅是示例性说明本申请的发明构思，并不构成对本申请技术方案的限定，凡是对本申请公开的技术细节所做的常规改变、替换或组合等，都与本申请具有相同的发明构思，都在本申请权利要求的保护范围之内。

Claims (10)

1.气体泄漏位置三维定位方法，其特征在于，包括：
建立基于基尔霍夫衍射远场声全息理论的声场模型；
在测量平面获取泄漏源声场信息；
处理获取的声场信息，得到泄漏源声场的声全息数据；
利用所述声场模型和所述声全息数据，确定泄漏源的三维位置信息。
2.根据权利要求1所述的所述气体泄漏位置三维定位方法，其特征在于，所述在测量平面获取泄漏源声场信息具体包括：
利用超声波传感器依照预设阵列规则在测量平面形成阵元间距相等的虚拟相控传感器阵列，所述超声波传感器包括至少一个参考传感器和至少一个扫描传感器；
利用所述扫描传感器逐一获取所述虚拟相控传感器阵列每一个阵元位置的声场信息，利用所述参考传感器同步获取其所在位置的声场信息。
3.根据权利要求2所述的气体泄漏位置三维定位方法，其特征在于，所述处理获取的声场信息包括：
利用互功率谱算法处理所述扫描传感器和所述参考传感器获得的声场信息，消除扫描时间间隔达到采集时间一致性，得到泄漏源测量平面的声全息数据。
4.根据权利要求3所述的气体泄漏位置三维定位方法，其特征在于，所述声场模型的表达式为：
Figure FDA0003238664420000011
其中，U(∈,η)表示泄漏源平面的声压分布，被定义为空间波函数h(x,y)，M、N为虚拟相控传感器阵列(M*N)的行数和列数；h*(x,y)为空间波函数h(x,y)的共轭函数，包含虚拟相控传感器阵列第m行、第n列位置的全息信息；rmn是虚拟相控传感器阵列第m行、第n列位置到泄漏源之间的距离；ΔxΔy是虚拟相控传感器阵列第m行、第n列传感器的覆盖面积；C为全息常数；z0是泄漏源到测量平面的垂直距离；
Figure FDA0003238664420000012
为复变函数；k为波数；m≤M,n≤N。
5.根据权利要求3所述的气体泄漏位置三维定位方法，其特征在于，所述处理获取的声场信息具体包括：
对参考传感器的信息xr(t)和第i个扫描传感器的信息xi(t)进行傅里叶变换，得到如下表达式：
Figure FDA0003238664420000021
其中，
Figure FDA0003238664420000022
和
Figure FDA0003238664420000023
分别表示xr(t)和xi(t)的傅里叶变换，e-j2πft为复变函数；
进一步得到如下表达式：
Figure FDA0003238664420000024
xi(t)和xr(t)的互功率谱表达为：
Figure FDA0003238664420000025
其中，
Figure FDA0003238664420000026
是xr(t)的共轭谱，
Figure FDA0003238664420000027
为复变函数；
第i个扫描传感器采集的信号与参考传感器同时采集的参考信号之间的相位差
Figure FDA0003238664420000028
表达为：
Figure FDA0003238664420000029
其中，Im[R(f)]表示R(f)的虚部，Re[R(f)]表示R(f)的实部；
虚拟相控传感器阵列中第一个扫描传感器采集到的信号为X1(t)，则，第i个传感器采集到的信号Xi(t)表达为：
Figure FDA00032386644200000210
其中，
Figure FDA00032386644200000211
是X1(t)与同时接收的参考信号之间的相位差，
Figure FDA00032386644200000212
是Xi(t)与同时接收的参考信号之间的相位差，
Figure FDA00032386644200000213
是复变函数。
6.根据权利要求1所述的气体泄漏位置三维定位方法，其特征在于，所述确定泄漏源的三维位置信息包括：
利用声全息数据和声场模型进行泄漏源声场重建计算，将计算结果显示为泄漏源图像信息；
以设定重建距离间隔，在不同距离上多次进行声场重建计算，显示多个泄漏源图像信息；
确定清晰度最高的泄漏源图像信息对应的泄漏源位置信息为泄漏源三维位置信息。
7.根据权利要求6所述的气体泄漏位置三维定位方法，其特征在于，所述泄漏源图像信息包括二维图像信息或三维图像信息。
8.根据权利要求2所述的气体泄漏位置三维定位方法，其特征在于，所述预设阵列规则包括阵元间距、阵元数量、阵列形状。
9.根据权利要求1所述的气体泄漏位置三维定位方法，其特征在于，所述确定泄漏源的三维位置信息包括全息显示所述泄漏源位置。
10.气体泄漏位置三维定位系统，用于实现权利要求1-8任一项所述的气体泄漏位置三维定位方法，其特征在于，包括：
信息采集组件，包括至少一个扫描传感器和至少一个参考传感器，设置用于获取泄漏源在测量平面的声场信息；
线性移动平台，用于控制所述扫描传感器和所述参考传感器形成虚拟相控阵列；
数据处理组件，其中嵌入有声场模型，设置用于传输并处理所述声场信息，获取声全息数据，利用建立的声场模型和处理得到的声全息数据，计算确定泄漏源的三维位置信息。
CN202111011720.4A 2021-08-31 2021-08-31 气体泄漏位置三维定位方法、定位系统 Active CN113639934B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111011720.4A CN113639934B (zh) 2021-08-31 2021-08-31 气体泄漏位置三维定位方法、定位系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111011720.4A CN113639934B (zh) 2021-08-31 2021-08-31 气体泄漏位置三维定位方法、定位系统

Publications (2)

Publication Number Publication Date
CN113639934A true CN113639934A (zh) 2021-11-12
CN113639934B CN113639934B (zh) 2024-03-26

Family

ID=78424594

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111011720.4A Active CN113639934B (zh) 2021-08-31 2021-08-31 气体泄漏位置三维定位方法、定位系统

Country Status (1)

Country Link
CN (1) CN113639934B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117629534A (zh) * 2023-11-28 2024-03-01 北京航天试验技术研究所 气体泄漏检测方法、装置及气体泄漏样本采集系统

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0042945A1 (en) * 1980-06-25 1982-01-06 International Business Machines Corporation Apparatus for and method of storing characteristics of an acoustic field
US7739899B1 (en) * 2005-04-29 2010-06-22 Iowa State University Research Foundation, Inc. Leak detection using structure-borne noise
CN102445307A (zh) * 2011-07-20 2012-05-09 南京帝坝工程科技有限公司 单井地下水流速流向及水库渗漏点测量方法及其测量装置
CN102790931A (zh) * 2011-05-20 2012-11-21 中国科学院声学研究所 一种三维声场合成中的距离感合成方法
WO2013136472A1 (ja) * 2012-03-14 2013-09-19 中国電力株式会社 チューブリーク検知装置、及びチューブリーク検知方法
US20160091387A1 (en) * 2014-09-30 2016-03-31 Uchicago Argonne, Llc Acoustic building infiltration measurement system
US20170016797A1 (en) * 2015-07-15 2017-01-19 Electronics And Telecommunications Research Institute Apparatus and method of monitoring gas based on variation in sound field spectrum
CN107255225A (zh) * 2017-05-18 2017-10-17 哈尔滨理工大学 基于加权修正参数p范数的管道泄漏高精度声学定位方法
CN108444653A (zh) * 2018-03-09 2018-08-24 北京理工大学 一种面阵型超声波气体泄漏点定位检测装置及方法
CN109813501A (zh) * 2019-01-30 2019-05-28 郑州大学 测量气体管道泄漏位置的方法、装置及系统
CN109885945A (zh) * 2019-02-26 2019-06-14 哈尔滨工程大学 一种半空间环境下的边界元法近场声全息变换方法
CN110045333A (zh) * 2019-04-12 2019-07-23 上海工程技术大学 一种基于卡尔曼滤波的声源三维定位方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0042945A1 (en) * 1980-06-25 1982-01-06 International Business Machines Corporation Apparatus for and method of storing characteristics of an acoustic field
US7739899B1 (en) * 2005-04-29 2010-06-22 Iowa State University Research Foundation, Inc. Leak detection using structure-borne noise
CN102790931A (zh) * 2011-05-20 2012-11-21 中国科学院声学研究所 一种三维声场合成中的距离感合成方法
CN102445307A (zh) * 2011-07-20 2012-05-09 南京帝坝工程科技有限公司 单井地下水流速流向及水库渗漏点测量方法及其测量装置
WO2013136472A1 (ja) * 2012-03-14 2013-09-19 中国電力株式会社 チューブリーク検知装置、及びチューブリーク検知方法
US20160091387A1 (en) * 2014-09-30 2016-03-31 Uchicago Argonne, Llc Acoustic building infiltration measurement system
US20170016797A1 (en) * 2015-07-15 2017-01-19 Electronics And Telecommunications Research Institute Apparatus and method of monitoring gas based on variation in sound field spectrum
CN107255225A (zh) * 2017-05-18 2017-10-17 哈尔滨理工大学 基于加权修正参数p范数的管道泄漏高精度声学定位方法
CN108444653A (zh) * 2018-03-09 2018-08-24 北京理工大学 一种面阵型超声波气体泄漏点定位检测装置及方法
CN109813501A (zh) * 2019-01-30 2019-05-28 郑州大学 测量气体管道泄漏位置的方法、装置及系统
CN109885945A (zh) * 2019-02-26 2019-06-14 哈尔滨工程大学 一种半空间环境下的边界元法近场声全息变换方法
CN110045333A (zh) * 2019-04-12 2019-07-23 上海工程技术大学 一种基于卡尔曼滤波的声源三维定位方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
RUJIA WANG等: "Optimization of Fixed Microphone Array in High Speed Train Noises Identification Based on Far-Field Acoustic Holography", 《ADVANCES IN ACOUSTICS AND VIBRATION》, pages 1894918 *
李磊 等: "基于超声波虚拟相控阵列的气体泄漏成像方法", 传感技术学报, vol. 32, no. 5, pages 677 - 680 *
李磊;刘庆辉;杨宽;田勇志;张斌;牧凯军;: "基于超声波虚拟相控阵列的气体泄漏成像方法", 传感技术学报, vol. 32, no. 5, 15 May 2019 (2019-05-15), pages 677 - 680 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117629534A (zh) * 2023-11-28 2024-03-01 北京航天试验技术研究所 气体泄漏检测方法、装置及气体泄漏样本采集系统

Also Published As

Publication number Publication date
CN113639934B (zh) 2024-03-26

Similar Documents

Publication Publication Date Title
Zuo et al. Damage identification for plate-like structures using ultrasonic guided wave based on improved MUSIC method
JP5496192B2 (ja) 音響場を再構成するための方法
CN109444861B (zh) 一种平面声呐阵阻抗特性近场声全息法校准测量方法
Quaegebeur et al. Dispersion-based imaging for structural health monitoring using sparse and compact arrays
US7739899B1 (en) Leak detection using structure-borne noise
Quaegebeur et al. Correlation-based imaging technique using ultrasonic transmit–receive array for non-destructive evaluation
Yang et al. Functional delay and sum beamforming for three-dimensional acoustic source identification with solid spherical arrays
CN102243306A (zh) 舰船目标宽带雷达散射特性的水声测量方法
CN103033816A (zh) 基于圆弧扫描转换的合成孔径聚焦超声成像实现方法
US6958950B2 (en) Method of visualizing sound fields of individual sound sources using acoustic holography
CN113868583B (zh) 一种子阵波束聚焦的声源距离计算方法及系统
CN111812581B (zh) 基于原子范数的球面阵列声源波达方向估计方法
Sun et al. Array geometry calibration for underwater compact arrays
US8371172B2 (en) Method and system for predicting acoustic fields based on generalized moving frame acoustic holography
Zhao et al. Design of low-complexity 3-D underwater imaging system with sparse planar arrays
CN113639934B (zh) 气体泄漏位置三维定位方法、定位系统
Li et al. Three-dimensional localization of gas leakage using acoustic holography
Varray et al. Nonlinear radio frequency image simulation for harmonic imaging: Creanuis
CN105572230B (zh) 一种用于裂纹类缺陷定量识别的极性加权矢量全聚焦成像方法
Bao et al. A novel multi-site damage localization method based on near-field signal subspace fitting using uniform linear sensor array
Comesana et al. An introduction to virtual phased arrays for beamforming applications
CN103048641A (zh) 具有联合约束优化形式的矢量阵稳健聚焦处理方法
US5341809A (en) Ultrasonic flowmeter
CN116952356A (zh) 基于浅海环境水下声全息技术的近场辐射噪声测量方法
Dolmatov et al. Frequency-domain imaging algorithm for ultrasonic testing by application of matrix phased arrays

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