CN113639680A - 基于sh导波频散测量管道厚度的方法 - Google Patents
基于sh导波频散测量管道厚度的方法 Download PDFInfo
- Publication number
- CN113639680A CN113639680A CN202110849350.5A CN202110849350A CN113639680A CN 113639680 A CN113639680 A CN 113639680A CN 202110849350 A CN202110849350 A CN 202110849350A CN 113639680 A CN113639680 A CN 113639680A
- Authority
- CN
- China
- Prior art keywords
- wave
- frequency
- expression
- thickness
- domain
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B17/00—Measuring arrangements characterised by the use of infrasonic, sonic or ultrasonic vibrations
- G01B17/02—Measuring arrangements characterised by the use of infrasonic, sonic or ultrasonic vibrations for measuring thickness
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Length Measuring Devices Characterised By Use Of Acoustic Means (AREA)
Abstract
本公开提供了基于SH导波频散测量管道厚度的方法,包括对柱状单层弹性介质外壁上环形剪切源激发的SH波场进行求解;通过近似方法求得SH波场的频散曲线近似解析解;根据近似解析解推导给出管道厚度及横波速度反演公式;通过实轴积分法进行SH波场时域全波波列的数值模拟,得到全波波列数据;通过频率波数法利用全波波列数据获取频散曲线;将相邻两阶频散曲线上的所有连续数据点代入管道厚度及横波速度反演公式进行计算,得到反演结果。本公开估测结果更加准确;无需事先已知被测材料声速,降低测量要求,提高测量方法的适用性;对管道材料无要求,能广泛适用于各种材质的管道测量;可测量厚度范围较大,准确度较高。
Description
技术领域
本公开涉及管道厚度测量技术领域,尤其涉及基于SH导波频散测量管道厚度的方法。
背景技术
随着工业现代化进程的加快,无损检测技术越来越受到大家的关注。在众多设备、管道、器件的制造和维修时,必须测量其厚度以确认规格大小、腐蚀情况、磨损情况,以保证产品质量和生产安全。其中管道的腐蚀检测在工业生产和运输中尤为重要。在管道使用年限较长时,由于内外部环境的原因可能造成管道腐蚀,导致管道壁逐渐变薄,进入事故高发期。因此需要测量管道厚度,确定管道腐蚀情况,及时发现风险,保证管道安全。
根据测量原理的不同,常用的测厚方法主要分为:(1)超声波测厚。在被测材料表面激发超声波,传播过程中遇到另一界面会发生反射,通过测量间隔时间并与已知的被测物体的波速相结合计算被测物体厚度。(2)磁性测厚。当导磁材料的表面覆盖有非导磁薄层时,测量该材料磁阻,结果会随表面的非导磁薄层厚度变化而发生变化。利用此方法便可反演出该导磁材料的表面覆盖的非导磁薄层厚度。(3)涡流测厚。当金属表面含有覆盖层时,将高频通电线圈放置在其表面,会在金属中产生涡流,此涡流的磁场又反作用于高频通电线圈,改变其阻抗值。通过测量其阻抗的变化值可反演金属表面覆盖层厚度。(4)同位素测厚。材料对于辐射的吸收与散射会随着材料厚度变化而发生变化,通过测量对辐射的吸收与散射情况可反演薄板、薄层类材料厚度。
但以上每种测厚方法都有其局限性。磁性测厚仅适用于导磁材料;涡流测厚仅适用于金属材料;同位素测厚仅适用于较薄的被测材料。超声测厚适用性较广,但该方法通常需要事先已知被测材料声速,但在实际生产生活中,无法确保所有被测材料的声速全部已知且准确。部分弹性材料的声速受压力,温度影响较大,部分孔隙材料的声速受其孔隙度、弯曲度等孔隙参数影响较大,若无法已知被测物体的准确声速,使用传统超声测厚法会产生较大误差。
发明内容
本公开的目的是要提供一种基于SH导波频散测量管道厚度的方法,可以解决上述现有技术问题中的一个或者多个。
根据本公开的一个方面,提供一种基于SH导波频散测量管道厚度的方法,包括以下步骤,
对柱状单层弹性介质外壁上环形剪切源激发的SH波场进行求解;
通过近似方法求得SH波场的频散曲线近似解析解;
根据近似解析解推导给出管道厚度及横波速度反演公式;
通过实轴积分法进行SH波场时域全波波列的数值模拟,得到全波波列数据;
通过频率波数法利用全波波列数据获取频散曲线;
将相邻两阶频散曲线上的所有连续数据点代入管道厚度及横波速度反演公式进行计算,得到反演结果。
在可能的实施方式中,对柱状单层弹性介质外壁上环形剪切源激发的SH波场进行求解包括,
设χ为弹性介质的剪切横波势,则χ满足:
式中,cs为弹性介质的横波速度;
采用柱坐标系求解,在柱坐标系下,当r=r1时,内表面处切向应力分量为0;当r=r2时,外表面处的切向应力分量与环形剪切源的切向应力分量一致,r1和r2分别为柱状单层弹性介质的内径和外径;
根据环形剪切源振动方式,弹性介质表面的已知切向应力表达式如下,σrθ(r,z,t)=σrθ(z)·f(t) (2)
式中,σrθ(z)为切向应力分量,f(t)为声源脉冲的函数;
边界表达式如下,
r=r1,σrθ=0 (3)
r=r2,σrθ=σrθ(z)·f(t) (4)
利用分离变量法解方程(1),得到χ的Fourier积分解为,
式中,A1和B1是待求未知量,I0(z)和K0(z)分别是修正Bessel函数,kz是轴向波数,ω是角频率,F(ω)声源脉冲的频率域函数,m为径向虚波数,ks为横波波数。
在柱坐标系下,剪切横波势χ,切向位移uθ和切向应力分量σrθ的关系如下,
式中,μ是弹性介质的剪切模量;
根据边界表达式,联立以上方程式(3)-(8),得到矩阵方程如下,
式中,等号左端的2×2的矩阵的行列式记为D(kz,ω),则aij(i,j=1,2)式行列式元素表达式如下,
a11=-m2I2(mr2) (10)
a12=-m2K2(mr2) (11)
a21=-m2I2(mr1) (12)
a22=-m2K2(mr1) (13)
式中,I2(z)和K2(z)是修正Bessel函数。
在可能的实施方式中,通过近似方法求得SH波场的频散曲线近似解析解包括以下步骤,
方程(9)有非零解的充分必要条件是系数行列式为0,即
D(kz,ω)=0 (14)
得到环形剪切源激发的SH波长的频散方程如下,
对方程(15)进行r1→∞的近似求解,当I2(z)和K2(z)中的z→∞,修正Bessel函数Iv(z)和Kv(z)的渐近表达式如下,
将公式(16)带入频散方程(15),可得表达式如下,
对表达式(17)进行整理可得表达式如下,
式中,m为径向虚波数,其表达式如下,
设参数h为管道厚度和α为横波径向波数,h和α满足以下表达式,
将表达式(19)和表达式(20)带入表达式(18)得到以下方程,
eiαh-e-iαh=0 (21)
sin(αh)=0 (22)
αh=nπ,(n=0,1,2……) (23)
整理得到表达式如下,
式中,c为阶数为n的角频率为ω时的SH导波的相速度,cs为横波速度。
在可能的实施方式中,根据近似解析解推导给出管道厚度及横波速度反演公式包括以下步骤,
设频散曲线第n阶模式上的两点坐标为(f1,c1)和(f2,c2),且ω=2πf,则带入表达式(24)得到方程组如下,
解方程组(25)得到,
设频散曲线第n+1阶模式上的两点坐标为(f3,c3)和(f4,c4),且ω=2πf,则带入表达式(24)得到方程组如下,
解方程组(27)得到,
整理可得柱状单层弹性介质的厚度反演公式如下,
同时可得柱状单层弹性介质的横波速度反演公式如下,
在可能的实施方式中,通过实轴积分法进行SH波场时域全波波列的数值模拟,得到全波波列数据包括以下步骤,
将方程(9)等号右侧的列矩阵分别代替行列式D(kz,ω)中的第一列和第二列,对方程(9)进行求解,可得,
Ni(kz,ω)(i=1,2) (31)
在频率波数域,外表面处的剪切位移μθ(r2,z,ω)的表达式为,
则剪切位移的全波表达式为,
选用余弦包络脉冲作为声源脉冲,声源脉冲的时域函数表达式如下,
式中,f0是声源中心频率,T是声源脉冲宽度;
表达式(36)对应的声源脉冲的频率域函数为:
在待测柱状单层弹性介质***内模拟全波,即可获取SH波场时域全波波列。
在可能的实施方式中,通过频率波数法利用全波波列数据获取频散曲线包括以下步骤,
将全波列数据从时间空间域变换到频率空间域,得到SH波场时域全波每一道所对应的频谱;
对频率空间域波场的每一列作同样的傅里叶变换,即将SH波场信息从频率空间域变换到频率波数域,得到频率波数域的能量映射图;
在频率波数域中对所有数据点进行重新采样,将坐标系由频率波数域转换到频率速度域,得到频率速度域的能量映射图;
在频率速度域能量映射图中提取频散曲线。
在可能的实施方式中,将相邻两阶频散曲线上的所有连续数据点代入管道厚度及横波速度反演公式进行计算时,将所有连续数据点全部代入反演公式,将得到的结果取中值。
在可能的实施方式中,环形剪切源为环状水平剪切源。
本公开提供技术方案与现有技术相比,存在以下有益效果:利用多阶SH导波反演柱状单层介质厚度及横波速度,SH导波在测量过程中不发生波形转换,使估测结果更加准确;无需事先已知被测材料声速,降低测量要求,提高测量方法的适用性;对管道材料无要求,能广泛适用于金属材料、非金属材料、弹性介质和孔隙介质的管道测量;可测量厚度范围较大,准确度较高。
另外,在本公开技术方案中,凡未作特别说明的,均可通过采用本领域中的常规手段来实现本技术方案。
附图说明
为了更清楚地说明本公开实施例的技术方案,下面将对实施例描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本公开的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本公开一实施例提供的基于SH导波频散测量管道厚度的方法的流程图。
图2为本公开一实施例提供的柱状单层弹性介质理论模型的示意图。
图3为声源脉冲的时域波形的波形图。
图4为本公开一实施例提供的多个接收器接收到的32道时域全波波列。
图5为本公开一实施例提供的32道SH波场时域全波波列每一道所对应的频谱。
图6为本公开一实施例提供的频率波数域的能量映射图。
图7为本公开一实施例提供的频率速度域的能量映射图。
图8为本公开一实施例提供的在频率速度域的能量映射图中提取频率速度点。
图9为本公开一实施例提供的选取频率速度数据点的曲线图。
具体实施方式
为使本公开实施例的目的、技术方案和优点更加清楚,下面将结合本公开实施例中的附图,对本公开实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本公开一部分实施例,而不是全部的实施例。基于本公开中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本公开保护的范围。
实施例:
如图1所示为本公开一实施例提供的基于SH导波频散测量管道厚度的方法,包括以下步骤:
步骤S11:对柱状单层弹性介质外壁上环形剪切源激发的SH波场进行求解。
具体的,环形剪切源可以采用环形水平剪切源,使用环形水平剪切源激发SH波的方法,对柱状单层弹性介质外壁上环形剪切源激发的SH波场进行求解。
参考说明书附图2,示出了柱状单层弹性介质理论模型,假设柱状单层弹性介质内外均为未知流体,环形剪切源用圆环模拟示意,并假设圆环表面紧贴柱状单层弹性介质外壁且与外壁耦合,有效激发SH波并传递到柱状单层弹性介质内。
在可选的实施例中,环形剪切源可以是声源换能器,由此,声源换能器的声源振幅与波形函数已知,可激发出环形剪切振动,环形剪切振动带动柱状单层弹性介质表面按已知规律沿周向振动。
在可选的实施例中,环形剪切源与柱状单层弹性介质均为轴对称模型且居中对称,r1和r2分别为柱状单层弹性介质的内径和外径。由此,能够保证模型整体具有轴对称性。环形剪切源在柱状单层弹性介质中单一激发SH波,无P波和SV波的产生,便于计算与测量。
具体的,对柱状单层弹性介质外壁上环形剪切源激发的SH波场进行求解包括以下步骤:
设χ为弹性介质的剪切横波势,则χ满足:
式中,cs为弹性介质的横波速度;
由于模型整体具有轴对称性,因此可以采用柱坐标系求解,在柱坐标系下,当r=r1时,内表面处切向应力分量为0;当r=r2时,外表面处的切向应力分量与环形剪切源的切向应力分量一致,r1和r2分别为柱状单层弹性介质的内径和外径;
根据环形剪切源振动方式,弹性介质表面的已知切向应力表达式如下,σrθ(r,z,t)=σrθ(z)·f(t) (2)
式中,σrθ(z)为切向应力分量,f(t)为声源函数;
边界表达式如下,
r=r1,σrθ=0 (3)
r=r2,σrθ=σrθ(z)·f(t) (4)
利用分离变量法解方程(1),得到χ的Fourier积分解为,
式中,A1和B1是待求未知量,I0(z)和K0(z)分别是修正Bessel函数,kz是轴向波数,ω是角频率,F(ω)是声源脉冲的频率域函数,m为径向虚波数,ks为横波波数。
在柱坐标系下,剪切横波势χ,切向位移uθ和切向应力分量σrθ的关系如下,
式中,μ是弹性介质的剪切模量;
根据边界表达式,联立以上方程式(3)-(8),得到矩阵方程如下,
式中,等号左端的2×2的矩阵的行列式记为D(kz,ω),则aij(i,j=1,2)式行列式元素表达式如下,
a11=-m2I2(mr2) (10)
a12=-m2K2(mr2) (11)
a21=-m2I2(mr1) (12)
a22=-m2K2(mr1) (13)
式中,I2(z)和K2(z)是修正Bessel函数。
步骤S12:通过近似方法求得SH波场的频散曲线近似解析解。
具体的,通过近似方法求得SH波场的频散曲线近似解析解包括以下步骤,
方程(9)有非零解的充分必要条件是系数行列式为0,即
D(kz,ω)=0 (14)
求解方程(14),得到的根为特定频率ω下的轴向波数kz,对应柱状单层弹性介质中沿轴向传播的导波,得到环形剪切源激发的SH波长的频散方程如下,
对方程(15)进行r1→∞的近似求解,当I2(z)和K2(z)中的z→∞,修正Bessel函数Iv(z)和Kv(z)的渐近表达式如下,
将公式(16)带入频散方程(15),可得表达式如下,
对表达式(17)进行整理可得表达式如下,
式中,m为径向虚波数,其表达式如下,
设参数h为管道厚度,参数α为横波径向波数,h和α满足以下表达式,
将表达式(19)和表达式(20)带入表达式(18)得到以下方程,
eiαh-e-iαh=0 (21)
sin(ah)=0 (22)
αh=nπ,(n=0,1,2……) (23)
整理得到表达式如下,
式中,c为阶数为n的角频率为ω时的SH导波的相速度,cs为横波速度。
步骤S13:根据近似解析解推导给出管道厚度及横波速度反演公式。
具体的,根据近似解析解推导给出管道厚度及横波速度反演公式包括以下步骤,
设频散曲线第n阶模式上的两点坐标为(f1,c1)和(f2,c2),且ω=2πf,则带入表达式(24)得到方程组如下,
解方程组(25)得到,
设频散曲线第n+1阶模式上的两点坐标为(f3,c3)和(f4,c4),且ω=2πf,则带入表达式(24)得到方程组如下,
解方程组(27)得到,
整理可得柱状单层弹性介质的厚度反演公式如下,
由此,只需已知某一阶频散曲线上的两个频率速度点即可求出该阶频散曲线的n,若求出相邻两阶分别对应的n与n+1,则可反演出柱状单层弹性介质的厚度。
同时可得柱状单层弹性介质的横波速度反演公式如下,
由此,只需已知某一阶频散曲线上的两个频率速度点即可求出柱状单层弹性介质的横波速度。
步骤S14:通过实轴积分法进行SH波场时域全波波列的数值模拟,得到全波波列数据。
具体的,通过实轴积分法进行SH波场时域全波波列的数值模拟,得到全波波列数据包括以下步骤,
将方程(9)等号右侧的列矩阵分别代替行列式D(kz,ω)中的第一列和第二列,对方程(9)进行求解,可得,
Ni(kz,ω)(i=1,2) (31)
将全波声场接收器紧贴柱状单层弹性介质的外表面,在声源函数已知的情况下,可以通过实轴积分法计算外表面处的剪切位移μθ(t)的瞬态波形。在频率波数域,外表面处的剪切位移μθ(r2,k,ω)的表达式为,
则剪切位移的全波表达式为,
选用余弦包络脉冲作为声源脉冲,声源脉冲的时域函数表达式如下,
式中,f0是声源中心频率,T是声源脉冲宽度;
表达式(36)对应的声源脉冲的频率域函数为,
在待测柱状单层弹性介质***内模拟全波,即可获取SH波场时域全波波列的数值。
在可选的实施例中,采用声源中心频率f0=100kHz,声源脉冲宽度为T=1×10-4s为声源参数的声源脉冲,在柱状单层弹性介质内模拟全波。声源脉冲的时域波形参考说明书附图3。
在本实施例中,选用32道全波波列进行模拟,在可选的实施例中,也可选用8道、16道或64道等全波波列进行模拟。
在可选的实施例中,可以选取模型参数进行SH波场时域全波波列的数值模拟。具体模拟参数如下:
横波速度cs=3.2km/s;密度ρ=7.7g/cm3;内径r1=0.1m;厚度=r2-r1=0.05m。
此时,32道SH波场时域全波波列数据的模拟结果参考说明书附图4,说明书附图4展示了多个接收器接收到的32道时域全波波列,第1道接收器源距为0.1m,之后n道依次增加0.03m,即第一道源距0.1m,第二道源距0.13m,第三道源距0.16m,以此类推。
步骤S15:通过频率波数法利用全波波列数据获取频散曲线。
具体的,通过频率波数法利用全波波列数据获取频散曲线包括以下步骤,
将全波列数据从时间空间域变换到频率空间域,得到32道SH波场时域全波波列每一道所对应的频谱,参考说明书附图5所示,其中,X轴为频率,Y轴为道数,Z轴为频谱所对应的能量大小,由说明书附图5可知,频谱能量随着道数的增加(即源距变大)而逐渐减小,且每一道能量均集中在100kHz附近,与环形剪切源中心频率f0=100kHz相符;
对频率空间域波场的每一列作相同的傅里叶变换,将SH波场信息从频率空间域变换到频率波数域,得到频率波数域的能量映射图,参考说明书附图6,由说明书附图6可知,能量集中在声源中心频率附近,且在激发频段内有多阶SH波能量;
在频率波数域中对所有数据点进行重新采样,将坐标系由频率波数域转换到频率速度域,得到频率速度域的能量映射图,参考说明书附图7,由说明书附图7可知,能量集中在声源中心频率附近,且集中在四条曲线上;
在频率速度域能量映射图中提取频散曲线,由于在在频率波数域中,波场的能量最大值在D(kz,ω)=0附近,这正是环形剪切源激发的柱状单层弹性介质中SH波频散方程,在频率速度域能量映射图中能量越大越集中的位置越接近频散曲线的准确值。因此在频率速度域能量映射图中针对横轴的每一个固定频率点,寻找该频率下声场最大值所对应的速度值,对每一个固定频率点进行相同操作,提取频散曲线的结果参考说明书附图8所示。
在可选的实施例中,在对频率空间域波场的每一列作相同的傅里叶变换时,由于傅里叶变换的积分为无穷,对于从时间空间域变换到频率空间域,需要采集足够多的点直至时域波形的幅度趋近于0,对于从频率空间域变换到频率波数域,需要采集足够多道数的数据直至信号幅度由于衰减趋近于0。但在实际反演过程中,可能无法采集足够的点使数据趋近于0,因此本公开中的傅里叶变换均进行补零措施,采用2048点的傅里叶变换。
在可选的实施例中,在频率波数域中对所有数据点进行重新采样时,利用公式v=f/k,其中,v表示相速度,f表示频率,k表示波数。
步骤S16:将相邻两阶频散曲线上的所有连续数据点代入管道厚度及横波速度反演公式进行计算,得到反演结果。
具体的,在选取数据点时,需要以能量映射图为参考,在激发频段内能量集中的曲线上取点,以保证选取的频率速度点的准确性。具体的,可以选取较为完整的相邻两阶频散曲线上的所有的连续数据点。
在可选的实施例中,参考说明书附图9,将相邻两阶频散曲线上的所有连续数据点代入管道厚度及横波速度反演公式进行计算时,将说明书附图9中所示单个选取框内的所有连续数据点全部代入反演公式,将得到的结果取中值。由此,中值与平均值相比更能反应数据组的一般水平,在可能存在极端值的数据组中,中值更加准确。
在具体模拟参数为:横波速度cs=3.2km/s;密度ρ=7.7g/cm3;内径r1=0.1m;厚度h=r2-r1=0.05m时,得到管道厚度与横波速度的计算结果如表1与表2所示。
表1管道厚度计算结果
速度(km/s) | n阶反演 | 误差 | n+1阶反演 | 误差 |
3.2 | 3.1897 | 0.3218% | 3.2003 | 0.0094% |
表2横波速度计算结果
通过表1与表2中的计算结果,可以观察到n阶数据所得到的结果误差均大于n+1阶数据所得到的结果。这是由于解析解是通过近似得到,在频率较低时存在误差,随着频率的增大而误差减小。由于n阶数据与n+1阶数据相比,频率较小,因此n阶数据所得到的结果误差均大于n+1阶数据所得到的结果误差。
在可选的实施例中,可以对以上初步反演结果进行优化,优化方法步骤如下:
根据初步反演,表1中可得到和h。由于初步反演误差结果不大,且阶数n为整数,因此将三者结合结果四舍五入可求出n与n+1的值,即所取两阶频散曲线的阶数。由于n阶的数据点误差大于n+1阶,因此舍去n阶的数据,利用n+1和的值,求得最终反演结果h。
本公开中基于SH导波频散测量管道厚度的方法,在初步反演时使用多阶是因为无法确定所提取到的频散曲线为哪一阶,结合相邻两阶可以避免这一问题。而已得到初步反演结果且阶数为整数,即可进行反推,确定较高一阶的阶数。继而便可单独使用较高一阶的数据点进行反演,以减小误差。
优化后的计算结果如表3所示:
厚度(m) | 反演 | 误差 | 速度(km/s) | 反演 | 误差 |
0.05 | 0.0500 | <0.0001% | 3.2 | 3.2003 | 0.0094% |
表3优化后的管道厚度与横波速度计算结果
表3中的计算结果展示了柱状单层弹性介质模型下的管道厚度与横波速度最终反演的结果。从表3中可以得到,最终反演的结果,准确度较高。该结果也证实了本公开提供的基于SH导波频散测量管道厚度的方法的正确性。
本公开提供的基于SH导波频散测量管道厚度的方法,利用多阶SH导波反演柱状单层介质厚度及横波速度,SH导波在测量过程中不发生波形转换,使估测结果更加准确;无需事先已知被测材料声速,降低测量要求,提高测量方法的适用性;对管道材料无要求,能广泛适用于金属材料、非金属材料、弹性介质和孔隙介质的管道测量;可测量厚度范围较大,准确度较高。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分的方法。
最后应说明的是:以上实施例仅用以说明本公开的技术方案,而非对其限制;尽管参照前述实施例对本公开进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本公开各实施例技术方案的精神和范围。
Claims (8)
1.基于SH导波频散测量管道厚度的方法,其特征在于,包括以下步骤:
对柱状单层弹性介质外壁上环形剪切源激发的SH波场进行求解;
通过近似方法求得SH波场的频散曲线近似解析解;
根据近似解析解推导给出管道厚度及横波速度反演公式;
通过实轴积分法进行SH波场时域全波波列的数值模拟,得到全波波列数据;
通过频率波数法利用全波波列数据获取频散曲线;
将相邻两阶频散曲线上的所有连续数据点代入管道厚度及横波速度反演公式进行计算,得到反演结果。
2.根据权利要求1所述的基于SH导波频散测量管道厚度的方法,其特征在于,所述对柱状单层弹性介质外壁上环形剪切源激发的SH波场进行求解包括,
设χ为弹性介质的剪切横波势,则χ满足:
式中,cs为弹性介质的横波速度;
采用柱坐标系求解,在柱坐标系下,当r=r1时,内表面处切向应力分量为0;当r=r2时,外表面处的切向应力分量与环形剪切源的切向应力分量一致,r1和r2分别为柱状单层弹性介质的内径和外径;
根据环形剪切源振动方式,弹性介质表面的已知切向应力表达式如下,
σrθ(r,z,t)=σrθ(z)·f(t) (2)
式中,σrθ(z)为切向应力分量,f(t)为声源函数;
边界表达式如下,
r=r1,σrθ=0 (3)
r=r2,σrθ=σrθ(z)·f(t) (4)
利用分离变量法解方程(1),得到χ的Fourier积分解为,
式中,A1和B1是待求未知量,I0(z)和K0(z)分别是修正Bessel函数,kz是轴向波数,ω是角频率,F(ω)是声源脉冲的频率域函数,m为径向虚波数,ks为横波波数;
在柱坐标系下,剪切横波势χ,切向位移uθ和切向应力分量σrθ的关系如下,
式中,μ是弹性介质的剪切模量;
根据边界表达式,联立以上方程式(3)-(8),得到矩阵方程如下,
式中,等号左端的2×2的矩阵的行列式记为D(kz,ω),则aij(i,j=1,2)式行列式元素表达式如下,
a11=-m2I2(mr2) (10)
a12=-m2K2(mr2) (11)
a21=-m2I2(mr1) (12)
a22=-m2K2(mr1) (13)
式中,I2(z)和K2(z)是修正Bessel函数。
3.根据权利要求2所述的基于SH导波频散测量管道厚度的方法,其特征在于,所述通过内径趋于无穷的近似方法求得SH波场的频散曲线近似解析解包括以下步骤,
方程(9)有非零解的充分必要条件是系数行列式为0,即
D(kz,ω)=0 (14)
得到环形剪切源激发的SH波长的频散方程如下,
对方程(15)进行r1→∞的近似求解,当I2(z)和K2(z)中的z→∞,修正Bessel函数Iv(z)和Kv(z)的渐近表达式如下,
将公式(16)带入频散方程(15),可得表达式如下,
对表达式(17)进行整理可得表达式如下,
式中,m为径向虚波数,其表达式如下,
设参数h为管道厚度,参数α为横波径向波数,h和α满足以下表达式,
将表达式(19)和表达式(20)带入表达式(18)得到以下方程,
eiαh-e-iαh=0 (21)
sin(αh)=0 (22)
αh=nπ,(n=0,1,2......) (23)
整理得到表达式如下,
式中,c为阶数为n的角频率为ω时的SH导波的相速度,cs为横波速度。
5.根据权利要求4所述的基于SH导波频散测量管道厚度的方法,其特征在于,所述通过实轴积分法进行SH波场时域全波波列的数值模拟,得到全波波列数据包括以下步骤,
将方程(9)等号右侧的列矩阵分别代替行列式D(kz,ω)中的第一列和第二列,对方程(9)进行求解,可得,
Ni(kz,ω)(i=1,2) (31)
在频率波数域,外表面处的剪切位移μθ(r2,z,ω)的表达式为,
则剪切位移的全波表达式为,
选用余弦包络脉冲作为声源脉冲,声源脉冲的时域函数表达式如下,
式中,f0是声源中心频率,T是声源脉冲宽度;
表达式(36)对应的声源脉冲的频率域函数为:
在待测柱状单层弹性介质***内模拟全波,即可获取SH波场时域全波波列。
6.根据权利要求5所述的基于SH导波频散测量管道厚度的方法,其特征在于,所述通过频率波数法利用全波波列数据获取频散曲线包括以下步骤,
将全波列数据从时间空间域变换到频率空间域,得到SH波场时域全波每一道所对应的频谱;
对频率空间域波场的每一列作同样的傅里叶变换,即将SH波场信息从频率空间域变换到频率波数域,得到频率波数域的能量映射图;
在频率波数域中对所有数据点进行重新采样,将坐标系由频率波数域转换到频率速度域,得到频率速度域的能量映射图;
在频率速度域能量映射图中提取频散曲线。
7.根据权利要求6所述的基于SH导波频散测量管道厚度的方法,其特征在于,所述将相邻两阶频散曲线上的所有连续数据点代入管道厚度及横波速度反演公式进行计算时,将所有连续数据点全部代入反演公式,将得到的结果取中值。
8.根据权利要求1所述的基于SH导波频散测量管道厚度的方法,其特征在于,所述环形剪切源为环状水平剪切源。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110849350.5A CN113639680B (zh) | 2021-07-27 | 2021-07-27 | 基于sh导波频散测量管道厚度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110849350.5A CN113639680B (zh) | 2021-07-27 | 2021-07-27 | 基于sh导波频散测量管道厚度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113639680A true CN113639680A (zh) | 2021-11-12 |
CN113639680B CN113639680B (zh) | 2022-08-16 |
Family
ID=78418507
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110849350.5A Active CN113639680B (zh) | 2021-07-27 | 2021-07-27 | 基于sh导波频散测量管道厚度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113639680B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2097583A1 (en) * | 1992-06-01 | 1993-12-02 | Paul J. Latimer | Hydrogen Damage Confirmation with Emats |
CN101819032A (zh) * | 2010-04-23 | 2010-09-01 | 北京工业大学 | 基于sh波的工业锅炉水垢厚度检测***及方法 |
CN102980942A (zh) * | 2012-11-20 | 2013-03-20 | 中国石油天然气集团公司 | 一种金属管道检测装置及方法 |
CN106018552A (zh) * | 2016-05-09 | 2016-10-12 | 南京航空航天大学 | 一种用于平板减薄缺陷定量化检测的无参考sh导波方法 |
CN108490079A (zh) * | 2018-03-19 | 2018-09-04 | 哈尔滨工业大学 | 一种基于超声换能器的波束成形方法 |
US20200150304A1 (en) * | 2018-11-08 | 2020-05-14 | Baker Hughes, A Ge Company, Llc | Casing wall thickness detection from higher order shear-horizontal mode signals |
US20200173773A1 (en) * | 2016-08-11 | 2020-06-04 | Guided Ultrasonics Ltd | Determining a thickness of a region of wall- or plate-like structure |
CN112903157A (zh) * | 2021-01-19 | 2021-06-04 | 吉林大学 | 基于纵向模态超声导波的圆管型结构的应力监测方法 |
-
2021
- 2021-07-27 CN CN202110849350.5A patent/CN113639680B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2097583A1 (en) * | 1992-06-01 | 1993-12-02 | Paul J. Latimer | Hydrogen Damage Confirmation with Emats |
CN101819032A (zh) * | 2010-04-23 | 2010-09-01 | 北京工业大学 | 基于sh波的工业锅炉水垢厚度检测***及方法 |
CN102980942A (zh) * | 2012-11-20 | 2013-03-20 | 中国石油天然气集团公司 | 一种金属管道检测装置及方法 |
CN106018552A (zh) * | 2016-05-09 | 2016-10-12 | 南京航空航天大学 | 一种用于平板减薄缺陷定量化检测的无参考sh导波方法 |
US20200173773A1 (en) * | 2016-08-11 | 2020-06-04 | Guided Ultrasonics Ltd | Determining a thickness of a region of wall- or plate-like structure |
CN108490079A (zh) * | 2018-03-19 | 2018-09-04 | 哈尔滨工业大学 | 一种基于超声换能器的波束成形方法 |
US20200150304A1 (en) * | 2018-11-08 | 2020-05-14 | Baker Hughes, A Ge Company, Llc | Casing wall thickness detection from higher order shear-horizontal mode signals |
CN112903157A (zh) * | 2021-01-19 | 2021-06-04 | 吉林大学 | 基于纵向模态超声导波的圆管型结构的应力监测方法 |
Non-Patent Citations (1)
Title |
---|
刘红梅: ""柱状多层Biot双相孔隙介质中的SH波"", 《中国优秀博硕士学位论文全文数据库(硕士)基础科学辑》 * |
Also Published As
Publication number | Publication date |
---|---|
CN113639680B (zh) | 2022-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103238065B (zh) | 用超声波进行管线壁性能测试的***和方法 | |
US10704901B2 (en) | Measuring wall thickness loss for a structure | |
JP6438957B2 (ja) | 欠陥を監視するためのシステム及び方法 | |
Plona et al. | Axisymmetric wave propagation in fluid‐loaded cylindrical shells. II: Theory versus experiment | |
CN1515876A (zh) | 分析双相流体流动的方法和*** | |
CN110702785B (zh) | 频散Lamb波多项式时频域模态分解和缺陷定位方法及装置 | |
CN110469782B (zh) | 一种基于自适应滤波器的管道泄漏定位装置 | |
CN105181793B (zh) | 基于超声波双频信号测量两相流含气率的方法 | |
Benstock et al. | The influence of surface roughness on ultrasonic thickness measurements | |
KR102204747B1 (ko) | 탄성파의 신호 에너지를 이용한 손상 탐지 방법 | |
CN110346004B (zh) | 一种双声道超声时差法的流量测量数据融合方法 | |
CN108469298B (zh) | 一种驻波管矢量水听器校准低频修正方法 | |
JP2002520582A (ja) | 音響型流量計における信号の交差計測 | |
CN111189930B (zh) | 基于管体的裂纹检测方法、***、设备及存储介质 | |
Gu et al. | The effect of internal pipe wall roughness on the accuracy of clamp-on ultrasonic flowmeters | |
CN113316705A (zh) | 用于确定细长或延伸结构的厚度的方法和*** | |
CN111350948A (zh) | 一种基于波束形成的管道泄漏位置计算方法 | |
Sun et al. | Technology of locating loose particles inside sealed electronic equipment based on Parameter-Optimized Random Forest | |
JP2017161406A (ja) | 伝搬定数の取得方法、吸音材の吸音率の計算方法、吸音率評価装置 | |
CN104748704A (zh) | 薄壁结构超声共振测厚频谱分析内插校正方法 | |
Dobson et al. | The scattering of torsional guided waves from Gaussian rough surfaces in pipework | |
CN114778690A (zh) | 一种增材制件气孔缺陷的激光超声定量检测方法 | |
CN113639680B (zh) | 基于sh导波频散测量管道厚度的方法 | |
CN105115552B (zh) | 基于侧音测相的超声波流量测量方法及装置 | |
CN113108853A (zh) | 一种提高低流速流体流量测量精度的方法 |
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 |